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^': 

^ Aims. Our purpose is to provide reliable stellar parameters for a significant sample of eclipsing binaries, which are representative of a whole 

}--i dwarf and metal-poor galaxy. We also aim at providing a new estimate of the mean distance to the SMC and of its depth along the line of sight 

C/2 , for the observed field of view. 

^ ^^ Methods. We use radial velocity curves obtained with the ESO FLAMES facility at the VLT and light curves from the OGLE-II photometric 
survey. The radial velocities were obtained by least-squares fits of the observed spectra to synthetic ones, excluding the hydrogen Balmer lines. 

y—( Results. Our sample contains 23 detached, 9 semi-detached and 1 overcontact systems. Most detached systems have properties consistent with 

^ stellar evolution calculations from single-star models at the standard SMC metallicity Z = 0.004, though they tend to be slightly overluminous. 

^^ The few exceptions are probably due to third light contribution or insufficient signal-to-noise ratio. The mass ratios are consistent with a flat 

CO ■ distribution, both for detached and semi-detached/contact binaries. A mass-luminosity relation valid from ~4 to ~18 Mq is derived. The 

C^ ■ uncertainties are in the ±2 to ±1 1 % range for the masses, in the ±2 to ±5 % range for the radii and in the ±1 to ±6 % range for the effective 

\l^ temperatures. The average distance modulus is 19. 1 1 ± 0.03 (66.4 ± 0.9 kpc). The moduli derived from the V and from the / data are consistent 

l/^ within 0.01 mag. The 2 a depth of the SMC is, for our field, of 0.25 mag or 7.6 kpc under the assumption of a gaussian distribution of stars 

f^ along the line of sight. Three systems show significant apsidai motion, one of them with an apsidai period of 7.6 years, the shortest known to 

^^ date for a detached system with main sequence stars. 

>: 

• "^ Key words, stars: early type - stars: binaries: eclipsing - stars: binaries: spectroscopic - stars: fundamental parameters - galaxies: Magellanic 
|r%] Clouds - distance scale 

1 . Introduction (2) precise distance moduli (DMi) derived from the luminosi- 

ties calculated from the combination of the absolute radii with 
Since the late 1990s, the usefulness of extragalactic eclips- ^^^ ^^^^^^^ temperatures. Until a purely geometrical distance 
ing binaries has been emphasized in a number of papers. The determination is feasible, Paczyriski (2001) considers that de- 
reader can notably r efer to the ex cellent reviews from Clausen ^^^^^^ gg^ ^^.^ ^j^^ ^^^^ promising distance indicators to the 
(12003 and Guinan ^^. The two major contributions Magellanic Clouds. Besides, Wyithe & Wilson ^ here- 
of eclipsing binaries (hereafter EBs) to astrophysics are to pro- ^^^^^ ^^q^) remarked that semi-detached EBs are even more 
vide (1) fundamental mass and radius measurements for the promising, since their parameters are better constrained, 
component stars, allowing to test stellar evolution models, and 

The renewal of interest in extragalactic EBs, especially EBs 

Send offprint requests to: -p. North i" *e Magellanic Clouds, has been stimulated by the release of 
* Based on observations made with the FLAMES-GIRAFFE multi- ^ huge number of light curves as a byproduct of automated mi- 
object spectrograph mounted on the Kuyen VLT telescope at ESO- crolensing surveys (EROS, MACHO, OGLE) with 1-m class 
Paranal Observatory (Swiss GTO programme 072.A-0474A; PI: P telescopes. As photometry is only half of the story, high reso- 
North) lution spectrographs attached to 4-m class or larger telescopes 
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had to be used to obtain reliable radial velocity (RV) curves. 
Four B-type EB systems belonging to the Large Magellanic 
Cloud (LMC) were accurately characterized in a series of pa- 
pers from Guinan et al. ( fT998] l. Ribas et al. (2000 2002) and 
Fitzpatrick et al. (120021 120031) . More recently, from high res- 
olution, high S /N spectra obtained with UVES at the ESO 
VLT, the analysis of eight more LMC systems was presented by 
Gonzalez et al. ( 120051 ). A corner stone in the study of EBs in the 
SMC was set up with the release of two papers from Harries et 
al. (120031 hereinafter HHH03) and Hilditch et al. (2005, here- 
inafter HHH05) giving the fundamental parameters of a total 
of 50 EB systems of spectral types O and B. The spectroscopic 
data were obtained with the 2dF multi-object spectrograph on 
the 3.9-m Anglo- Australian Telescope. To our knowledge this 
is the first use of multi-object spectroscopy in the field of extra- 
galactic EBs. Recently, even the distance to large spiral galax- 
ies were measured on the basis of two EBs in M31 (Ribas et 
al. 120051 Vilardell et al. l20T0l ) and one in M33 (Bonanos et al. 
120061 ). 

The huge asymmetry, between the number of EB light 
curves published so far and the very small number of RV 
curves, is striking. If one considers the SMC, the new OGLE-II 
catalog of EBs in the SMC (Wyrzykowski et al. 2004) con- 
tains 1350 light curves and currently only 50 of these systems 
have moderately reliable RV curves. This paper reduces a little 
this imbalance by releasing the analysis of 28 more EB sys- 
tems plus revised solutions for 5 systems previously described 
by HHH03 and HHH05. The RV measurements were derived 
from muti-object spectroscopic observations made with the 
VLT FLAMES facility. 

Another strong motivation for increasing the number of 
fully resolved binaries is to settle the problem of the distri- 
bution of the mass ratio of detached binaries with early B pri- 
maries. Recently, two papers were published supporting two di- 
ametrically opposed conclusions: van Rensbergen et al. (2006 ), 
whose work is based on the 9th Catalogue of Spectroscopic 
Binaries (Pourbaix et al. 120041 ). support the view that the q- 
distribution (where q is the mass ratio) follows a Salpeter-like 
decreasing power law; however, from the examination of the 
homogeneous sample of the 21 detached systems characterized 
by HHH03 and HHH05, Pinsonneault & Stanek (2006) draw 
the conclusion that the proportion of close detached systems 
with mass ratio q > 0.87 far outnumbers what can be expected 
from either a Salpeter or a flat (7-distribution (the "twins" hy- 
pothesis). Finally, let us mention that the ^'-distribution of semi- 
detached (i.e. evolved) systems is no more settled, the statistics 
strongly depending on the method used to find the mass ra- 
tios, i.e. from the light curve solution or from SB2 spectra (van 
Rensbergen et al. 2006). 

Although the controversy about the characteristic distance 
to the SMC seems to be solved in favour of a mid position 
between the "short" and the "long" scales, distance data and 
line of sight depth remain vital for comparison with theoreti- 
cal models concerning the three-dimensional structure and the 
kinematics of the SMC (Stanimirovic et al. 120041 ). 

Our contribution provides both qualitative and quantita- 
tive improvements over previous studies. Thanks to the VLT 
GIRAFFE facility, spectra were obtained with a resolution 



three times that in HHH03/05's study. Another strong point 
is the treatment of nebular emission. The SMC is known to 
be rich in H 11 regions (Fitzpatrick 119851 Torres & Carranza 
119871 ). Consequently, strong Balmer lines in emission are very 
often present in the spectra of the binary systems under study. 
Therefore, it appeared rapidly that it was essential to find a con- 
sistent way to deal with this "third component" polluting most 
double-lined (SB2) spectra. 

We present the observations in Section 2. The reduction of 
the spectroscopic data and the interpretation of both photomet- 
ric and spectroscopic data are described in Section 3, where the 
errors are also discussed in detail. The individual binary sys- 
tems are described in Section 4, while the sample as a whole is 
discussed in relation with the SMC properties in Section 5. 

2. Observations 

The targets, astrometry included, were selected from the first 
OGLE photometric catalog. The GIRAFFE field of view (FoV) 
constrained to choose systems inside a 25'-diameter circle. 
Other constraints were / < 18 mag, at least 15 well-behaved de- 
tached light curves (for the SMC field) and finally seven bump 
cepheids in the FoV (for another program). The positions in 
the sky of 33 objects studied in this paper are shown in Fig.[T] 
The epoch, exposure time, air mass, seeing and age of moon 
for each of the 16 CCD frames are gathered in Table[T] 

The relation between our own 1-33 labeling and the OGLE 
names can be found in Table |3] which lists the basic parame- 
ters of the systems. The coordinates are from Wyrzykowski et 
al. (120041 ). The orbital periods and epochs of the primary mini- 
mum are close to those listed by Wyrzykowski et al. (120041 ). but 
the periods were improved as far as possible using the radial 
velocity curves determined in this work. That was worth the 
eff'ort, since spectroscopic observations were performed more 
than three years after the last photometric ones. Since the times 
of the photometric minima are quite sharply defined, the uncer- 
tainty on the period (mentioned between parentheses in Table 
13) is based on the uncertainty of the spectroscopically defined 
epoch of the primary minimum. The latter is quite precise for 
circular orbits; for eccentric orbits, it is less accurate, because 
the Kepler equation had to be solved and the solution is af- 
fected by the uncertainty on the eccentricity. We have decided 
here to adopt the dynamic definition of the primary and sec- 
ondary components, rather than the photometric one. In other 
words, the primary component is always the more massive one, 
and the primary minimum always corresponds to the eclipse of 
the primary by the secondary component. As a consequence, 
it may happen that the so-called primary minimum is not the 
deeper one. Figure |2] shows the histogram of the periods. The 
strong observational bias in favour of short periods is conspic- 
uous. 

For all but two binaries, the light curves come from 
the new version of the OGLE-II catalog of ecUps- 
ing binaries detected in the SMC (Wyrzykowski et al. 
120041 ). This catalog is based on the Difference Image 
Analysis (DIA) catalog of variable stars in the SMC 
(see http : //sirius . astrouw . edu . pl/~ogle/ogle2/ 
smc_ecl/index.htinl). The data were coUected from 1997 
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Fig.l. Positions of the 33 binaries in the 1° x 1° surrounding field centered on (J2000) 00H9"20^ -73°13'37.3". They lie 
alongside the southwestern part of the SMC bar (oblique line) at ~0.5° of the SMC optical center. Image retrieved from the 
STScI Digitized Sky Survey (POSS-II/UK Schmidt Telescope, Red) (copyright ©1993-5 Anglo-Australian Observatory.) 



to 2000. The systems 4 121084 and 5 100731 were selected 
from the first version of the catalog (standard PSF photometry) 
but for an unknow reason they do no appear anymore in the 
new version. Nevertheless, they were retained as there is no 
objective reason to exclude them. 

The DIA photometry is based on /-band observations (be- 
tween 202 and 312 points per curve). B and V light curves were 
also used in spite of a much poorer sampling (between 22-28 
points/curve and 28-46 points/curve in B and V respectively). 
To give an idea of the accuracy of the OGLE photometry, the 
objects studied in this paper have an average / magnitude and 
scatter in the range 15.083 ±0.009 to 18.159±0.047. These val- 
ues were calculated from the best-fitting synthetic light curves. 
For the two other bands, we get 14.701 ±0.011 to 18.090±0.025 
for B and 14.966±0.009 to 18. 173±0.023 for y. The quality of 
an observed light curve can be better expressed by comparing 
the depth of the primary minimum A/mini to the average RMS 
scatter a-]. These ratios are shown in Table [161 This permits 
to classify the light curves in five categories: low (A/mini /ct/ < 
10), low-to-medium (10 < A/mini/cr/ < 20), medium (20 < 
A/mini /cT/ < 30), medium-to-high (30 < Al^nini /o-i < 40) and 
high (A/mini/cT/ > 40) quality. According to this scheme, most 
/-band light curves (58%) belong to the low-to-medium and 



medium quality categories, one-third (33%) in the medium-to- 
high and high quality categories, and the remaining 9% in the 
poor quality category. This classification scheme is not useful 
for the other bands, the low sampling being the limiting factor 

VLT FLAMES/GIRAFFE spectroscopy was obtained by 
one of us (FR) during eight consecutive nights from November 
16 to 23, 2003. The spectrograph was used in the low resolution 
(LR2) Medusa mode: resolving power R - 6400, bandwidth 
A/l = 603 A centered on 4272 A. The most prominent absorp- 
tion lines in the blue part of early-B stars spectra are: He, He i 
^4026, H(5, Hei ^4144, Hy, Hei ^4388, and Hei /14471. For 
late-O stars. He ii /14200 and He ii /14542 gain in importance. 
Two fields, one in the SMC and one in the LMC, were alter- 
natively observed at a rate of four exposures per night with an 
integration time of 2595 s for all but one epoch. Therefore, 16 
spectra per target were obtained, with a total of 104 targets in 
the SMC and 44 in the LMC. The LMC SB2 systems are being 
analyzed and will be presented in another paper 

Beside the spectra of the objects, 21 sky spectra were ob- 
tained for each exposure in the SMC. The parameters related 
to the spectroscopic observations are gathered in Table [T] The 
observed signal-to-noise ratios (S /N) were determined for each 
smoothed spectrum (see Section lTTI ) in the continuum between 
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4195 and 4240 A. For each object, two values are presented in 
Table [TT] the highest and lowest S /N values for an exposure 
of 2595 s. Not surprisingly, the short exposure of 707 s (due to 
a technical problem) was useful for the brightest objects only. 
For a given binary with fgxp = 2595 s, the ratio of the highest 
S/N to the lowest S /N is ~2. 

3. Data reductions and analysis 

3.1. Spectroscopic data reduction 

The basic reduction and calibration steps including ve- 
locity correction to the heliocentric reference frame 
for the spectra were performed with the GIRAFFE 
Base Line Data Reduction Software (BLDRS) (see 
http://girbldrs.sourceforge.net). Sky subtrac- 
tion, a critical step for faint objects, was done as follows: 
for each epoch an average sky spectrum was computed from 
the 21 sky spectra measured over the whole FoV. For a given 
epoch the sky level was found to vary slightly across the 
field, but interpolating between spectra was not considered 
a valuable alternative. Local sky variations with respect to 
the average spectrum are given in Table |8] The values Asky 
are read as follows: for example, the sky position labeled 
S19 is on average (i.e. over all epochs) ~20% brighter than 
the mean (i.e. over all sky positions) sky spectrum. The 
variations were found to be between about -11 and -1-35%. 
Normalization to the continuum, cosmic -rays removal and 
Gaussian smoothing {FWHM = 3.3 pix) were performed with 
standard NOAO/PyRAF tasks. 

The first 60 A of the spectra, i.e for wavelengths between 
3940 and 4000 A, were suppressed. The reason is that below 
~4000A the S /N is getting very poor and therefore there is no 
reliable way to place the continuum. Furthermore, the region 
around He was found to be strongly contaminated by the inter- 
stellar Can H and K absorption lines. The last few A (above 
4565 A) were equally found unusable because of a strongly 
corrupted signal. 



3.2. Analysis 

For historical reasons, the analysis has been made in essentially 
two steps. 

First, RG did a complete analysis of all systems, using the 
KOREL code (Hadrava 1995 2004) to obtain both the radial 
velocity curves and the disentangled spectra of the individ- 
ual components. Then, the simultaneous analysis of light and 
RV curves was made with the 2003 version of the Wilson- 
Devinney (WD) Binary Star Observables Program (Wilson & 
Devinney 197^, Wilson [19791 [1990) via the PHOEBE inter- 
face (Prsa & Zwitter [20051 1 ■ However, simulations performed 
following the referee's request regarding this early version of 
the work, showed that the amplitude Kp of the RV curve and 
the mass ratio q = MslM-p - Kp/K^ were not recovered with 
the expected robustness. More details about these simulations 
are given below (Section r3.5.1b . Although, on average, the right 
values are recovered, one particular solution may be off by as 
much as five percent (one sigma) or ten percent (two sigma), 
which was deemed unsatisfactorju. 

Thus, a second, almost independent analysis was made by 
PN, using a least-squares fit of synthetic binary spectra to ob- 
served unsmoothed spectra for the RV determinations. The lat- 
ter technique seemed much more robust, according to the same 
simulation: the Kp and K^ amplitudes are recovered to better 
than one percent - at least for the particular binary that was 
simulated - and the dispersion of the values of the small eccen- 
tricity (e ~ 0.03) is no larger than two percent. Small system- 
atic errors may result from temperature or rotation mismatch, 
but they remain smaller than the uncertainties of the previous 
analysis. 

We took the opportunity of this new analysis to define a 
more objective determination of the effective temperatures of 
the components. Instead of a visual fit of a spectrum close to 
quadrature, we used a least-squares fit of synthetic binary spec- 
tra to all observed spectra falling out of eclipses. Then, the error 
on the effective temperatures could be naturally defined as the 
RMS scatter of the results. More details are given in the next 
sections and in the following discussion of individual binaries. 
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Fig. 2. Histogram of periods of our sample of 33 eclipsing bi- 
naries in 0.5 day bins. 



3.3. Pliotometry: quality check 

The quality of the /-band light curves was discussed in Section 
|2| Despite the high range of RMS scatter, we can expect a very 
accurate determination of the out-of-eclipse /-magnitude be- 
cause of the large number of data points (~300). This is not 
the case with the B- and V-bands. A much poorer sampling 
can lead to erroneous zero-level computation in the light-curve 
analysis step and result in wrong colour-index determination. 
Therefore, it is necessary to perform a quality check of the pho- 
tometric data. This was done in the form of colour-colour dia- 
grams of our sample. Figure |3]presents the three diagrams that 
can be obtained from the three colour indices (B - V)'', (V - /)'• 
and (B - /)''. These are the values at quadratures, i.e an aver- 
age value characterizing a "hybrid" star of intermediate proper- 



' This should not be interpreted as a criticism of the KOREL code, 
but only as a warning that this code should be used in a very careful 
way. See also Subsection l3.5.4l for a remark about the other method. 
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Table 1. Spectroscopic observations: epochs, heliocentric JuUan dates, exposure times and sky conditions. The value of the air 
mass is given at mid exposure. The value of the seeing is an average between the start and end values. The moon age is the 
number of days elapsed since the last new moon. 



Date of observation (start) 


HJD 


'exp 


Air mass 


Seeing 


Age of Moon 




(-2450000) 


(s) 




(") 


(d) 


2003-11-16X00:39:59.373 


2959.5423 


2595 


1.526 


0.77 


21.00 


2003-11-16X04:43:39.664 


2959.7115 


2595 


1.744 


0.53 


21.16 


2003-11-17X00:35:51.283 


2960.5394 


2595 


1.526 


1.27 


21.96 


2003-11-17X03:55:11.067 


2960.6778 


2595 


1.645 


0.95 


22.09 


2003-11-18X00:27:41.157 


2961.5336 


2595 


1.529 


0.80 


22.94 


2003-11-18X04:40:08.506 


2961.7090 


2595 


1.755 


0.96 


23.12 


2003-11-19X00:41:56.657 


2962.5435 


2595 


1.519 


1.00 


23.97 


2003-11-19X05:06:05.545 


2962.7269 


2595 


1.845 


0.67 


24.16 


2003-11-20X00:27:32.704 


2963.5335 


2595 


1.524 


1.01 


25.02 


2003-11-21X00:46:40.291 


2964.5358 


707 


1.518 


0.98 


26.13 


2003-11-21X01:00:00.631 


2964.5560 


2595 


1.512 


0.83 


26.14 


2003-11-21X05:39:31.454 


2964.7501 


2595 


1.998 


0.74 


26.35 


2003-11-22X00:20:45.386 


2965.5287 


2595 


1.523 


0.94 


27.23 


2003-11-22X04:32:44.400 


2965.7037 


2595 


1.779 


0.92 


27.43 


2003-11-23X00:44:54.938 


2966.5454 


2594 


1.514 


0.68 


28.40 


2003-11-23X04:27:51.663 


2966.7002 


2595 


1.778 


1.02 


28.58 



ties with respect to the two components of a particular system. 
Not surprisingly, most binaries are found on a relatively narrow 
linear strip. For any diagram, the scatter of the objects is low 
because the reddening hne is almost parallel to the sequence. 
For example, the ratio Eb-v/Ev-i = 0.81, determined from Eq. 
[TTlin Section [3. 121 is close to 0.69, the slope of the sequence 
(B - V)'* vs. {V - Z)''. Nevertheless, four outliers appear, which 
are marked with open symbols. In principle, there is a possibil- 
ity to restore a bad colour index, as illustrated by the example 
of the system 5 261267: Inspecting the three light curves, one 
can suspect that the cause of the discrepancy lies in a poorly 
sampled V-band light curve in the out-of-eclipse domain. The 
other two light curves (B and 7) seem more rehable. Therefore, 
only the B - I colour-index is reliable for this system. But the 
two other indices, B - V and V - I, can safely be interpolated 
from the B - I value under the assumption that the system lies 
on the linear strip. The method is illustrated by the dashed lines 
in the diagrams. Of course, this reconstruction of two bad in- 
dices from a good one is not possible, unless only one of the 
three light curves is unreliable (either B or V). In this particular 
example, the situation is not so clear-cut because the recon- 
structed indices (B - V = -0.129 and V - / = -0.139) would 
imply too low a V value. Therefore, all four outliers will be 
excluded in the final estimate of the distance modulus. 

On the upper right diagram of Fig. [3] red crosses indicate 
the intrinsic colours computed below (Section[TTT]). Their po- 
sitions, slightly below the regression line of the sequence, is 
entirely compatible with the reddening arrow and the observed 
colours, which inspires confidence regarding the colour ex- 
cesses determined in Section [3.11l 
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Fig. 3. Observed (reddened) colour-colour diagrams of our 
sample. These are the brightest out-of-echpse values from 
Table |3] They were determined from the best light-curve fits. 
The low scatter is due to a reddening line almost parallel to the 
main sequence (see upper right panel). Note that some values 
(open symbol) are doubtful as they lie off the linear trend. The 
dashed lines illustrate how to recover correct B -Y and Y - I 
indices from doubtful ones, under the assumption that the B-l 
index is reliable, for the system 5 261267. 



3 A. Synthetic spectra 

Except for the first two steps of the analysis, i.e. (a) the simul- 
taneous disentangling of the composite spectra and retrieving 



of the RV curves through the KOREL code (Hadrava 119951 
120041) and (b) the simultaneous analysis of the light and RV 



p. North et al.: VLT spectroscopy of eclipsing binaries in the SMC 



curves through the WD/PHOEBE code, the search for the pa- 
rameters of a binary reUes heavily on synthetic spectra. Indeed, 
the systemic velocity, the projected rotational velocities, the ra- 
tio of radii and the primary temperature are found by compar- 
ing the observed spectra (disentangled and composite as well) 
with a library of synthetic spectra. Actually, two libraries were 
used. For objects with effective temperature 15 000 K < Tetf < 
30 000 K, we used the BSTAR2006 Ubrary recently released by 
Lanz & Hubeny ('2007'). For the few objects above 30 000 K, we 
took the OSTAR2002 library previously released by the same 
authors (|2003] |. Both librai-ies are available at the TLUSTY 
Web site (verb+http://nova.astro.umd.edu/index.html-i-). 

The grids with a metallicity suitable for the SMC, 
Z = 0.004, i.e. one-fifth of the solar metallicity, were cho- 
sen. Concerning the BSTAR2006 library, we took the grid 
with a microturbulent velocity of 2 km s"'. Both grids of 
spectra are based on NLTE line-blanketed model atmo- 
spheres. The temperature step is 1000 K below 30000 K 
and 2500 K for early-O stars. We restrained to surface 
gravities 3.25 dex < log^ < 4.75 dex (0.25 dex step) . 
The spectra were convolved with the appropriate rotational 
profiles (Vrotsin/ = 30, 40, 50, 75, 100, 150, 200, 250 
and 300 km s ') and with a Gaussian instrumental pro- 
file (resolution of 0.67 A) by mean of the program R0TIN3 



3.5. Radial velocities 

As explained above (Section 13.21 ). the radial velocities were 
determined using two independent methods. As a first step, a 
disentangling method was used, which allowed to give a first 
estimate of the parameters of all systems. As a second and final 
step, a least-squares method was used. Each step is described 
in turn below. 



3.5.1 . First Step: spectrum disentangling 



Simon & Sturm (1 19941 ) were the first to propose a method al- 
lowing the simultaneous recovery of the individual spectra of 
the components and of the radial velocities. Another method 
aimed at the same results, but using Fourier transforms to 
save computing time, was proposed almost simultaneously by 
Hadrava ( 1995 ). The advantages of these methods are that they 
need no hypothesis about the nature of the components of the 
binary system, except that their individual spectra remain con- 
stant with time. Contrary to the correlation techniques, no tem- 
plate is needed. In addition to getting at once the radial ve- 
locities and orbital elements, one gets the individual spectra 
of the components ("disentangling"), with a signal-to-noise ra- 
tio which significantly exceeds that of the observed composite 
(http://nova.astro.umd.edu/Synspec43/synspec-frames-rotin.html)spectra. For instance, in the case of a binary system with two 



Beside a grid of continuum normalized spectra, a grid of flux 
spectra was generated for colour indices calculation through 
synthetic photometry. 

Formally, a normalized synthetic composite spectrum is 
computed at a given orbital phase (p from the normalized 
^-component spectra F'}^^^ ^, the radial velocities of both 
components Vj{4>,p) and the light dilution factors tj^cXnA'- 



KxnM) 



Yj ij,cX^M)FlcinA * 6{c\nA- Vj(4>,p)) (1) 



6 is the Dirac (J-function and p is the set of orbital parameters. 
The spectra are expressed in the logarithmic wavelength scale. 
The light dilution factors were calculated from light-curve 
modeling (Ilijic et al. 120041 ): 



f-pj 



(0) 



1 



l+X^frS 



and {sA'P) = 1 - 4,^(0) (2) 



epW 



where ej((p), the flux of component /' in a given passband 
normalized by the value at quadrature, accounts for possible 
out-of-eclipse variability due to departure from sphericity and 
reflection effects (Fitzpatrick et al. 120031 ): Xi is the ratio of the 
continuum monochromatic luminosities, i.e. the ratio of the 
mean surface brightnesses 7j ,; times the ratio of radii. 



e,(4>) = 



Am 



, r ^S,A Js,A I Rs 

and i,,! = = — 

Lpj Jpj \Rp 



(3) 



We checked that X.a can safely be considered constant through 
the 4000-4580 A spectral interval and can be safely identified 
with the luminosity ratio in the Cousins B-band of effective 
wavelength /letr = 4360 A. 



components of equal brightness and observed 16 times with 
S /N = 50, the signal to noise ratio of each disentangled spec- 
trum would be S/N = 4 x 50 x 0.5 = 100 (the factor 0.5 being 
due to the fact that there are two components). That is impor- 
tant, because the nature of the components can then be deter- 
mined safely. In this work we use this advantage to determine 
the effective temperatures of some components, but for brighter 
binaries observed at higher resolution and S/N, it would also 
be possible to determine photospheric abundances. Other de- 
tails about these techniques and their applications (including 
abundance determinations) can be found in e.g. Hensberge et 
al. (12000b . Pavlovski & Hensberge (120051) and Hensberge & 
Pavlovski (!2007l) . 

The radial velocities were determined from the lines of He i 
(/14471, ^4388, ^4144, /J4026) only. We preferred to avoid the 
H Balmer lines (as did Fitzpatrick et al. 2002) because (1) 
their large width make them more sensitive to systematics due 
e.g. to wrong placement of the continuum, and (2) because 
of moderate to strong nebular emission polluting most sys- 
tems (only 6/33 systems were found devoid of any emission). 
Consequently, four regions with a width of 80 A centered on 
the four Hei lines were cut from each spectrum of the series. 
Since the KOREL code makes use of the Fourier transform of 
the spectra, both edges of each spectral region were fixed to 
1 by hammering of the signal to 1 with a cosine bell function 
(Hanning window). KOREL was run with out-of-eclipse spec- 
tra only, although the line-strength factors, i.e. the contributions 
of each component to the system continuum, could be obtained 
in principle as results of the KOREL analysis. However, most 
of our spectra were found to have a too low S /N to provide re- 
liable results. Therefore, the selection was performed using the 
out-of-eclipse phase ranges given by the light curves. The pe- 
riod P was taken from Table[3] which are slightly improved val- 
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ues relative to those of Wirzykowski et al. (120041 1. as explained 
in SectionlJl A first estimate of the epoch of periastron passage 
To, of the eccentricity e and of the longitude of periastron wq 
were determined from the light curves. In the case of eccentric 
systems, a first solution was found neglecting apsidal motion. 
The only orbital parameters allowed to converge were the pri- 
mary semi-amplitude Kp and the mass ratio q. For each system, 
KOREL was run with a grid of values (Kp, q). The solution with 
the minimum sum of squared residuals as defined by Hadrava 
(I2004I I was retained as the best solution. For eccentric systems, 
a second run was performed letting Kp, q, Tq and u free to con- 
verge (e being determined by photometry). It is important to 
notice that the four spectral regions were analyzed simultane- 
ously, i.e in a single run of KOREL. Each region was weighted 
according to the S IN of each Hei line (weight oc (S /Nfy. To 
circumvent the difficulty to measure S /Ns inside the lines, they 
were estimated from the values calculated with the GIRAFFE 
Exposure Time Calculator of the ESO. The calculated values 
were then normalized to the measured value between 4195 and 
4240 A. The non-Keplerian correction and Rossiter effect were 
calculated from the WD/PHOEBE solution. 

Beside the simultaneous retrieving of RV curves, orbital pa- 
rameters and disentangled spectra, the KOREL code is able to 
disentangle spectra for a given orbital solution (Kp, q, Tq and w 
fixed). A final run of KOREL with this mode was then used to 
disentangle the regions around the Balmer and He ii 4200 and 
4542 lines. Indeed, He n lines and a number of Si m-iv lines are 
very useful to set the temperature of hot components. 



3.5.3. Second step: least-squares RV determination 

On the basis of the first analysis, we selected two synthetic 
spectra from the OSTAR2002 and BSTAR2006 libraries for the 
two components of each system, with the parameters closest to 
the estimations. A chi-square was computed as the quadratic 
sum of the differences between the observed spectrum and 
the composite synthetic one, for arbitrary radial velocities. 
However, we did not use the complete spectra: the hydrogen 
Balmer lines were suppressed because of their large width 
and because they are mixed with nebular emission in a num- 
ber of cases. A SuperMongo (Lupton & Monger [20001 1 pro- 
cedure implementing the amoeba minimization algorithm was 
used, letting the two radial velocities and the blue intensity ra- 
tio free to converge. The radial velocities are essentially con- 
strained by five He i fines (.14009, 4026, 4144, 4388 and 4471). 
Convergence was generally fast and robust, in the sense that 
the results did not depend on the initial guess values. Some 
iterations were necessary, however, to clearly identify the pri- 
mary and secondary components, so that the right model be 
attributed to the right component. 

A preliminary analysis of the radial velocities was then 
performed using an interactive code (Lucke & Mayor [1980]l, 
which allowed to assess the quality of the RV curves (espe- 
cially the RMS scatter of the residuals) and obtain first orbital 
elements. 



3.5.4. Simulations 



3.5.2. Testing the robustness by simulations 



In order to test possible biases on the determination of Kp and 
q by KOREL, we have simulated ten sets of nine out-of-eclipse 
composite spectra of the system 5 266131. We used the fitted 
radial velocity curves to shift the synthetic spectra of each com- 
ponent and add them at the observed phases, using the adopted 
luminosity ratio in the B band. We used the synthetic spectra 
with parameters closest to the observed ones, from the grid of 
Munari et al. (12005 1) . A Gaussian noise was added to the com- 
posite synthetic spectra, such that the signal-to-noise ratio var- 
ied from 37 to 71 (see Table fTTTi and assuming that the S IN ra- 
tio is inversely proportional to the seeing given in Table [1] The 
KOREL code was run on each of these ten simulated datasets, 
and the averages of the resulting Kp and q values were com- 
puted. We found Kp - 228 kms"' and q - 0.877 instead of 
the input values 223.4 + 4.8 km s"' and 0.889 ± 0.039 respec- 
tively. The simulation then gives K^ - Kplq - 260 kms"' 
instead of the input 263.0 + 3.9 kms '. Thus, the parameters 
obtained from the simulated spectra agree with the input value 
to within Icr, so that no significant systematic error is to be 
feared. However, the RMS scatter of the Kp and q parameters 
proved disappointingly large, about 10.8 kms"' and 0.044 re- 
spectively. This means that the uncertainty on the RV amplitude 
reaches about 5%, which translates into 15% for the masses. 



The same ten sets of nine composite synthetic spectra, de- 
scribed above, was used to test the least-squares method of RV 
determination. The results proved very encouraging, since they 
follow distributions, the sigma of which amount to only 0.8 and 
0.6% of the means, for the amplitudes Kp and Ks respectively. 
The sigma of the eccentricity distribution is 1.8% of e = 0.036, 
and the argument of the periastron has cr(w) = 11°. 

The effect of a mismatch was tested by using template 
spectra with effective temperatures smaller by 3300 K, re- 
spectively 4080 K for the primary and secondary components, 
compared to the temperatures used to build the artificial "ob- 
served" spectra (the projected rotational velocities were also 
smaller by about 40 kms"' ). The amplitudes changed by 
+ 1.1, resp. +0.3% only for the primary and secondary com- 
ponents. Increasing the temperatures by 2700, resp. 1860 K 
(and the Vrotsin/ by 40 kms"' ) lead to relative differences 
AKplKp = -1.1% and AKslKs = -1.4%. Thus, the mismatch 
that can be expected will not induce systematic errors much 
larger than about one percent, which is in general smaller than 
the random error. In view of its robustness, we have adopted 
the least-squares technique for RV determination, rather than 
the results of the KOREL code. However, we are aware that the 
above comparison between the two techniques may not be quite 
fair, because the KOREL code recovers the individual spectra 
from the data, while the least-squares fit uses external template 
spectra. In that sense the advantage of the least-squares fit may 
prove somewhat artificial. 
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3.6. Apsidal motion 

The WD code allows to determine the time derivative oj of the 
argument of the periastron. That possibility was used for all ec- 
centric systems but one (4 175333), the latter having the small- 
est eccentricity of all. The systems for which a significant ap- 
sidal motion was found in that way were examined further by 
subdividing the photometric data into four consecutive time se- 
ries, and examined using an interactive version of the EBOP16 
code (Etzel I1980I I. The a> angle (if precise enough) and the 
e cos cl) quantities ajusted for each time series were then exam- 
ined for systematic variation. That allowed to better visualize 
the effect of apsidal motion on the light curves, and to assess 
better its significance. 



3.7. Wilson-Devinney analysis 

3.7.1. First step 

For each system, a preliminary photometric solution had been 
found (before taking thre radial velocities into account) by the 
application of the method of multiple subsets (MMS) (Wilson 
& Biermann'1976). The groups of subsets used were essentially 
the same as those advocated by Wyithe & Wilson (2001 , 2002) 
(Table |2l). That allowed to provide fairly precise values of e 
and ll) that were then introduced in the KOREL analysis. Then, 
all three light curves and both RV curves provided by KOREL 
were analyzed simultaneously using the WD code. That does 
not imply, however, that photometric and spectroscopic data 
were analyzed in a really simultaneous way, since results from 
the preliminary Ught curve analysis were used in the KOREL 
analysis; it is rather an iterative analysis. The / light curve is 
the most constraining one, thanks to the large number of points, 
but the B and V light curves are very important too, since they 
provide accurate out-of-eclipse B and V magnitudes. The mass 
ratio q was fixed to the value found by KOREL. The semi- 
major orbital axis a, treated as a free parameter, allows to scale 
the masses and radii. In a first run, the temperature of the pri- 
mary was arbitrarily fixed to 26 000 K. Second-order parame- 
ters like albedos and gravity darkening exponents were fixed to 
LO. Metallicities [M/H] were set at -0.5. The limb -darkening 
coefficients were automatically interpolated after each fit from 
the van Hamme tables (van Hamme|T993]l. The code needs an 
estimation of the standard deviations of the observed curves in 
order to assign a weight to each curve. These cr's were calcu- 
lated from the sums of squares of residuals of the individual 
curves, as advocated by Wilson and van Hamme (I2004I I. These 
values were refined for subsequent runs. 

A fine tuning run was performed with the primary temper- 
ature found after analyzing the observed spectra. The standard 
uncertainties on the whole set of parameters were estimated in 
a final iteration by letting them free to converge. 

The standard procedure described above is sufficient for 
symmetric light curves only. For systems displaying a small de- 
pression before the primary minimum, as it is occasionally the 
case with semi-detached systems, it is necessary to introduce 
a cool spot on the primary component. Obviously, this step is 
performed after obtaining the symmetric best-fit solution. The 



spot is characterized by four parameters, i.e. colatitude, longi- 
tude, angular radius and temperature factor As the observed 
feature can be described by a large number of combinations of 
the four parameters related to the spot (high degeneracy), the 
spot was arbitrarily put on the equator of the primary (i.e. co- 
latitude of 7r/2) and the three other parameters were optimized 
alternately following the MMS. In case of high propensity to 
diverge, one of the three free parameters was set to an arbitrary 
value, the MMS being performed on the two remaining param- 
eters. 

In this first step, the WD analysis was performed using 
the photometric convention, according to which the primary 
star is the one that is eclipsed near phase zero, i.e. the star 
with the higher mean surface brightness in a given passband 
(Jpj > 7s,^)- It followed that in some cases q may not necessar- 
ily be < 1 . In this paper, we have finally adopted the dynamic 
convention ^ < 1 in order to avoid confusion. 

For the detached systems, the orbit was considered as cir- 
cular when the eccentricity given by the WD code was smaller 
than its estimated error. 

Table 2. Method of multiple subsets: groups of parameters al- 
lowed to converge for each model. Porb, to, Tp and q are fixed; 
e is fixed to for clearly non-eccentric systems; a may be in- 
cluded in any subset as it does not correlate with any other pa- 
rameters. See text for treatment of asymmetric light curves. 



Subset 1 


Subset 2 




Subset 3 


Detached systems 


(e) 
A(p 

J B.V.I 


i 
Sis 




Ts 

np 

(Wo) 




semi-detached and overcontact 


systems 




A0 

, B.V.I 

Ts 









3.7.2. Second step 

Both photometric and RV curves were analyzed simultane- 
ously, fixing the effective temperature of the primary compo- 
nent to the spectroscopic value (see below for the determination 
of the latter). For semi-detached and contact systems, there is 
no need to fix any other parameter For detached systems, how- 
ever, the ratio of radii is very poorly constrained by photom- 
etry alone when the eclipses are partial, which is the case of 
all detached systems in our sample. Therefore, we adopted the 
brightness ratio determined by spectroscopy, and fixed the po- 
tential of the primary, Qp, to a value such that the brightness 
ratio in the blue band matched the spectroscopic one within the 
uncertainties. The potential depends on both radius and mass, 
but the latter is constrained by the RV curve, so that fixing a po- 
tential is equivalent to fixing a radius. In some cases it was not 
possible to reproduce the spectroscopic brightness ratio with- 
out degrading the photometric fit, so we gave priority to the 
latter 
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Table 3. Basic parameters of the observed eclipsing binaries: identifying number (Fig. 1), OGLE identifying code, coordinates, 
orbital period, epoch of primary minimum (see text), infrared and (B - V) and (V - I) colour indices. The difference between 
/(DIA) and /(DoPhot) photometry is indicated as well. The "q" superscript denotes values at a quadrature (the brightest in case of 
asymmetric light curve), a, 6 are from Wyrzykowski et al. ( 120041 ). while Porb and to are updated values based on both photometry 
and radial velocities (see text). The uncertainty on the period is given between parentheses as the value of the last digit or, in a 
few cases, of the last two digits. P, V, (B - V)'* and (V - 1)'^ were determined from the light-curve solutions. The colour indices 
between parentheses are unreliable: they correspond to the outliers in Fig. [3] 



id 


OGLE 


a (J2000) 


6 (J2000) 


^orb 


to (HJD 


'dia 


A/'' 

"^DIA-DoP 


(B - V)i 


(V-D'f 




object 


(hms) 




o / // 


) 


(d) 


-2450 000) 


(mag) 


(mag) 


(mag) 


(mag) 


1 


4 110409 


00:47:00.19 


-73 


18 


43.1 


2.973170(4) 


619.49136 


15.840 


-0.002 


-0.036 ± 0.007 


0.006 ± 0.004 


2 


4 113853 


00:47:03.76 


-73 


15 


19.8 


1.320757(4) 


620.90811 


17.340 


-0.008 


-0.002 + 0.006 


0.045 ± 0.003 


3 


4 117831 


00:47:31.74 


-73 


12 


02.2 


1.164566(2) 


621.36981 


17.799 


0.003 


-0.115 ±0.004 


-0.102 ±0.003 


4 


4 121084 


00:47:32.05 


-73 


09 


08.4 


0.823722(1) 


624.39596 


16.959 


-0.000 


-0.084 ± 0.004 


-0.052 ± 0.004 


5 


4 121110 


00:47:04.60 


-73 


08 


40.1 


1.111991(1) 


622.29034 


17.003 


0.009 


-0.090 ± 0.005 


-0.076 ± 0.002 


6 


4 121461 


00:47:24.69 


-73 


09 


35.5 


1.94670 


624.3954 


17.926 


0.003 


-0.060 ± 0.005 


-0.081 ± 0.005 


7 


4 159928 


00:48:13.53 


-73 


19 


30.8 


1.150460(2) 


621.13880 


16.704 


-0.000 


-0.079 ± 0.002 


-0.059 ± 0.002 


8 


4 160094 


00:48:10.17 


-73 


19 


37.1 


1.699634(66) 


620.04883 


17.125 


-0.000 


-0.163 ±0.002 


-0.176 ±0.002 


9 


4 163552 


00:47:53.24 


-73 


15 


56.5 


1.545811(2) 


620.73188 


15.771 


0.008 


-0.045 ± 0.003 


-0.030 ± 0.002 


10 


4 175149 


00:48:34.80 


-73 


06 


52.6 


2.000375(3) 


623.85898 


14.970 


-0.003 


-0.218 ±0.007 


-0.238 ± 0.004 


11 


4 175333 


00:48:15.38 


-73 


07 


05.3 


1.251126(9) 


622.86576 


17.732 


-0.009 


-0.121 ±0.003 


-0.122 ±0.003 


12 


5 016658 


00:49:02.93 


-73 


20 


55.9 


1.246158(2) 


466.70225 


17.446 


-0.005 


-0.123 ±0.004 


-0.133 ±0.002 


13 


5 026631 


00:48:59.84 


-73 


13 


28.8 


1.411680(1) 


465.98392 


16.242 


-0.001 


-0.137 ±0.003 


-0.168 ±0.001 


14 


5 032412 


00:48:56.86 


-73 


11 


39.7 


3.607857(1) 


464.67202 


16.318 


0.005 


-0.033 ± 0.006 


-0.023 ± 0.002 


15 


5 038089 


00:49:01.85 


-73 


06 


06.9 


2.389426(2) 


468.55092 


15.256 


-0.002 


(-0.216 ±0.011) 


(-0.162 ±0.007) 


16 


5 095337 


00:49:15.34 


-73 


22 


05.8 


0.904590(1) 


466.18186 


17.090 


0.004 


-0.118 ±0.052 


-0.154 ±0.030 


17 


5 095557 


00:49:18.07 


-73 


21 


55.3 


2.421185(21) 


466.80139 


17.440 


-0.012 


-0.100 ±0.003 


-0.098 ± 0.003 


18 


5 100485 


00:49:19.86 


-73 


17 


55.6 


1.519124(1) 


467.15922 


17.150 


-0.008 


-0.138 ±0.008 


-0.162 ±0.003 


19 


5 100731 


00:49:29.33 


-73 


17 


57.9 


1.133344(3) 


467.82186 


17.378 


0.000 


-0.111 ±0.003 


-0.140 ±0.002 


20 


5 106039 


00:49:20.00 


-73 


13 


37.3 


2.194069(5) 


465.38253 


16.695 


-0.005 


-0.050 ± 0.033 


-0.016 ±0.019 


21 


5 111649 


00:49:17.19 


-73 


10 


24.5 


2.959578(3) 


470.15054 


16.726 


0.007 


0.008 ± 0.009 


0.032 ± 0.003 


22 


5 123390 


00:49:22.66 


-73 


03 


42.8 


2.172917(41) 


464.12108 


16.203 


-0.001 


-0.162 ±0.015 


-0.171 ±0.008 


23 


5 180185 


00:50:02.63 


-73 


17 


34.4 


5.491165(95) 


469.37759 


17.321 


0.034 


(-0.115 ±0.011) 


(0.029 ± 0.009) 


24 


5 180576 


00:50:13.44 


-73 


16 


33.1 


1.561124(2) 


466.91033 


17.607 


0.006 


-0.071 ±0.009 


-0.036 ± 0.003 


25 


5 185408 


00:50:24.52 


-73 


14 


56.0 


1.454991(2) 


466.28931 


17.524 


-0.006 


-0.122 ±0.003 


-0.126 ±0.002 


26 


5 196565 


00:50:30.17 


-73 


07 


38.2 


3.942732(12) 


468.26098 


16.942 


0.003 


(no data) 


(no data) 


27 


5 261267 


00:51:35.04 


-73 


17 


11.5 


1.276632(2) 


464.97658 


16.833 


-0.005 


(-0.129 ±0.006) 


(-0.099 ± 0.004) 


28 


5 265970 


00:51:28.13 


-73 


15 


17.6 


3.495685(54) 


465.39125 " 


16.226 


0.005 


-0.136 ±0.005 


-0.164 ±0.002 


29 


5 266015 


00:51:16.82 


-73 


13 


01.9 


1.808925(2) 


465.10449 


15.964 


0.002 


-0.086 ± 0.004 


-0.072 ± 0.003 


30 


5 266131 


00:51:35.81 


-73 


12 


44.8 


1.302945(22) 


465.50898 


17.119 


0.001 


-0.103 ±0.016 


-0.104 ±0.004 


31 


5 266513 


00:50:57.49 


-73 


12 


30.3 


1.107510(2) 


467.15823 


18.066 


0.005 


-0.035 ± 0.005 


0.009 ± 0.003 


32 


5 277080 


00:51:11.68 


-73 


05 


20.3 


1.939346(4) 


465.96082 


16.070 


0.002 


(-0.019 ±0.013) 


(-0.171 ±0.003) 


33 


5 283079 


00:50:58.67 


-73 


04 


35.8 


1.283583(1) 


466.92376 


17.422 


-0.002 


-0.073 ± 0.004 


-0.013 ± 0.002 



Incomplete observations for one of the eclipse. 



3.8. Systemic velocity and projected rotational 
velocities 

3.8.1. First step 

The component spectra of the four regions centred on the He i 
lines were normalized with the help of the KORNOR program 
(Hadrava l2004l l. The systemic velocity Vy was found from the 
disentangled spectra of the four regions centred on the Hei 
lines. The observed spectra were cross-correlated via the IRAF 
fxcor task against synthetic spectra computed for the estimated 
Teff, \ogg and Viotsin/. The Vy values gathered in Table [TT] 
were obtained as the 5 /A^- weighted averages of the individual 
velocities calculated for each line. 



The projected rotational velocities, V^^'j sin/, were tenta- 
tively measured by calibrating the FWHMs of the He i lines 
against a grid of (FWHM, Vjot sin /) values obtained from syn- 
thetic spectra (Hensberge et al. .2000) . The FWHMs were com- 
puted from Gaussian or Voigt profiles fitting via the IRAF splot 
task. The Vrot sin / values retrieved by this method were of- 
ten found unsatisfactory when comparing observed and syn- 
thetic spectra retrospectively. The problem proved to lie in 
the high sensitivity of the FWHM measurement to the con- 
tinuum placement. Therefore, a synchronous rotational veloc- 
ity was assumed for most circular binaries unless profile fit- 
ting proved this hypothesis wrong. Anyway, this assumption 
is certainly justified for short-period systems, i.e. binaries with 
(R/a) > 0.25 (North & Zahn l2003l l where the ratio (R/a) is 
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the star radius divided by the separation. In the case of ec- 
centric systems, pseudo-synchronization was assumed (Mazeh 
12008 ' Eq. 5.1). For a given star, its pseudo-synchronous ro- 
tational velocity is computed from its radius and the pseudo- 
synchronization frequency of the binary. This equilibrium fre- 
quency, close to the orbital periastron frequency, is given in Hut 
(flMTTl . 

3.8.2. Second step 

Contrary to the first step, when the KOREL code was used, we 
do not need to define the systemic velocity a posteriori here. 
The least-squares method directly provides "absolute" radial 
velocities (i.e. not only relative ones), even though mismatch 
might bias them by a few kms"'. Thus the systemic velocity 
naturally flows from the WD analysis, which includes the RV 
curves. 

As in the first step, rotational velocities were derived from 
the assumption of synchronous (for circular orbits) or pseudo- 
synchronous (for eccentric orbits) spin motion. No clear depar- 
ture from this assumption could be seen on the spectra. 

3.9. Spectroscopic luminosity and ratio of radii 
3.9.1. First step 

As mentioned above, and as emphasized repeatedly by 
Andersen (e.g. 1 19801 1 and rediscovered by Wyithe & Wilson 
(I200n hereafter WWOl), the ratio of radii k = Rs/Rp of an 
EB displaying partial eclipses is poorly constrained by its light 
curve. The ratio of monochromatic luminosities X^i is equally 
not well recovered in fitting light curves of simulated EBs (i.e. 
EBs with previously known parameters). On the contrary, the 
surface brightness ratio and consequently the derived effec- 
tive temperature ratio T^f^/T^ff is, in general, reliably recov- 
ered. The sum of the radii "LRj/a is also very well constrained. 
The poor constraining of k is very well illustrated by Fig. 3 in 
Gonzalez et al. (l2005l l. 

Since our sample comprises only systems with partial 
eclipses, X.^, then k must be determined in order to find 
reliable radii and surface gravities. We followed the procedure 
described in Gonzalez et al. (I2005I I. The ratio of the monochro- 
matic luminosities can be expressed by Eq. l4l (Hilditch i200 1 i : 



equivalent width (EW) of line / for component j. The 'obs' in- 
dex means that the EWs are measured in an observed (out of 
eclipse) spectrum of the binary. These "apparent" EWs are then 
normalized by the true values measured in synthetic component 
spectra. The dependence of the true EWs on the effective tem- 
perature, surface gravity, metallicity and helium abundance (Y) 
is emphasized. 

In principle, from the analysis of a set of four lines, know- 
ing the masses, metallicities and helium abundances, it should 
be possible to derive a purely spectroscopic solution for the two 
radii and the two effective temperatures. If we further assume 
that the sum of the radii is known from the light-curves analysis 
as well as the temperature ratio, the analysis of only two lines 
is in principle sufficient to determine a mixed photometric- 
spectroscopic ratio of radii. However, since EWs are sensitive 
to a possible continuum misplacement, we preferred to fit the 
observed line profiles with a synthetic composite spectrum or to 
determine the ratio of the EWs of two different lines in the same 
component. The latter methods are more reliable than blind ap- 
plication of Eq.|4]to estimate the effective temperatures, at least 
in case of spectra with moderate S /N. 

For a given chemical composition, true undiluted EWs de- 
pend on both the effective temperature and the surface gravity 
of the stars. Moreover, this dependence is not always mono- 
tonic even if we restrict to late O to early B stars, the He i lines 
having a maximum strength at ~20 000 K. Consequently, in or- 
der to avoid the hassle of working with non-explicit equations, 
for a given line, Eq. |4] was solved with the photometric tem- 
perature ratio and the true EWs values corresponding to the 
photometric \ogg and a first guess of T^^.. X.b is then used with 
T^g/T^g and a first guess of ABCb to compute the ratio of radii 
k. Combining k and "ERj/a, the new Rj and log gj values are 
obtained straightforwardly. The small error introduced in the 



chain Xb 



Ri 



log^j 



fP 

-'eff 



Wj 



because of 



using approximate values for the true EWs could be removed 
after iterating one more time. 

Nevertheless, this method is not very efficient when the ob- 
served EWs have large uncertainties as in case of a composite 
spectrum of low S /N. In this case, a more pragmatic approach 
consists in optimizing both Xs and T^g in a single step by look- 
ing for the best-fitting synthetic composite spectrum for a given 
pair (£b, T^g) and the Uij/a and T^g/T^g constraints. 



£a 






RsY 
Rp) 



wJ(^p„,loggp,Zp,yp)^ 

WfiTlg,loggs,Zs,Ys)) 



rpP 



lQ-OA(BC'',iTl„)-BCliTl„)) 



(4) 



where Rj is the radius, T^^^ the effective temperature and BC^ 
the bolometric correction in the /l-band of component j. In our 
case, A = B and the bolometric correction in the B-band is 
given by BCBiT.g) = BCviT^s) - (B - V)o(T,g), where the 
bolometric correction in the V-band is interpolated from val- 
ues given by Lanz & Hubeny ( 120031 12007 ) and the intrinsic 
colour index is computed from synthetic photometry. Wj is the 



3.9.2. Second step 

Here the luminosity ratio in the blue is simply one of the three 
parameters determined by the non-linear least-squares algo- 
rithm amoeba, the other two parameters being the effective 
temperatures (see more details below, Section [3.10l ). So the lu- 
minosity ratio is determined in a very homogeneous way, and 
an error estimate naturally arises through the RMS scatter of 
the resulting values. This does not guarantee, however, that the 
results are free from any bias. In particular, one may suspect 
that, in the temperature regime where the strength of all lines 
(H and He ones) vary in the same way with temperature, some 
degeneracy may arise between the temperatures and the lumi- 
nosity ratio. Since that temperature regime spans roughly from 
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22000 to 30000, it means that the luminosity ratio of the ma- 
jority of systems may be fragile. Nevertheless, a posteriori ex- 
amination of the resulting HR diagrams does not confirm this 
fear, even though a few systems fail to match the evolutionary 
tracks. 



3.10. Effective temperatures 

Once reliable Viot sin / and log g values have been found, a way 
for setting the temperature of the primary must be found (the 
temperature of the secondary is a by-product via the photomet- 
ric temperature ratio). 

For late O and early B stars, the Hy and Y\6 Balmer lines 
are far better temperature indicators than He i lines (Huang & 
Gies 20051 I2005I I. Therefore, the most direct way to determine 
the effective temperatures of both components of a given sys- 
tem would consist in calibrating the equivalent widths mea- 
sured on the normalized disentangled spectra with those ob- 
tained with a library of synthetic spectra. Unfortunately, this is 
not always possible because of the high proportion of systems 
contaminated by Hy and H6 nebular emission lines. Thus, most 
spectra of individual components are not reliable around the 
Balmer lines. 



3.10.1. First Step 

A safer method consists in comparing an observed composite 
spectrum close to quadrature with a synthetic composite spec- 
trum computed at the same orbital phase. The spectra retained 
for the temperature determination are those with 0.21 < (p < 
0.29 or 0.71 < (p < 0.79. This method is quite sensitive to the 
continuum placement. A low S /N and/or strong emission lines 
can hinder a reliable profile fitting. 

Another method is a variant of the traditional spectral type 
vs. temperature calibration. In the traditional method, the line 
strengths ratio of two lines are measured and compared to the 
values obtained from a series of reference spectra whose spec- 
tral types are known. The effective temperature is then found 
via a spectral type - Teff calibration scale. As emphasized by 
HHH03/05, this technique is efficient for O-Bl stars but far 
less straightforward for later types. Above this limit, the rela- 
tive strength of the Hen 4542 and Sim 4553 lines is a reliable 
tool, as is the relative strength of the Siiv 4089-4116 and He: 
4121 lines. For temperatures below ~29-30000 K, the prob- 
lem is the lack of exploitable metallic lines. Unfortunately, the 
faint Sin lines are totally undetectable. For later B stars, the 
only detectable metallic line is Mg ii 4481, but this line is often 
severely buried in the noise for most components. 

3.10.2. Second step 

The method is the same, qualitatively speaking, as that of the 
first step, consisting in fitting a composite synthetic spectrum 
to the observed one near quadrature. However, the least-squares 
fit method allows in principle to use all out-of-eclipse spectra, 
and provides a much more objective estimate of the tempera- 
tures. Since the radial velocities are known, the only parameters 



which have to be fit are the effective temperatures of both com- 
ponents and the blue luminosity ratio. As mentioned in Section 
13.9.21 the fit is quite robust at both ends of the temperature 
range of our sample: at the cool end, the He i lines increase in 
strength with temperature, while the H i Balmer lines decrease; 
in addition, the Mg ii line decreases very fast. At the hot end, 
both H and He i lines decrease with temperature, but the He ii 
lines begin to appear. In the intermediate range, all lines vary 
more or less in parallel, which may lead to degeneracy when 
the S/N ratio is poor 

For all systems showing a significant nebular emission in 
the core of the hydrogen Balmer lines, we simply removed a 
4 A wavelength interval centered on the emission line, in both 
observed and synthetic composite spectra. But, contrary to the 
synthetic spectra used for RV determinations, here we include 
the H Balmer lines in the fit, except for their very centres. 

The RMS scatter of the fitted effective temperatures is typ- 
ically of the order of 1000 K, and is even smaller than that for 
one third of the sample. Although formally, the error bar on 
reft should be set to that scatter divided by the square root of 
the number of spectra, we chose to put it equal to the scatter 
itself. Indeed, visual examination of the observed and model 
spectra show that the temperature effect is often very subtle, so 
we feel this choice is more realistic. 

The fit proved to depend somewhat on the normalization 
of the spectra. The latter were first normalized using an auto- 
matic procedure with fixed continuum regions. Then, another 
automatic procedure was used, which corrected the first nor- 
malization with the help of a pair of synthetic spectra with pre- 
liminary stellar parameters. That normalization resulted in a 
slightly higher continuum and was found satisfactory in gen- 
eral, except for the bluer end of the spectrum. A final normal- 
ization was made by hand, which was adopted in most cases, 
but not in all because the continuum proved sometimes too 
high. The temperature determination was run on the automati- 
cally normalized spectra as well as on the manually normalized 
ones. The results were found to depend little on the normaliza- 
tion, which was not unexpected since the last two normaliza- 
tions did not differ much from one another 

The average temperatures were computed on all out-of- 
eclipse spectra on the one hand, and on a selection of those 
for which the RV difference is greater than 300 km s ' (250 or 
even 200 km s ' for longer period systems) on the other hand. 
The selection often resulted in a smaller scatter of the temper- 
ature, though not in every case. The average temperature was 
weighted with the inverse of the chi-square provided by the 
amoeba procedure. 

3.11. Synthetic photometry and reddening 

Intrinsic (B - V) colour indices are needed for two purposes: the 
computation of the Eb-v colour excess for a given system and 
the computation of the B-band bolometric corrections BCb of 
the individual components. In the first case, (B - V)o is phase- 
dependent and characterizes the binary as a whole, while in the 
second case (B- y)^ is the usual (constant) colour characteristic 
of a given star. Both types of colour indices were computed 
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from synthetic photometry, i.e. from synthetic stellar spectra 
and the response functions of the filters. The general formula 
for the phase-dependent {B - V)o of a binary is given by 



(B - V)o(</>) = -2.5 log 4 TT 7 + Cb-v 

JSvF'l(cp)dAjSBdA 



(5) 



where 5x is the response function of the X-band filter (Bessell 
I1990I I. F^^ is the synthetic composite flux spectrum at a given 
epoch and Cb-v - 0.606 mag is the zero point (Bessell et 
al. 11998b . The synthetic composite spectra (i.e. the theoretical 
unreddened flux received on Earth) were calculated from the 
synthetic component spectra F^. j (i.e. the surface fluxes), the 
radial velocities Vj((f),p) and the radii Rj: 

1 ^ 
FU<f>) - -T^Yj'^yiA^SicXnA-Vjic^,?)) 

" (t) Z "^J^u * ^(^ 1" -^ - ^Z-^' -P)) ^6) 

j=p 

where wp - 1 and ws - {Rs/Rp)"- Substituting F'j^^ in Eq. ^ 
by the expression given by Eq. ^, one sees that Eq. (|5]l does 
not depend on the distance d to the system. 

A similar procedure was used to compute the (V-I)o colour 
indices, taking Cy-i = 1.268 and the appropriate response 
functions S v and 5 /. This index is needed for the determination 
of the distance modulus from the /-band photometric observa- 
tions. 

Alternatively, the intrinsic colour index (B - V)q could be 
estimated via a colour-temperature relation from the literature. 
This kind of calibration being established with stars from the 
solar neighborhood (e.g. Flower 1996), the coefficients of the 
fit are in turn representative of stars with a solar metallicity. 
Consequently, synthetic photometry with Z - 0.004 spectra 
was considered more reliable for our objects. Both (B- V)o and 
(V - /)o were computed for a grid of temperatures and surface 
gravities. Because the fluxes are not given beyond 7500 A in 
the OSTAR2002 library (which contains fluxes for stars hotter 
than 30 000 K), we had to complement the flux distribution of 
stars hotter than 30 000 K by the appropriate Kurucz fluxe^ in 
order to be able to compute the (V - /)o index. We chose the 
grid with a metallicity of -0.5. Because of that inhomogene- 
ity, the (V - /)o indices of binary systems hosting components 
with Teff > 30 000 K may be slightly less reliable than those of 
the other systems (the systems 4 1 10409, 4 121084, 4 175149, 
5 32412, 5 38089 and 5 266015 are in this case). For a given 
object, the colour indices were linearly interpolated from these 
grids. 

To estimate the uncertainties, however, the following first- 
order approximations valid for early-B stars of intermediate 
logg values were used (see Section [3.13l ). 

(B - V)o ~ 1 .768 - 0.455 log T^g for 15 000 < Teff < 32 OP) 

^ 0.228 - 0.1 13 log Teff for 32 000 < T^g < 37 500 

(y - /)o ~ 1 .958 - 0.502 log T^« for 15 000 < T^s < 32 000 

=s 0.502 - 0.178 log Jeff for 32 000 < T^e < 37 500 



Table 4. Intrinsic colour indices at quadrature computed by 
synthetic photometry. Because of the too short wavelength 
range provided in the OSTAR2002, (V - /)o values of systems 
with a component hotter than 30 000 K were computed using 
Kurucz red fluxes (see text). 



Object 


(B-V)l 


(^-^)o 




(mag) 


(mag) 


4 110409 


-0.243 ± 0.006 


-0.252 ± 0.007 


4 113853 


-0.196 ±0.007 


-0.206 ± 0.008 


4 117831 


-0.179 ±0.004 


-0.192 ±0.004 


4 121084 


-0.271 ± 0.004 


-0.294 ± 0.004 


4 121110 


-0.248 ± 0.006 


-0.269 ± 0.006 


4 121461 


-0.207 ± 0.007 


-0.228 ± 0.008 


4 159928 


-0.247 ± 0.006 


-0.263 ± 0.007 


4 160094 


-0.238 ± 0.009 


-0.260 ± 0.009 


4 163552 


-0.232 ± 0.008 


-0.246 ± 0.009 


4 175149 


-0.271 ± 0.002 


-0.288 ± 0.002 


4 175333 


-0.191 ±0.006 


-0.211 ±0.007 


5 016658 


-0.203 ± 0.005 


-0.218 ±0.005 


5 026631 


-0.242 ± 0.003 


-0.257 ± 0.003 


5 032412 


-0.283 ± 0.001 


-0.303 ±0.001 


5 038089 


-0.273 ± 0.001 


-0.292 ± 0.001 


5 095337 


-0.242 ± 0.009 


-0.261 ±0.010 


5 095557 


-0.216 ±0.003 


-0.233 ± 0.003 


5 100485 


-0.215 ± 0.004 


-0.232 ± 0.005 


5 100731 


-0.217 ± 0.004 


-0.233 ± 0.004 


5 106039 


-0.204 ± 0.003 


-0.210 ±0.003 


5 111649 


-0.169 ±0.002 


-0.168 ±0.002 


5 123390 


-0.254 ± 0.004 


-0.274 ± 0.004 


5 180185 


-0.182 ±0.004 


-0.188 ±0.005 


5 180576 


-0.221 ± 0.009 


-0.237 ±0.010 


5 185408 


-0.222 ± 0.004 


-0.239 ± 0.004 


5 196565 


-0.200 ± 0.003 


-0.209 ± 0.004 


5 261267 


-0.253 ± 0.006 


-0.266 ± 0.007 


5 265970 


-0.214 ±0.001 


-0.226 ±0.001 


5 266015 


-0.265 ± 0.004 


-0.280 ± 0.004 


5 266131 


-0.239 ± 0.006 


-0.258 ± 0.007 


5 266513 


-0.194 ±0.011 


-0.208 ±0.012 


5 277080 


-0.244 ± 0.004 


-0.256 ± 0.004 


5 283079 


-0.234 ± 0.004 


-0.252 ± 0.004 



Once the colour excess Eb-v = (B - V) - (B - V)o is deter- 
mined, one must make an assumption about the value of the ex- 
tinction parameter :Rv = Ay /Eb-v- This parameter is assumed 
equal to the standard value 3.1+0.3 for each system. It is worth 
mentioning that this parameter sufifers from a rather large un- 
certainty. For the SMC bar, Gordon et al. (2003) propose the 
mean value "Rv = 2.74 ±0.13 from a sample of four stars with 
Rv ranging from 2.4+0.3 to 3. 3 ±0.38. Nevertheless, because of 
the small size of Gordon et al.'s sample, the more conservative 
standard value was retained for this paper Moreover, the con- 
tribution ofRv to the total error budget for the DMi is low. The 
diff'erence Adm between a DM calculated with Gordon et al. 's 
value and the standard value is given by Adm = -0.36 Eb-v- 
Thus, with Eb-v ~ 0.1, the DM is slightly increased by 0.036 
mag ifRv = 2.74 is adopted. 



see http://kurucz.harvard.edu/grids.html 
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In order to determine the distance modulus from the /-band 
data, the ratio of the absorptions in the infrared and optical 
bands Ai/Ay is calculated from the relationship, 



Al 1 

-^ = a(x) + b(x) ■— 
Ay Kv 



or A,^EB-v{a(x)'Rv + b{x)) (8) 



with X = l//leff- The coefficients a{x) and b{x) are computed 
from the polynomials given in Eq. 3a-b from Cardelli et al. 
(11989V An effective wavelength of 7980A is taken for the 
Cousins /-band filter (Bessell 1990 ). Hence, it follows that 
AilAv ~ 0.600. That value contrasts with the value of 0.479 
given by Cardelli et al. ( 119891) in their Table 3 for the / band, 
because they adopt another effective wavelength. 

3.12. Distance modulus: V- or 1 -band approach? 

There are two ways of determining the distance modulus of a 
particular system, depending on whether one takes visual or 
infrared data: 



5 log d- 5 



yq _ Ml - Al 




/q _ m; - A] 




I^- Ml + {V-I)l- 


- 0.600 A^ 



(9) 



The safest choice is to take the expression which leads to the 
value with the smallest uncertainty. The /-band magnitude 
at quadrature is far more reliable than the V-band value be- 
cause there is typically ~10x more points in the /-band light 
curve. Moreover, the contribution of the interstellar absorption 
is lower at infrared wavelengths. Indeed, it is easy to demon- 
strate that the variance of the visual absorption is always larger 
than the sum of the contributions due to the infrared absorption 
and the intrinsic (V - /) colour: 



cr^, 



> (T ,^ + cr' 



(V-Dl 



with 



(v-'« 



(B-V)l 



(10) 



Consequently, the /-band calculation must be preferred, at 
least for system for which there is a reliable way to obtain the 
(y - /)o index. We have seen that this is the case of systems 
where both components have T^s ^ 30 000 K, while for hotter 
binaries, the synthetic spectra of the OSTAR2002 library had to 
be complemented by Kurucz models beyond 7500 A. In spite 
of this slight inhomogeneity, we have adopted the (V - /)o in- 
dices computed in this way, rather than refrain from computing 
the / distance modulus for these hot binaries. 

In order to check the purely synthetic intrinsic (V-I) index, 
both V- and /-band formulations of Eq.|9]can be used to derive 
the following relation between the two colour excesses: 



Ev-i 



= 0.400 -Ry 



ifi-y 



A simple modification of this equation gives, 
(V - /)o = (y - /) - OAOO'RvEb-v 



(11) 



(12) 



This is a semi-empirical dereddening, as Eb-v is calculated 
from a synthetic intrinsic {B - V) index. The (V - /)o indices so 
obtained agree within 0.01 - 0.02 mag for 4 of the 6 hot bina- 
ries; the differences reaches 0.03 mag for the system 5 32412, 



but this is less than the estimated error on the distance modu- 
lus. For 5 38089, the difference is 0.06 mag, but this system has 
abnormal colours and so will be excluded from the statistics of 
the distance moduli. 

3.13. Uncertainties 

3.13.1. Distance modulus and related parameters 

The calculus of the uncertainties in a number of parameters 
relies on some assumptions and simplifications. These points 
are discussed in this paragraph. The uncertainty in the distance 
modulus of a particular binary cr5i(,g^_5 is determined via the 
standard rules of the propagation of errors for independent vari- 
ables. 



'51ogd-5 



cr' 



K 



+ ^(v-/)3 + 



(o.eoo'RvfL 



'(B-Vyi + ""(B-Vt 



') 



(13) 



This expression was easily derived from the infrared formu- 
lation of Eq. |9] The uncertainty in the / magnitude (cr/q < 2 
mmag) and the E^y cr^ term are negligible and consequently 
omitted from Eq.[T3] It is worth noting that formally the intrin- 
sic colour indices (B - V)^ and (V - /)q and the visual absolute 
magnitude of the system M^ are not independent variables, as 
they depend ultimately on the effective temperatures and radii 
of the stars. Nevertheless, for the sake of simplicity these vari- 
ables were considered as independent in the following develop- 
ment. The uncertainty on M^ can be expressed as a function of 
the uncertainties on the absolute magnitudes of the components 



M^ and M^ 



(^IQ-O AMI + lO-O AMiy j-^ 



2(10-"'"''-^^)' (14) 



The absolute visual magnitude of a component is found via 
the absolute bolometric magnitudes Mboi and the visual bolo- 
metric correction BCy, 

MviR, reff,Z) = MUR, T,ff) - BCv(T,ff,Z,\ogg) (15) 

There is now an important point emphasized by Clausen 
(I2OOOI 1. Mboi and BCy depend both on the effective tem- 
perature of the star and thus are not independent variables. 
Consequently, from the bolometric corrections calculated by 
Lanz and Hubeny (l2003l |2007] |. we find that for late-O and 
early-B stars the bolometric correction can be given by 

BCv ~ a - p log Teff with 

a = 21.72 ±0.28 mag and /? = 5.51 ± 0.07 mag (16) 



with the associated uncertainty 



P 



In 10 






(17) 



Combining Eq. [15] and [16] and expressing Mboi as a func- 
tion of effective temperature Teff and radius R, the visual abso- 
lute magnitude My can be written in turn as a function of T^s 
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and R (see Clausen [20001 for more details). It follows that the 
uncertainty in My is given by 



1 



Ml 



(In lOf 



5cruV ((10 -P)cr^j 



Ri 



^eff 



(18) 



Thus the uncertainty in Mboi ('10' factor) is partially can- 
celled by the uncertainty in BCy {'P' factor). 

The uncertainties in the intrinsic colour indices at quadra- 
ture, (B - V)o and {V - /)o, are estimated via the approximation 
for early-B stars given in Section l3.11l Since the colour indices 
relate, not to a single star, but to a binary system, one defines an 
"equivalent" effective temperature Teff.syst of the system (which 
is the effective temperature of a single star with the same colour 
index as the system) by inverting Eq.|7] Then, a rough estimate 
of the error on the intrinsic colour index is obtained by propa- 
gating the error on Teff^syst, which is identified with that on T^^, 
the temperature of the primary. In general, the error on the ab- 
solute magnitude of the system will outweigh by far the other 
terms on the right side of Eq. [13] 

3.13.2. Masses and radii 

Masses and radii depend heavily on the radial velocity semi- 
amplitudes Kj. The WD code does not provide directly the er- 
rors on the masses and radii, but only those on the semi-major 
axis a, on the mass ratio q, the inclination /, the effective tem- 
perature of the secondary and the potentials. It is in principle 
possible to derive the errors on the masses from those on a and 
q using the third Kepler law (which gives the total mass) and 
appropriate propagation formulae. However, the resulting er- 
rors tend to appear underestimated, so we preferred to use the 
errors on Mp^s sin / given the bina code that was used for a 
preliminary interpretation of the RV curves alone. The errors 
on the masses are then obtained through the formula 



cr 



Mps 



' Ml 2 sin-' i 



sin"' / 



3 cos / 



O"; 



and may be twice larger than those obtained from WD. Still, 
they have to be considered rather as lower limits to the real 
uncertainties, because they are based only on the scatter of the 
radial velocities around the fitted curve. Systematic errors may 
arise, however, from the choice of the template spectra used to 
obtain the RV values, and from the unavoidable fact that these 
synthetic spectra can never perfectly match the real ones. 

The relation between the uncertainties in the masses and 
pertaining variables can be found in Hilditch (I2001I I. for exam- 
ple. The absolute radius of component j is obtained from 



Rj — a Kj (S,-, k) 



(20) 



where a is the semi-major axis and rj the relative radius of 
component /'. The relative radii are considered as functions of 
two independent variables, the sum Z^ = rp -H rs and the ratio 
k = rs/rp of the relative radii. The sum of the radii is obtained 
reliably from the analysis of the light curves. The ratio of the 
radii is obtained either from the spectroscopic B luminosity ra- 
tio, in the case of detached systems (unless the minima are so 



deep that the photometric estimate proves accurate enough), or 
from the light curve in the case of semi-detached and contact 
systems. The individual radii are thus given by 



1 



rp 



1 +k 
with the associated uncertainties 

\2 



2r and rs —"Lr — rp 



2 ( Q-P.s Y 



(Tv + 



(l+k)^ 



(21) 



(22) 



where ap - I and as - k. E,. is robustly determined from the 
light curve. For the sake of simplicity, the variance in 1,,- was 
taken from the EBOP solution of the / curve, by scaling the 
error on rp obtained by fixing the ratio of radii. 

The spectroscopic ratio of radii k is obtained by inverting 
Eq.a 

Rp ^B \Tp, 
where 



10' 



,0.2ASCb 



(23) 



BCb = BCv -(B- V)o ~ 21.72 - 5.51 log T^g 
-(1.768 -0.455 log Teff) 



ABCb = BCliT^g) - BC^(rfjf) 



-5.055 log I :^ 



^eff 



eff^ 



Therefore, the ratio of radii can be written 



■T-S 

/ii/z I eff 
ell 



^1/2 



0.99 



(24) 
(25) 

(26) 



(27) 



Neglecting the uncertainties on the empirical coefficients in the 
expression for the bolometric correction, the error on k is then: 



m » ife) --' 






CTt^s 
■'eff 



(28) 



(19) 4. The individual binaries 



Each system is discussed thoroughly in this section. We give 
details concerning the light-curve solution, the radial-velocity 
solution, the temperature and luminosity-ratio determinations 
and the characteristics of the spectra. Also discussed are the 
positions of the components in the mass-surface gravity plane 
and the temperature-luminosity (HR) diagram. Review of the 
distances and collective properties of the whole sample of 33 
binaries follows in Section |5] Except where otherwise stated, 
when referring to the light curve of a specific system, it means 
the /-band light curve. 

The /-band light curves and the best-fit solutions are 
shown in Figs [30ll32] The RV curves are shown in Figs |33} 
|35] The mass-log^ diagrams are shown in Figs [36ll38] The HR 
diagrams are shown in Figs [39pn The parameters found from 
the WD/PHOEBE analysis are given in Tables [TT] (orbital pa- 
rameters) and [T2] ( temperature ratios, potentials and luminosity 
ratios). The astrophysical parameters of the primary and sec- 
ondary components are given in Tables[T3]andfT4l respectively. 
Finally, the RMS scatters of the light curves and RV curves are 
summarized in Tables [T6]and[T7] 
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4.1. 4 110409 

With a difference of 0.04 mag in the brightness level between 
phase 0.25 and phase 0.75, this semi-detached system displays 
the most asymmetric light curve among all the systems stud- 
ied in this paper. The light curve is bright (Z'' < 16 mag) and of 
high quality, with a low RMS scatter combined with a deep pri- 
mary eclipse (A/mini/o"/ ~ 65). This EB-type light curve shows 
a relatively strong depression occurring just before the primary 
minimum. Actually, this is strong evidence for absorption by 
a gas stream stemming from the (inner) LI Lagrangian point 
and seen in projection against the primary surface (HHH05). 
As a consequence, the use of a "simple" symmetric model for 
the light-curve fit is not satisfactory, resulting in a rather poor 
fit despite the intrinsic quality of the observations. Therefore, 
this solution was subsequently improved by adding a cool spot 
on the equator of the primary component (see Section 13.7b . 
The parameters of the spot are: a colatitude of n/l rad (fixed), 
a longitude of 0.569 rad, an angular radius of 0.3 rad and a 
temperature factor of 0.6, i.e. the effective temperature of the 
spot is 0.6 that of the rest of the stellar surface. Although this 
new synthetic light curve gives a far more satisfactory fit, the 
O - C curve reveals that this system is certainly more com- 
plex than this "one circular cool spot" model. Actually, there 
are still some discrepancies at the bottom of the eclipses and 
just after the secondary minimum. Nevertheless, this model is 
certainly sufficient to set reliably the inclination, the brightness 
ratio of the components and the maximum out-of-eclipse flux. 
On the finding chart, the image of this star is slightly elon- 
gated in the EW direction, suggesting a blend with another, 
fainter star which would lie 1" or slightly farther away to the 
West. Nevertheless, no clear sign of a 3rd light is seen in the 
lightcurve. 

The RV curves are well constrained with 1 1 out-of-eclipse 
spectra and notably observations close to phase 0.75. This sys- 
tem was previously studied by HHH05. There is significant dif- 
ferences between their RV parameters and ours. Our RV semi- 
amplitudes are 135 and 259 km s"', to be compared to their 
values of 160 and 247 km s"'. Beside having lower S /N and 
resolving power than us, in this particular case the discrepancy 
is certainly due to their admitted lack of observations close to 
the quadratures. Consequently, our value for the mass ratio, 
q = 0.52, is certainly more secure than theirs (0.65). 

We found a spectroscopic B luminosity ratio of 1.45. This 
is higher than the photometric value (1.29), perhaps because 
of the large distortion of the Roche lobe filling companion. 
Interestingly, the brighter, i.e. primary, component has lower 
monochromatic luminosities than the secondary component: 
even though the primary has a higher bolometric luminosity, 
it emits mostly in the UV part of the spectrum, so that its / 
luminosity, for instance, is lower than for the secondary. 

The most interesting parts of the disentangled spectra of 
both components are presented in Fig. |4] As a consequence of 
the low B luminosity of the primary, the spectrum of the latter is 
the noisier of the pair. Not surprisingly, in both spectra the most 
prominent features are the Hi and He: lines. It is tempting to 
identify a number of features in the primary spectrum with the 
Cii 4267, On 4276-4277, Sim 4553, Siiv 4089 and Siiv 4116 



lines. Nevertheless, both the lack of positive identification of 
the He ii 4542 line for a ~ 14 Mq star and the noisy profile of the 
He I lines mean that one must be careful in not over-interpreting 
a spectrum of rather low quality. The better secondary spectrum 
displays cleaner features. The Hei 4471 and Mgii 4481 lines 
allow to secure the temperature of the secondary. 

By fixing the photometric temperature and B luminosity ra- 
tios, a least-squares fit of the 1 1 out-of-eclipse spectra provided 
a primary temperature very close to 32 500 K, that is to say 
7000 K more than what was determined by HHH05. 

Both mass-log g and the HR diagrams are typical of a mas- 
sive Algol-type binary. The brighter and more massive compo- 
nent of the system appears to be close to the zero-age main se- 
quence (ZAMS), while the secondary component is larger and 
far more luminous than a non-evolved star of the same mass. 

4.2. 4 113853 

The best fit was obtained with a semi-detached model. Because 
of a moderate RMS scatter (0.017) combined with shallow 
eclipses (~0.18 mag) if not purely ellipsoidal variations, the 
light curve of this binary is one of the poorest of the whole 
sample. This low amplitude is due to a low inclination (~60°). 
The O -C curve reveals that the profiles of the eclipses are not 
perfectly symmetrical. The quality of the data is not sufficient, 
though, to trace the possible astrophysical cause of this asym- 
metry. On the finding chart, the star seems fairly well isolated. 

Despite only 7 out-of-eclipse spectra, the RV curves are 
rather well constrained with observations close to both quadra- 
tures. 

The S /N of the composite spectra are low (25 - 65) 
and there is a sizeable nebular emission in the Balmer lines. 
Because of the lack of metallic lines and the severe contamina- 
tion of the Balmer lines by nebular emission which hinder the 
disentangling procedure, the disentangled spectra of the com- 
ponents were not used. The least-squares fit was performed, 
letting both temperatures and the B luminosity ratio free to 
converge. It provided a temperature ratio remarkably close to 
the photometric one, and a spectroscopic B luminosity ratio 
of 0.68, also in perfect agreement with the photometric ratio 
(0.68). Thus, fixing the temperature and luminosity ratios to 
the photometric values was not needed to determine a reliable 
temperature of the primary. 

Both the mass-log g and the HR diagrams show an evolved 
system with a primary component seemingly half-way between 
the ZAMS and the terminal-age main sequence (TAMS). On 
the HR diagram, the primary lies much higher than the evolu- 
tionary track corresponding to its mass. Whether this is due to 
a temperature overestimate (linked e.g. with an underestimated 
sky background) or to some evolutionary effect remains to be 
examined. Besides, the distance modulus perfectly agrees with 
the currently accepted value for the SMC. 

4.3. 4 117831 

This faint system has a low-to-medium quality light curve of 
the EA type. There is a slight ellipsoidal variation between the 
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Fig. 4. Sections of the disentangled spectra of the primary and secondary components of the binary 4 1 10409. The scale on the 
vertical axis is adapted to the size of the features in each section. The red line represents the best fit synthetic spectrum shifted to 
the systemic velocity. Beside prominent H i and He i lines, Mg ii and Si in lines are visible in the spectrum of the secondary. The 
spectrum of the primary is far less convincing; nevertheless Si iv lines seem to be present next to H(5. 



eclipses and the latter have a similar depth (~0.4 mag). This is 
a close detached system with similar components. The finding 
chart suggests a possible slight blend with a fainter star located 
some 1" to the East of the system. No clear sign of a 3rd light, 
however, is seen in the lightcurves. 

The RV curves are well constrained with 12 out-of-eclipse 
spectra and observations close to phase 0.25 and phase 0.75. 
The mass ratio close to one (q - 0.98) is indicative of a binary 
with 'twin' components. 

The disentangled spectra are of too low quality to see any 
useful metallic line. The Mgii 4481 and Cii 4267 are barely 
visible. Disentangling of the Balmer lines was hindered by the 
strong emission. A first least-squares fit provided both tempera- 
tures and a poorly constrained spectroscopic B luminosity ratio 
of 0.97 ± 0.1 1. The WD code converged to a higher luminos- 
ity ratio (1.17 + 0.05), even when we tried to minimize it by 
fixing the potential of the primary. The temperature of the pri- 
mary was finally set by a fit where the ratio of temperatures 
was fixed to the photometric value, and the luminosity ratio 
assumed equal to one. The small number of photometric data 
in the minima, especially the primary one, probably makes the 
photometric luminosity ratio unreliable and explains why the 
radius of the primary component appears slightly smaller than 
that of the secondary one. 

According to the mass-log g diagram, the age of the system 
is about 50 Myr, assuming the standard SMC metallicity Z = 
0.004. The positions of both components in the HR diagram 
agree to within the error bars with the evolutionary tracks. 

This system was studied by Wyithe et al. 120021 (see Table 
|5]l. Their results not being constrained by spectroscopy, it is not 
surprising that they found a very different solution. They con- 



sidered this system as a semi-detached binary with a photomet- 
ric mass ratio of 0.157. Our spectroscopic results completely 
rule out that model. 



4A. 4 121084 

This system displays deep eclipses (> 0.6 mag) of similar 
depth. A slight ellipsoidal variation is visible. This is clearly a 
close detached system with slightly distorted twin components. 
No clear sign of crowding is seen on the finding chart, except 
possibly with very faint neighbour stars. 

The RV curves are well constrained with 9 out-of-eclipse 
observations regularly distributed around the quadratures. 

The composite spectra are polluted by strong nebular emis- 
sion in both Hy and H6 lines. Nevertheless, the widely sepa- 
rated Balmer lines allow a reliable temperature and luminos- 
ity ratio determination. The disentangled spectra are useful to 
confirm the rather high Vrot sin / values of the components. Not 
surprisingly, no metallic lines are visible because of the moder- 
ate S /N combined with fast rotational velocities. The potential 
of the primary was fixed so that the luminosity ratio given by 
the WD code matches the spectroscopic one. The temperature 
of the primary was obtained by fixing the temperature ratio to 
the photometric one. 

Both stars lie on the ZAMS, both in the mass-log g and HR 
diagrams. On the HR diagram, however, they are clearly more 
luminous and hotter than their expected positions for a metal- 
licity Z = 0.004. They would better agree with the ZAMS and 
evolutionary tracks for Z - 0.001, as many other systems do. 
Moving the representative points to the their expected positions 
for Z = 0.004 would require a 2000 K decrease in effective 
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temperature; that seems large, but the residuals between the ob- 
served and synthetic composite spectra show only very subtle 
changes. Only a modest systematic effect might be responsible. 



-0.11 



4.5. 4 121110 

The medium-to-high quality light curve shows a deep (~0.5 
mag) primary ecUpse. A slight ellipsoidal variation is visible 
between the eclipses. This is again a close detached system 
with slightly distorted components. No star closer than 3" is 
seen on the finding chart, except for a very faint one lying about 
2" away to the SW. 

The RV curves are well constrained with 1 1 out-of-eclipse 
spectra. 

There is strong nebular emission in both Balmer lines. The 
spectroscopic B luminosity ratio (0.415 + 0.047) nicely agrees 
with the photometric one (0.424), without any need for fixing 
the potential of the primary. The temperature of the primary 
was fitted after fixing the temperature and luminosity ratios 
to their photometric values, as ususal. The Sim 4553 line is 
clearly visible on the disentangled spectrum of the primary. The 
lack of Mg 11 4471 confirms the relatively high temperature of 
the primary. The spectrum of the secondary is too noisy for the 
identification of metallic lines. 

On the mass-log g diagram, the stars match an isochrone 
corresponding to about 7-8 Myr In the HR diagram, the po- 
sitions of both components are above the Z = 0.004 evolution- 
ary tracks but are consistent with the lower metallicity ones 
(Z = 0.001). Increasing the helium content would also help to 
reconcile their positions with the evolutionary tracks, unless a 
systematic effect raises the apparent effective temperatures. 

4.6. 4 121461 

This is an eccentric system with two (relatively) widely sepa- 
rated components. Both eclipses are very similar in depth and 
width. With P ~ 17.9 mag, this is one of the faintest sys- 
tems in our sample. Nevertheless, the finding chart indicates 
no crowding problem whatsoever No significant apsidal mo- 
tion was found on the basis of photometry. An analysis with the 
EBOP code shows that the loq value depends critically on the 
e sin Wo quantity, which is poorly constrained, while the more 
robust e cos loq quantity is such that cos loq ^ 0.75 and thus 
does not constrain wo very tightly. Fig. |5] suggests a marginal 
decrease of e cos wq with time which, if real, could only be due 
to gravitational perturbations from a third body, because d) <Q, 
while pure tidal effects always result in <ii > 0. We have as- 
sumed no apsidal motion. 

This system is one of the two in our sample having 16 out- 
of-eclipse spectra, of which 15 were used (the 8th one has too 
poor SNR). There is a sufficient number of observations close 
to the quadratures. 

The composite spectra are very noisy. The disentangled 
spectra are of very low quality, with no exploitable metallic 
lines. There is strong nebular emission. Letting the temperature 
of both components free to converge (together with the B lumi- 
nosity ratio) resulted in very uncertain values, so we fixed the 
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Fig. 5. Variation with time of the e cos co quantity of the sys- 
tem 4 121461, determined with the EBOP code on the basis 
of the / magnitudes grouped in four time intervals. The solid 
line corresponds to a hardly significant negative apsidal mo- 
tion, while the dashed horizontal line indicates no apsidal mo- 
tion, fore = 0.185 and wq = 222.8°. 



temperature ratio to the photometric one in order to estimate 
the temperature of the primary. The rotational velocities were 
fixed to the pseudosynchronized values. A luminosity ratio of 
0.95 + 0.06 was found on the basis of the spectra, which agrees 
well with the photometric one (0.91 + 0.03) obtained without 
fixing the potential of the primary. The photometric luminosity 
ratio was adopted, which results in almost identical radii for the 
components. This results in a slightly lower surface gravity on 
the secondary component than on the primary, because of the 
mass ratio, but this difference is not significant. 

On the mass-log g diagram, this system lies close to the 
ZAMS but might be up to 15 - 20 Myr old. The positions of 
the stars in the HR diagram agree, within the error bars, with 
the evolutionary tracks, although they tend to lie too high, as is 
the case of other systems. 

4.7. 4 159928 

This system presents low-to-medium quality light curves of 
the EW type with eclipses of unequal depth. The RMS scat- 
ter of the light curve is rather low, but the minima are not very 
deep (A/mini ~ 0.25 mag). It can be inferred from this that the 
system comprises close, strongly distorted components of un- 
equal brightness. The best-fitting model corresponds to a semi- 
detached binary with a low inclination, close to 60°. The find- 
ing chart shows a well isolated target, except for a quite faint 
neighbour at about 2" to the NNE. 
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Fig. 6. Example of observed spectra close to a quadrature. Both spectra are at the same scale. The red line represents the com- 
posite synthetic spectrum of the binary system, i.e. the addition of the scaled and velocity-shifted synthetic spectra of the two 
components for the corresponding orbital phase. Beside a low S /N, the spectrum of 4 121461 shows strong nebular emission in 
the Balmer lines. From the best fit, the B luminosity ratio was found to be ~0.95 with a primary temperature of ~22 000 K. The 
spectrum of 5 123390 points to a B luminosity ratio of 0.59 and a primary temperature of 27 000 K. 



There are only 8 out-of-eclipse spectra, but they are close 
to the two quadratures. 

A spectroscopic luminosity ratio of 0.45 was found in the 
B band, which is higher than the photometric value (0.36). The 
temperature of the primary was determined by fixing the tem- 
perature and luminosity ratios to the photometric values, after 
suppression of the very strong nebular emission lines. The Si iii 
4553 line is the only metallic line clearly visible in the disentan- 
gled spectrum of the primary. The spectrum of the secondary 
shows no exploitable metallic line. The Mgii 4481 is barely 
visible. The synchronized values for y,ot sin / are close to 200 
km s ' and therefore all but the strongest lines are buried in the 
noise. 

Both the mass-log g and HR diagrams show positions typi- 
cal of a semi-detached system with an evolved secondary com- 
ponent. The primary is slightly overluminous relative to the 
evolutionary track of a single star. 



of the primary was fixed to a value that implies a luminosity ra- 
tio close to the spectroscopic one. No metallic lines are visible 
in the very noisy disentangled spectra. The nebular emission is 
strong in both H i lines. 

Despite the moderate quality of the photometric and spec- 
troscopic data, the positions of both stars fall right on the 
ZAMS in the mass-log g diagram. In the HR diagram, their 
position agree well with the evolutionary tracks, though they 
appear slightly overluminous. 

A notable characteristic of this system is its fast apsidal mo- 
tion (1> - 9.8 ± 1.9 ° yr"'. Figure |7] shows the ecos w product 
as a function of time, as obtained using the EBOP code. The 
solid line represents the WD solution, which appears consis- 
tent with the EBOP results, even though the latter would be 
compatible with a faster apsidal motion coupled with a slighly 
smaller eccentricity. Further discussion of this result is deferred 
to Section 5. 



4.8. 4 160094 

This detached system of moderate eccentricity presents low- 
to-medium quality light curves with rather shallow eclipses 
(A/mini ~ 0.20 mag). Except for a few very faint neighbours, 
the target seems free from crowding on the finding chart. 

There are 1 1 out-of-eclipse spectra. Both quadratures are 
well covered by the observations. 

The Teff of the primary was determined together with that of 
the secondary and with the B luminosity ratio. The temperature 
of the secondary proved rather ill-defined, so the photometric 
temperature ratio was used to define it, as usual. The potential 



4.9. 4 163552 

This is one of the five systems with /^ < 16 mag, displaying 
a high quality light curve of the EB type. The eclipses have 
very similar depths (~ 0.4 mag), indicating that the tempera- 
ture ratio is very close to unity. Some faint neighbours are seen 
on the finding chart about 2" from the target, and it is difficult 
to judge whether closer neighbours might lurk within the rel- 
atively large spot left by this bright system. According to the 
light curve, there is a substantial third-light contribution to this 
system: we find ijj ~ 0.12, while Graczyk (.2003. hereafter 
dG03) found ^3/ ~ 0.26. A significant ellipsoidal variation is 
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Fig. 7. Variation with time of the ecosw quantity of the 4 
160094 system, determined with the EBOP code on the basis 
of the / magnitudes grouped in four time intervals. The solid 
line corresponds to the parameters listed in Table [TTI 



indicative of tidally distorted components, though the system is 
still in a detached configuration. 

There are only 9 out-of-eclipse spectra, but the RV curves 
are well constrained near quadratures. 

The observed spectra of this binary were roughly corrected 
for the presence of a third component, by subtracting a constant 
from the normalized spectra. The value of the constant was ad- 
justed until the efifective temperatures obtained from the least- 
squares fit to the corrected spectra give a reasonable match to 
the evolutionary tracks at the metallicity Z = 0.004 in the HR 
diagram. A constant of 0. 1 only proved sufficient to this pur- 
pose, which is widely dififerent from dG03's photometric es- 
timate of the third light, {t,^b - 0.26. On the other hand, our 
WD analysis resulted in l-i^B - 0.14, which is more consistent 
with our rough spectroscopic estimate. Actually, a first attempt 
to analyze this system resulted in highly overluminous com- 
ponents. This was interpreted as a clear sign that the observed 
spectra were severely contaminated by the object responsible 
for the photometric third light. 

The disentangled spectrum of the primary shows a faint 
Mgii 4481 line and a noisy Sim 4553 line. These lines were 
used to constrain the temperature of the primary. The Si in 4553 
line is visible in the spectrum of the secondary too. 

This system was studied by HHH05, dG03 (see his Fig. 
6) and WW02. Only dG03 took the third-Hght contribution 
into account. Not surprisingly, our photometric solution is very 
close to theirs. 

Since we adjusted the amount of third light so that the pri- 
mary component have a realistic luminosity, the HR diagram 



shows a primary close to the stellar evolutionary track of a 9.6 
Mq star, and a secondary close to the evolutionary track of a 
9.1 Mq star Their positions are clearly far from the zero-age 
main sequence. From the mass-log g diagram, the age of the 
system is estimated to be 20 Myr As the position of the pri- 
mary was inferred from the expected position in the HR dia- 
gram, the proposed solution is not entirely reliable and thus is 
not used for the determination of the mean distance modulus of 
the SMC (see SectionlS!6ll. 

4.10. 4 175149 

The medium-to-high quality light curve of this system is of the 
EB type; unfortunately, most of the right part of each eclipse 
is missing. That is due to a period very close to 2 days ex- 
actly. The minima are fairly deep (A/njini ~ 0.50 mag), well- 
defined and of unequal depth. There is a strong ellipsoidal vari- 
ation. The binary is a semi-detached system, with distorted 
components of different brightnesses. A slight depression oc- 
curring before the primary eclipse is indicative of a gas stream. 
However, we did not venture into adding a spot on the primary 
in order to mimick this effect, because we felt that the large 
gaps in the light curve already limit the reliability of the pro- 
posed solution. 

This system is close to the edge of the CCD used in the 
OGLE 11 survey, so it was also listed under the designation 
5 38079 in the OGLE database. Thanks to that circumstance, 
there are 574 photometric / magnitudes instead of less than 
300. An examination of the OGLE finding chart reveals that 
this binary is clearly blended. 

This is the second system with 16 out-of-eclipse spectro- 
scopic observations. Because of the 2-day period, all these ob- 
servations took place before the quadratures. Nevertheless, the 
RV curves seem to be sufficiently constrained. 

A spectroscopic B luminosity ratio of 1.39 was found, 
though with a large scatter (cr ~ 0.14). This is higher than the 
photometric value (1.18). The disentangled spectra are fairly 
good. The spectrum of the hot primary component displays 
the He ii 4200 and 4542 lines. The Si m 4553 and a faint Si iv 
4089 lines are equally visible. These lines provide strong con- 
straints for the temperature of the primary. The best-fitting pri- 
mary temperature was obtained after fixing the temperature and 
luminosity ratios to the photometric values. 

The mass-log g diagram shows the typical oblique orienta- 
tion of the segment connecting the components of an evolved 
binary. The HR diagram shows a highly overluminous primary, 
relative to the evolutionary track of a 11.8 Mo star, and an 
evolved secondary far more luminous than a main-sequence 
star of 7.8 Al©. One may wonder, whether the strange posi- 
tion of the primary is due to an unrecognized third light, both 
in photometry and spectroscopy, or to some evolutionary efifect. 
In any case, there is no obvious third light in the light curve. 

4.11. 4 175333 

This slightly eccentric system presents low-to-medium quality 
light curves of the EA type. This is one of the few systems with 



20 



P. North et al.: VLT spectroscopy of eclipsing binaries in the SMC 



Z'' > 17.5 mag. Consequently, the O - C curve shows a rela- 
tively high scatter. The minima are of unequal depth. A slight 
ellipsoidal variation is visible. This is clearly a detached sys- 
tem with components of unequal brightness. The target appears 
perfectly isolated on the finding chart. 

There are 14 out-of-eclipse spectra. The observations con- 
strain well the RV curves. 

This system was studied by WWOl. The fact that they con- 
sider the eclipses as total (while we consider them as partial) 
and their lack of spectroscopic constraints on their ratio of radii 
account for the differences between their solution and ours. 
The evidence for total eclipses does not appear compelling, 
so additional photometry would be needed to settle the issue. 
We found a spectroscopic B luminosity ratio of 0.55, slightly 
smaller than that finally adopted taken photometry into ac- 
count. The disentangled spectrum of the primary shows a noisy 
Mgii 4481 line. The spectrum of the secondary is too noisy to 
detect any metallic line. The temperature of the primary was 
fitted with the temperature and luminosity ratios fixed to the 
photometric values. The fit with both temperatures free to con- 
verge, together with the luminosity ratio, gave a rather large 
scatter of about 1600 K and a secondary temperature about 
900 K cooler. In spite of the partial eclipses, there was no need 
to fix the potential of the primary in order to find a B luminosity 
ratio that matches the spectroscopic value, so the photometric 
value of the luminosity ratio was adopted. 

On the mass-log^ diagram, both stars fall right on the 20 
Myr isochrone. On the HR diagram, however, both stars ap- 
pear significantly overluminous relative to their evolutionary 
tracks, suggesting that the effective temperatures may be over- 
estimated by at least 2000 K! Strangely enough, the colour ex- 
cess of this system appears to be small (E(B - V) = 0.07) and 
the distance modulus (18.6) clearly smaller than the accepted 
value for the SMC (~ 18.9 - 19.0). If the effective temperatures 
had indeed been overestimated, this would have implied both 
a too blue intrinsic colour and a too large intrinsic luminos- 
ity (hence a more negative absolute magnitude), se one would 
rather expect a large colour excess and a large distance modu- 
lus. 

4.12. 5 016658 

This close detached system presents medium quality light 
curves with eclipses of equal depths, and is composed of 
tidally distorted twin components. The finding charts reveals 
no crowding problem. 

There are 1 1 out-of-eclipse spectra. The RV curves are well 
constrained by the observations around phase 0.25. 

This system was studied by WWOl. As for the previous 
binary, the differences observed between their (photometric) 
solution and ours is due to their lack of a spectroscopic con- 
straint on the ratio of radii, and to the fact that they assume 
total eclipses. Evidence for the latter is not compelling, how- 
ever, and awaits further photometric measurements for confir- 
mation. A spectroscopic B luminosity ratio of ~0.60 was found, 
which guided the choice of the potential of the primary compo- 
nent in the WD analysis. The Mgii 4481 line is clearly visible 



on the disentangled spectrum of the primary. The spectrum of 
the secondary is too noisy to show any metallic line. The best- 
fitting primary temperature was obtained simultaneously with 
the temperature of the secondary, which appeared quite com- 
patible (within 300 K) with the photometric one (i.e. given the 
spectroscopic primary temperature and the photometric tem- 
perature ratio), and with the luminosity ratio. The photometric 
temperature ratio was adopted. 

On the mass-log g diagrams, both stars fall on the 30 Myr 
isochrone within the errors. On the HR diagram, the primary 
has a position compatible with its evolutionary track within er- 
rors, though it appears slightly too luminous. The secondary 
is slightly hotter than the primary, and is more overluminous; 
still, it remains compatible with its evolutionary track if the er- 
rors on both luminosity and mass are considered. 

4.13. 5 026631 

This system presents a medium-quality light curve of the EW 
type with minima of unequal depth. This is clearly a semi- 
detached system with strongly distorted components of un- 
equal brightness. It presents the second lowest inclination of 
the sample with / » 61°, implied by the rather small amplitude 
of the light curve. No blend is apparent on the finding chart, 
except for two or three very faint neighbours at about 2". 

There are only 8 out-of-eclipse spectra, but these ones are 
sufficiently constraining to get reliable RV curves. 

A spectroscopic B luminosity ratio of 0.74 was found, 
while the photometric value is 0.50. The disentangled spectra 
show no useful metallic lines. This is due notably to the high 
Vrot sin; values (~190 km s"'). The temperature of the primary 
was obtained by fixing the temperature and luminosity ratios to 
the photometric values. 

The mass-log^ and HR diagrams show the typical posi- 
tions for the components of a semi-detached system, with the 
primary near its expected evolutionary track and an overlumi- 
nous secondary. The primary is slightly overluminous relative 
to its track, as is often the case in this work, while the secondary 
is slightly below its track, a rare occurrence. 

This binary was studied by HHH05. Their primary temper- 
ature (25 500 K) and mass ratio (~1) differ noticeably from our 
values. 



4.14. 5 032412 

This wide, detached system presents medium-to-high quality 
Ught curves with minima of unequal depth, betraying compo- 
nents of unequal brightness. The target appears well isolated on 
the finding chart. 

There are 13 out-of-eclipse spectra. The RV curves are very 
well constrained and the RMS scatters are low. Interestingly, 
both the light and velocity curves indicate a negligible eccen- 
tricity, in spite of the small relative radii of the components, as 
if circularization had taken place during the protostellar phase. 
Note that this is the most massive system of our whole sample: 
its total mass reaches 30 Mq. 



p. North et al.: VLT spectroscopy of eclipsing binaries in tiie SMC 



21 



A spectroscopic B luminosity ratio of 0.55 was found. The 
disentangled spectra are of high quality (Fig. [8]l, even for the 
Hi lines. Beside the Hi and Hei lines, the following lines are 
visible in the spectrum of the primary: He ii 4200 and 4542 
(strong), Si iv 4089, Si iv 41 16, O n 4185, Si iv 4212, O ii 4276- 
7 and Si m 4553. An effective temperature of about 35 000 K 
was inferred from the best-fitting synthetic spectrum. 

The following metallic lines are visible in the disentangled 
spectrum of the secondary: Siiv 4089, Siiv 4116, On 4185, 
On 4190, Cii 4267, On 4276-7, On 4415-7, Mgii 4481 and 
Si m 4553. He ii 4542 is clearly visible too. Comparing the rel- 
ative depths of Mg ii 4481 with He 1 4472, He ii 4542 with Si iii 
4553, C II 4267 with O ii 4276-7, Si iv 4089 and Si iv 41 16 with 
He I 4121 allows to estimate an effective temperature close to 
31 000 K. Thanks to the good SNR of the spectra, fitting si- 
multaneously the temperatures of the components and the lu- 
minosity ratio resulted in a temperature ratio very close to the 
photometric one. Nevertheless, the adopted temperatures are 
those obtained by imposing the photometric ratio. 

The mass-log g diagram shows a very young binary with 
both components on the ZAMS. On the HR diagram, the po- 
sitions of both components agree fairly well with the stellar 
evolutionary tracks of 17.1 and 13.1 Me stars. However, the 
primary appears slightly overluminous relative to its track. 

4.15. 5 038089 

This bright detached system presents medium-quality light 
curves with eclipses of very similar depth. Therefore, the com- 
ponents are very similar. 

This system has abnormal colour indices (see Table [3] and 
Fig. [3J, which suggests the presence of a third light. An ex- 
amination of the OGLE finding chart reveals that this binary is 
clearly blended, though with much fainter stars. 

There are 1 1 out-of-eclipse spectra. The RV curves are very 
well constrained and the RMS scatter is remarkably low. 

A spectroscopic B luminosity ratio of 0.68 was found. The 
eclipses are not very deep, so we fixed the potential of the pri- 
mary to reproduce this luminosity ratio. The disentangled spec- 
tra, shown in Fig. |9] are very similar. Strong Heii 4200 and 
4552 lines are visible in both spectra. The following metallic 
lines are equally identifiable (Fig.©: O n 4076, Si iv 4089, Si iv 
41 16, O II 4185, O II 4190, C II 4267, O ii 4276,0 ii 4415-4417, 
Sim 4553. This wealth of lines allows to determine the tem- 
peratures of both components with a great accuracy. From the 
best-fitting synthetic spectra, we found 30400 K and 30 800 
K for the eff'ective temperature of the primary and secondary, 
respectively. This is very close (i.e. within 200 K) to the tem- 
peratures estimated from the composite spectra by imposing 
the photometric temperature and luminosity ratios; thus, one 
can safely conclude from this example that the two methods 
are equivalent. 

The high quality of the spectroscopic observations allowed 
to estimate the astrophysical parameters with a greater accu- 
racy than most systems in our sample. On the mass-log g di- 
agram, both components lie just above the 10 Myr isochrone; 
their respective positions suggest that the ratio of radii may be 



slightly underestimated. The HR diagram shows that both stars 
are significantly overluminous with respect to the evolution- 
ary tracks of 13.0 and 11.7 Al© stars. Invoking the blending 
of the binary with a third-light contributor does not seem to 
help much. No clear sign of a third light can be seen in the 
lightcurve; this is admittedly a weak argument, since the light 
curve of a detached system with weak proximity effects cannot 
constrain well a third light. But, in addition, a third light will 
not change much the relative radii of the components, and the 
eff'ective temperatures seem well constrained by the relative in- 
tensities of several lines, so that the luminosities should remain 
unaff'ected. Furthermore, the radial velocity curves are of such 
quality that it is difficult to imagine how the masses could be bi- 
ased otherwise than through the inclination angle /. A third light 
would make the photometric minima less deep and so, it would 
mimick a lower /. Since the RV amplitudes give the product 
M sin^ /, underestimating / is equivalent to overestimating M, 
while we would need the reverse to explain the HR diagram we 
see. A test with the EBOP code confirms this qualitative argu- 
ment: assuming a third light f^j = 0.2 changes the radius of the 
primary by -0.26% only (keeping the ratio of radii constant), 
so that the luminosity decreases by half a percent, while the in- 
clination is increased by more than three degrees, lowering the 
masses by 3.6%. That would imply evolutionary tracks with a 
luminosity about 0.05 dex lower (or 11%) in the HR diagram. 
Finally, the distance modulus of this system is very close to the 
one expected for the SMC. Therefore, for a third light to be the 
cause of the inconsistency, it can only be via the RV curves, 
the amplitude of which should be biased to small values by a 
stationary third spectrum. 

This binary was part of the first release of 10 bright SMC 
systems by HHH03. Comparing our results with theirs, we see 
that we have similar estimates for the mass ratio and tempera- 
tures of the stars, but markedly different values for the masses 
and radius of the secondary. A photometric solution was pro- 
posed by dG03, who do not recommend it for distance deter- 
mination because its components do not fall on their mass- 
luminosity relation. 

4.16. 5 095337 

This is a close but detached system, with tidally distorted com- 
ponents of unequal brightness. The O - C curve suggests that 
the primary eclipse is not perfectly reproduced by the synthetic 
light curve, but there is no obvious third-light contribution. The 
finding chart shows some blends with two or three fainter stars 
at r.'5 or so. 

There are 10 out-of-eclipse spectra. The RV curves are well 
constrained with observations close to both quadratures. 

A spectroscopic B luminosity ratio of 0.66 was found. The 
disentangled spectra are of low quality. Because of the low 
S /N and the high Viot sin / (~200 km s"'), no metallic line can 
be positively identified. There is some nebular emission in the 
Balmer lines. The best-fitting synthetic spectra allowed to esti- 
mate simultaneously the temperatures of the primary and of the 
secondary, the ratio of which agrees quite well with the photo- 
metric one. In addition, the spectroscopic B luminosity ratio 
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Fig. 8. Sections of the disentangled spectra of the primary and secondary components of the binary 5 032412, The scale on the 
vertical axis is adapted to the size of the features in each section. The red line represents a velocity-shifted synthetic spectrum 
of the star. Beside prominent H i and He i lines. He ii, O n, Mg ii, Si iii and Si iv lines are visible in both spectra. The strong He ii 
4542, next to a smaller Sim 4553 line, and the Hen 4200 line confirm the high temperature of the primary (~34900 K), The 
lower temperature of the secondary is confirmed by the stronger Mg ii 448 1 line and the smaller He ii lines. 



perfectly matches the photometric one, so that the potential of 
both the primary and the secondary were left free to converge. 
On the mass-log g diagram, both components define a seg- 
ment which is perfectly parallel and very close to the 10 Myr 
isochrone. On the HR diagram, however, both components 
appear strongly overluminous compared to the evolutionary 
tracks of 8,7 and 7,6 Al© stars. Decreasing the effective tem- 
peratures by about 1300 K would reconcile the luminosities 
with the tracks. However, this looks difficult. The emission in 
both Balmer lines was suppressed on a 4 A range centered on 
each emission line, and we verified that increasing that range 
to 8 A does not change the estimated temperature in a signif- 
icant way. Thus, either this system suffers from some bias on 
the RV curves, or its metallicity is closer to Z = 0.001 than to 
Z = 0.004. 

4.17. 5 095557 

This is the system with the highest eccentricity, displaying a 
medium quality light curve with minima of unequal depths. 
The target is perfectly isolated on the finding chart. 

There are 1 1 out-of-eclipse spectra. The RV curves are well 
constrained with observations close to both quadratures, but the 
fit is not very good and, unfortunately, most spectra are grouped 
in the phase interval with the smaller amplitude. 

A spectroscopic B luminosity ratio of ~0.5 was found when 
limiting the fit to the seven spectra for which the radial velocity 
difference ARV > 250 kms"'. This ratio increases to 0.63 if 
all eleven spectra are taken into account, so that this quantity is 
rather poorly constrained. The disentangled spectra are of very 



low quality. This is probably due partly to the low S /N of the 
observed spectra, and partly to some inaccuracies in the orbital 
parameters. The temperature of the secondary given by the fits 
to the composite spectra is about 1000 K higher than the photo- 
metric estimate, which was adopted. The pseudosynchronized 
values of Viot sin / were adopted. 

An apsidal motion is detected at the 5 cr significance level, 
and the WD result is confirmed by the variation of the e cos a> 
quantity as given by the EBOP code. The photometric data was 
divided into four sets and the fits were obtained by fixing the 
inclination, the ratio of radii and the relative radius of the pri- 
mary. The result is displayed in Fig.[10]and suggests the apsidal 
motion to be real. The adopted apsidal motion seems underes- 
timated on that figure, but the constraint imposed by the RV 
curves has to be kept in mind. On the mass-log g diagram, the 
two components lie right on the 30 Myr isochrone, in spite of 
the rather large uncertainty on the masses. On the HR diagram, 
both components appear clearly overluminous relative to their 
respective evolutionary tracks, as in many other systems. 

4.18. 5 100485 

This is a detached system with 'twin' components and a circu- 
lar orbit. The finding chart reveals no crowding problem, the 
closest neighbour being at about 2'.'5. 

There are 13 out-of-eclipse spectra. The RV curves are well 
constrained, with observations close to both quadratures and 
small residuals, 

A spectroscopic B luminosity ratio of 0,93 was found, 
which matches the photometric value (0.97) rather well. The 
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Fig. 9. Sections of the disentangled spectra of the primary and secondary components of the binary 5 038089. The scale on the 
vertical axis is adapted to the size of the features in each section. The red line represents a velocity-shifted synthetic spectrum of 
the star. Beside prominent H i and He i lines. He ii and a number of fainter metallic lines (C ii, O ii, Si iii, Si iv) are visible. The 
relative intensities of the He ii 4542 and Si iii 4553 are very useful to constrain the temperature of the two stars. 



disentangled spectra are very similar. The following metallic 
lines can positively be identified: Cii 4267, Mgii 4481 and 
Si m 4553. These lines allowed to estimate a temperature close 
to ~22000 K for each component. The least-squares fit, per- 
formed with both temperatures and the luminosity ratio free to 
converge, provided a temperature ratio very close to the photo- 
metric one. 

The two components lie right on the 20 Myr isochrone on 
the mass-log^ diagram, and they are compatible, within the 
errors, with their respective evolutionary tracks in the HR di- 
agram. Here again, however, they are slightly overluminous, 
unless their metallicity is low. 

4.19. 5 100731 

This system presents low-to-medium quality light curves of the 
EW type, with the smallest amplitudes among those in our sam- 
ple. This is another case of binary seen under an unfavourable 
inclination (~60°). The finding chart shows a perfectly isolated 
target. 

The most satisfying fit of the light curves was obtained with 
an overcontact model. 

There are only 8 out-of-eclipse spectra. Nevertheless, these 
observations are sufficiently constraining to get reliable RV 
curves. 

Spectroscopy gives a B luminosity ratio of ~ 0.44, while 
the photometric ratio amounts to ~ 0.51. As usual, the lat- 
ter was preferred, especially because of the strong constraints 
provided by a Roche-lobe-filling pair. The disentangled spectra 
have a low S /N, and a least-squares fit to the composite spectra 
with both temperatures free provided two temperatures close 



to 23 000 K. However, the scatter of the secondary temperature 
was high, so the photometric temperature ratio was adopted and 
fixed. 

The mass-log g diagram shows a secondary component 
more evolved than the primary, as expected. On the HR dia- 
gram, both components appear to be underluminous with re- 
spect to their evolutionary tracks, though by a little amount. 
This is an exceptional occurrence. 

4.20. 5 106039 

This system is a typical semi-detached one. There is a small de- 
pression occurring just before the eclipse of the primary, which 
is strong evidence for a gas stream. We did not attempt to model 
that stream with a cool spot on the primary, because the dis- 
torsion of the light curve remains relatively mild. The target 
appears perfectly isolated on the finding chart. 

There are 9 out-of-eclipse spectra. The RV curves are well 
constrained with observations close to both quadratures. 

A spectroscopic B luminosity ratio of 1.03 was found, in 
excellent agreement with the photometric value (1.01). The fol- 
lowing metallic lines are visible in the disentangled spectrum 
of the primary: C ii 4267, Mg ii 448 1 (faint) and Si iii 4553. A 
stronger Mg ii 448 1 line is equally visible in the spectrum of the 
secondary. From these lines, the primary temperature was esti- 
mated to be close to 25 500 K. The least-squares fit, performed 
with temperature and luminosity ratios fixed to the photometric 
values, provided a primary temperature ~ 800 K higher. 

The position of the primary component in the HR diagram 
is in fair agreement with the theoretical evolutionary track of 
a 8.6 Al© star, though it is shghtly overluminous. The evolved 
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Fig. 10. Same as Fig.|7] but for the system 5 95557. The filled 
black dots represent independent e cos o) determinations made 
on four successive times series containing about the same num- 
ber of photometric observations in the / band. The solid line is 
based on the parameters listed in Table [TTl while the dotted 
line corresponds to the fitted w + 2cr. The horizontal dashed 
line shows the result obtained from the whole time series under 
the assumption of a constant co. 



secondary component is overluminous with respect to the track 
of a single star of the same mass. 

4.21. 5 111649 

This is a detached system with very slightly distorted twin com- 
ponents. There is a group of bright stars close to the target on 
the finding chart, but they are remote enough (3.5 - 4") that no 
3rd light should be expected from them. 

There are 10 out-of-eclipse spectra. The RV curves are not 
well sampled because the period of 2.95955 days is very close 
to an integer number of days, but the scatter is small, thanks 
to the slow projected rotational velocities induced by the rela- 
tively long period. 

A spectroscopic B luminosity ratio of 0.89 was found, 
which was imposed (via the potential of the primary) to define 
the ratio of radii. Indeed, the small amplitude of the light curve 
prevents a purely photometric ratio to be well constrained. 
Noisy C II 4267 and Mg ii 448 1 lines are visible in the disentan- 
gled spectrum of the secondary. The temperature of the primary 
was determined from a fit with fixed temperature and luminos- 
ity ratios. According to both spectroscopy and photometry, the 
secondary appears marginally hotter than the primary, the tem- 
perature difference being about 2 cr. 



Both components have a very similar mass (~5.4 Mq), and 
according to the mass-log g diagram, the empirical mass con- 
trast appears a bit too high to match the ~ 80 Myr isochrone. 
Interestingly, both components lie right on the Terminal Age 
Main Sequence (TAMS), where evolution is so fast that a the- 
oretical lower limit to the mass ratio can be settled. Starting 
from the purely empirical surface gravities, and increasing their 
difference by two sigma (so that loggp - 1 cr = 3.72 and 
loggs + I cr - 3.80), one can read the corresponding masses 
along the logf - 7.9 isochrone. Assuming that both stars have 
not yet passed the "hook" that marks the end of the mains se- 
quence, the resulting mass ratio is ^ = 0.98. But it is quite 
possible that the primary has just passed the hook while the 
secondary has not, in which case q - 0.966. In both cases, 
the mass ratio is closer to one than the value directly obtained 
from the RV curves by 0.03 to 0.04. Therefore, this system can 
be considered as hosting real twins. 

The position of the primary in the HR diagram agrees al- 
most perfectly with the theoretical evolutionary track, espe- 
cially if the star has just evolved beyond the "hook". That of 
the secondary, however, is a bit too high, as if its mass was 
underestimated. 



4.22. 5 123390 

This slightly eccentric system {e = 0.042) presents low-to- 
medium quality light curves with a small amplitude. This is 
clearly a detached system with components of unequal bright- 
ness. The finding chart shows a neighbour at about 1'.'5 to 
the SW of the binary, which might have polluted the spectra 
slightly. 

There are 14 out-of-eclipse spectra. The RV curves have a 
small RMS scatter and are very well constrained by the obser- 
vations. 

A spectroscopic B luminosity ratio of 0.58 was found. The 
disentangled spectra are of fairly good quality. The following 
metallic lines are visible in both spectra: On 4076, Cii 4267, 
On 4276, On 4415-4417, Mgii 4481 and Sim 4553. For the 
primary, the best-fitting temperature for these lines is ~26 000 
K. The least-squares fit of the composite spectra gave temper- 
atures of 28400 and 26260 K for the primary and secondary 
respectively, but with a large scatter A plot of the fitted temper- 
atures versus the unnormalized chi-square shows that in some 
cases, the fit switched components, i.e. attributed the high tem- 
perature to the secondary and vice-versa, which partly explains 
the large scatter (see Fig.fTTTi. Strangely enough, the photomet- 
ric ratio of temperatures is close to one, so that the tempera- 
tures become, when fitted while keeping this ratio fixed, 27840 
and 28320 K for the primary and secondary component respec- 
tively. Taken at face value, however. Fig. [TT] rather suggests 
29000 and 25000 K. 

According to the mass-log § diagram, this binary is ~12 
Myr old. The position of the primary component in the HR dia- 
gram is a bit too high relative to its evolutionary track. The sec- 
ondary component is much too luminous, falling on the track 
of the primary ! If the purely spectroscopic temperatures were 



p. North et al.: VLT spectroscopy of eclipsing binaries in tiie SMC 



25 



35000 




• 










30000 


— ■ 










— 




• 


• •• 




• 




_ 




• 














• 


• 








• _ 




• 

• 


• 


• 


• 




• 


25000 


— • 
• 


• 
• • 








- 


20000 










123390 





0.0003 0.0004 



0.0005 



0.0006 0.0007 



Fig. 11. Temperature of the primary (black dots) and of the sec- 
ondary (red dots) component of the system 5 123390, obtained 
by fitting composite synthetic spectra to the observed ones, ver- 
sus the quality of the fit. The trends are roughly horizontal, at 
least for good fits, which inspires confidence. In a few cases, 
the components are exchanged (see text). 
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Fig. 12. Same as Fig. |7] but for the system 5 123390. The filled 
black dots represent independent e cos a> determinations made 
on four successive times series containing about the same num- 
ber of photometric observations in the / band. The open red 
dots refer to e cos co determinations based on only two time se- 
ries, but which are twice longer The solid line is based on the 
parameters listed in Table [TT] 



adopted, the primary would be even more overluminous, but 
the secondary would fall right on its track. 

There is a marginally significant apsidal motion of 4.75 ± 
1.63 ° yr"'. We show in Fig.[T2]the e cos co quantity found with 
the EBOP code on the basis of the / magnitudes. There is in- 
deed a slight trend corresponding to an increase of w with time. 



4.23. 5 180185 

The EA-type light curves of this system are irregularly sam- 
pled, because the period is close to 5.5 days, and there are only 
very few data in the primary minimum. This is a typical well- 
detached system with twin components and a small eccentric- 
ity. The finding chart reveals a slightly fainter neighbour about 
2" away from the target, i.e. far enough that its influence on the 
spectra may be considered negligible. 

This system has unreliable colour indices (see Table[3]and 
Fig.0. 

There are 12 out-of-eclipse spectra. The RV curves are 
excellent and well constrained by the observations, with the 
smallest residuals seen in our sample. The formal errors on 
the resulting amplitudes are smaller than one percent, allow- 
ing mass estimates to within 2-3%. 

A spectroscopic B luminosity ratio of 0.63 was found, 
which is close to the value given by the WD analysis and fi- 
nally adopted (0.71). This is the system with the longest period 



(P ~ 5.5 days) in our sample. The long period is responsible 
for the low Vrot sin / values (~40 km s"') compatible with syn- 
chronous rotation. As a consequence, the disentangled spectra 
of this binary show rather sharp lines. The following metallic 
lines are visible in both spectra: C ii 4267, Mg ii 4481 and Sim 
4553. The temperature of the primary was obtained as usual by 
fixing the temperature and luminosity ratios to the photometric 
values, in spite of the small number of points in the minima. 
The secondary has a slightly hotter temperature than the pri- 
mary. 

Both components are fairly well aligned on the 50 Myr 
isochrone in the mass-log g diagram, although the primary 
should be slightly more evolved. On the HR diagram, the sec- 
ondary falls right on its track, while the luminosity of the pri- 
mary appears too low. This might be due to an unreliable tem- 
perature ratio, because of the small number of photometric 
points in the minima. 

4.24. 5 180576 

This system presents low quality light curves of the EB type. 
The depths of both minima are rather low and the RMS scatter 
is high. This is a detached system with components of unequal 
brightness and a circular orbit. The finding chart shows a close 
neighbour, at about r.'5 to the NNW, which might have dis- 
torted the temperature estimate of the binary components. 
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There are 12 out-of-eclipse spectra. The RVs curves are 
rather good and well constrained by the observations. 

A spectroscopic B luminosity ratio of 0.42 was found. This 
is close to the value reached by the final WD analysis. The 
observed composite spectra are very noisy (25 < S /N < 64) 
and are contaminated by nebular emission. The disentangled 
spectra are of rather poor quality. This is partly due to the low 
reliability of the continuum placement. The C n 4267 and Mg ii 
448 1 lines are visible in the spectrum of the primary. The tem- 
perature of the primary, determined from a least-squares fit 
where the temperature and luminosity ratios were fixed, de- 
pends on the quality of the fit, as shown by Fig. [T3] We adopted 
the temperature of the primary corresponding to the best;(f', in 
view of the roughly hnear correlation between Teff and x^, but 
without attempting to extrapolate the relation to;i'^ = 0. 
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Fig. 13. Effective temperatures of the system 5 180576 ob- 
tained from a fit to the composite spectra, keeping the tempera- 
ture and luminosity ratios fixed. Black dots are for the primary, 
red dots for the secondary. Note that the two curves should be 
considered as one and the same, since their ratio is constant. 



On the mass-log g diagram, the components are perfectly 
aligned along the 15 Myr isochrone. The position of the pri- 
mary component in the HR diagram appears slightly too high 
relative to the corresponding evolutionary track, while the sec- 
ondary falls right on its track. 

4.25. 5 185408 

This system presents medium quality light curves of the EA 
type. The minima are of similar depth and there is some ellip- 
soidal variation. This is a typical detached system with closely 
similar components, with medium quality light curves and low 



amplitude minima of similar depth. The orbit is circular. The 
/-band light curve was cleaned from five outliers before the 
PHOEBE/WD analysis. The finding chart shows a well isolated 
target, the closest neighbour being at about 3" to the WNW. 

There are 12 out-of-eclipse spectra. The RV curves are ex- 
cellent, especially in view of the faintness of the object, and 
well constrained by the observations. 

A spectroscopic B luminosity ratio of ~0.64 was found. The 
WD analysis tended to raise this value and the potential of the 
primary had to be fixed to keep it close to the spectroscopic 
one. The C ii 4267, Mg ii 448 1 and Si iii 4553 lines are visible in 
both disentangled spectra. The temperatures were obtained by 
fixing the temperature and luminosity ratios to the photometric 
values. 

On the mass-log g diagram, the components are fairly well 
aligned on the 10 Myr isochrone. On the HR diagram, both 
components are significantly overluminous with respect to their 
respective evolutionary tracks, as is often the case. 



4.26. 5 196565 

This detached eccentric system presents medium quality light 
curves with eclipses of similar depth. The V and B data are 
missing in the OGLE database. The finding chart shows a 
neighbour of similar brightness as the target about 2" to the 
S, raising some concern regarding possible spectral pollution. 
The results do not confirm these fears, however. 

There are 13 out-of-eclipse spectra. The RV curves are very 
good and well constrained by the observations, in spite of a 
lack of observations where the expected amplitude reaches its 
maximum. 

A spectroscopic B luminosity ratio of ~0.41 was found, and 
the potential of the primary was adjusted and fixed in the WD 
analysis to maintain that value. The Cn 4267 and Mgii 4481 
lines are visible in both disentangled spectra. The least-squares 
fit of composite synthetic spectra to the observed ones provided 
a temperature ratio larger than that provided by the photometry, 
and a primary temperature of 20960 ± 600 K. Fixing the tem- 
perature ratio to the photometric value resulted in the slightly 
hotter primary temperature that has been adopted. 

The mass-log g diagram suggests that the secondary has a 
too large radius relative to the primary, as if the luminosity ra- 
tio was overestimated: while the primary lies close to the 40 
Myr isochrone, the secondary lies on the 50 Myr one. On the 
HR diagram, however, each component has a position com- 
patible with its respective evolutionary track within the error 
bars, though the primary is slightly overluminous while the sec- 
ondary is slightly underluminous. 

This system was studied by WWOl . Their relative radius of 
the primary (0.204) is very close to ours (0.200) but their radius 
of the secondary (0.1 16) is smaller than ours (0.150). The dif- 
ference is due to the fact that they consider the eclipses as total, 
while we consider them as partial. Additional accurate photom- 
etry in both minima would be welcome to settle the question. 
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4.27. 5 261267 

This is a typical semi-detached system with high quahty Ught 
curves and ecHpses of unequal depths. The O - C curve shows 
no detectable depression before the primary eclipse. Although 
the finding chart shows a few neighbours, they all lie beyond 
2" of the target. 

This system has peculiar colour indices (see Table [3] and 
Fig.[3]l. Since it lies near the edge of the CCD in the OGLE-II 
survey, it is also listed under the name 6 11806 in the corre- 
sponding database. Thus, there are more than 600 data points 
in the /-band light curve, instead of about 300. 

There are 10 out-of-eclipse spectra, which constrain the RV 
curves relatively well. 

Both photometric and spectroscopic ratios are very similar 
(~0.4). Because of the relatively high Vrot sin / (>150 km s"') 
and moderate S /N of the observations, there are no exploitable 
metallic lines in the disentangled spectra. The temperature of 
the primary was obtained, as usual, by least-squares fit to the 
composite spectra, after fixing the temperature and luminosity 
ratios to the values given by a preliminary WD analysis. 

The mass-log g and HR diagrams are typical of a massive 
Algol-type binary. The primary is overluminous relative to the 
evolutionary track of an isolated star of the same mass, as are 
other semi-detached systems like 4 113853 and 5 277080. 

4.28. 5 265970 

This slightly eccentric detached system has medium-to-high 
quality light curves of the EA type. The sampling of the light- 
curve is incomplete, due to an orbital period close to 3.5 
days. In particular, the depth of the secondary minimum is ill- 
defined. Therefore, the photometric temperature ratio and incli- 
nation are not very reliable. Actually, there is a correlation be- 
tween these two parameters, in the sense that an increase of in- 
clination implies a decrease of the temperature ratio. The find- 
ing chart shows a well defined target, but with an only slightly 
fainter neighbour about 2'.'6 to the NW. 

This system lies near the edge of the OGLE-II CCD, and 
so was measured also on the adjacent chip under the name 6 
17345, so that there are as many as 586 data points in the / 
band. The two data sets were merged after applying a small 
magnitude offset to each. Using the EBOP code, the fitted mag- 
nitude at quadrature and its error was defined for each set, then 
the mean magnitude at quadrature weighted by the inverse of 
the variance was computed. Finally, the appropriate offset was 
applied to each of the two sets to adjust it to this mean magni- 
tude. 

There are 10 out-of-eclipse spectra. The RV curves are 
quite good and well constrained by the observations. 

Because of the loose constraints on the light curves, the 
proposed solution for this system heavily relies on the spectro- 
scopic observations. Nevertheless, the proposed solution meets 
very well the spectroscopic and photometric constraints, so we 
consider it as close to reality. 

A spectroscopic B luminosity ratio of ~0.23 was found 
from the usual least-squares fit, which provides the tempera- 
ture of the primary with an excellent internal precision. The 



temperature of the secondary is much less certain, because of 
the small luminosity ratio. In order to maintain the luminos- 
ity ratio to the spectroscopic value in the WD analysis, one 
has to fix the potential of the primary to an appropriate value, 
and the temperature ratio converges to a smaller value than the 
spectroscopic one, but still compatible with it given the errors. 
Thus the photometric temperature ratio was adopted, and the 
temperature of the primary determined in the usual way. 

A number of metallic lines are identifiable in both disen- 
tangled spectra: Cii 4267, On 4276, On 4415-17, Mgii 4481 
and Sim 4553. There is no emission in the Balmer lines. 

The positions of the stars in the mass-log g diagram are not 
quite mutually consistent: the ratio of radii should be decreased 
in order to bring the two components on the same isochrone, 
which would correspond to about 26 Myr On the HR dia- 
gram, the primary component matches its theoretical evolution- 
ary track surprisingly well, while the secondary has a position 
consistent with its track within the error bar The spectroscopic 
constraints, which are strong, are well fulfilled, but additional 
photometric data would be useful to improve our solution. This 
system is especially interesting, because the primary is very 
close to the TAMS while the secondary is much less evolved. 

We applied the EBOP code on four subset of the total time 
series, after fixing all parameters to their average value, except 
inclination, e cos oj, e sin w, magnitude at quadrature and phase 
shift. No significant trend can be seen on Fig. [141 which does 
not prove, however, that ca remains constant with time. The less 
reliable e sin o) quantity does not differ significantly from zero 
according to the EBOP code. In the WD solution, we have ar- 
bitrarily imposed a small value w = 1.8 x 10 "* radday ' which 
roughly corresponds to the theoretical prediction. 

4.29. 5 266015 

This system presents high quality light curves of the EB type 
with minima of unequal depth. Its bright /-band light curve is 
perfectly well sampled. The small depression occurring just be- 
fore the primary eclipse is indicative of a semi-detached system 
with a secondary component filling its Roche lobe and pouring 
matter onto the primary component. However, the amplitude of 
this effect was judged too small to justify an attempt to model 
it through a spot. 

There are 10 out-of-eclipse spectra. The RV curves are well 
constrained by the observations. 

A spectroscopic B luminosity ratio of 0.61 was found, 
rather remote from the photometric value (0.49), as is often 
the case in semi-detached systems. The disentangled spectra of 
both components have a decent S /N (Fig. [TsTl. Nevertheless, 
the Vrot sin / being high (>150 km s"'), the small metallic lines 
are not very conspicuous. The Hen 4120 and 4542 lines are 
visible in the spectrum of the primary. The Si iii line appears 
in both spectra. This is one of the few systems with apparently 
no significant nebular emission lines, therefore the disentan- 
gled Balmer lines can be used to find the temperature of the 
primary. The spectral features of the primary point to a 32 000 
K star, which is confirmed by the usual least-squares fit. 
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Fig. 14. Same as Fig.|7j but for the system 5 265970. The hor- 
izontal dashed line shows the weighted average of the e cos w 
values. The solid line corresponds to w ^ 3.8° yr', the ex- 
pected theoretical value. 



The positions of the stars in the mass-log g and HR dia- 
grams are coherent with an evolved system having undergone 
mass exchange. The primary component is close to the track 
corresponding to a 15.6 Al© star. 

This system was studied by WW02. As with most semi- 
detached systems they studied, their photometric mass ratio 
proved to be unreliable. 

4.30. 5 266131 

The most striking feature of this slightly eccentric detached 
system (e ^ 0.04) is its huge apsidal motion, which explains 
the apparently very bad /-band O — C curve. An examination 
of the OGLE finding chart revealed that this binary is slightly 
blended, but a posteriori, this does not seem to have distorted 
the results. 

This system lying close to the edge of the CCD in the 
OGLE-II survey, it also exists under the name 6 22883, which 
doubles the number of data points in the / band. The magnitude 
at quadrature was determined using the EBOP code for each 
of the two data sets; the two resulting values agreeing within 
one thousandth of a magnitude, the two data sets were merged 
without applying any magnitude off'set. 

There are 10 out-of-eclipse spectra, which constrain quite 
well the RV curves. 

A spectroscopic B luminosity ratio of ~0.65 was found. As 
usual with the well detached systems, the potential of the pri- 
mary was fixed at a value such that this luminosity ratio is pre- 
served through the WD analysis. The S /N of the disentangled 



spectra are too low to allow any useful metallic line to be seen. 
Since the temperature ratio is well constrained by the photo- 
metric data, the temperature of the primary was determined by 
fixing this ratio at its photometric value. There is no nebular 
emission in the Balmer lines. 

The mass-log g diagram suggests that the ratio of radii has 
been slightly underestimated, since the secondary has a too 
large surface gravity compared to the primary. Still, the po- 
sitions of both components are compatible, within the errors, 
with an isochrone at 7 - 8 Myr. The HR diagram shows a good 
match between the positions of both components and the evo- 
lutionary tracks of single 9.0 and 7.7 Al© stars. However, both 
components are slightly overluminous relative to their respec- 
tive tracks. 

The apsidal motion amounts to w ^ 50 ° yr"\ as evidenced 
by the e cos a> values obtained with the EBOP code for four 
successive subsets of the whole time serie in the / band. That 
value is confirmed by the WD analysis. The run of e cos w ver- 
sus time is shown in Fig. [16] together with the best fit curve 
provided by the WD code. The reality of a very fast apsidal 
motion is beyond any doubt. It is further discussed in Section 
5. 
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Fig. 16. Same as Fig.|7l but for the system 5 266131. The solid 
line is the best fit provided by the WD code. 



4.31. 5 266513 

This a close detached system with similar components, accord- 
ing to its medium quality light curves with minima of simi- 
lar depth. The finding charts shows a slightly fainter star at 
1 .2 - 1 .3" to the W of the target. 
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Fig. 15. Sections of the disentangled spectra of the primary and secondary components of the binary 5 266015. The scale on the 
vertical axis is adapted to the size of the features in each section. The red line represents a velocity-shifted synthetic spectrum 
of the star. Beside prominent H i and He i lines. He ii, O n, Mg ii, Si iii and Si iv lines are visible in both spectra. The strong He ii 
4542, next to a smaller Si in 4553 line, and the He ii 4200 line confirm the high temperature of the primary (~32 000 K). Poor 
continuum placement is responsible for the apparent O - C mismatch of the Hd line. 



There are 10 out-of-eclipse spectra. The RV curves are well 
constrained by the observations, though the noise is rather large 
because this system is the faintest in our sample. 

A spectroscopic B luminosity ratio of 0.73 was found, but 
with a large scatter of 0.09. The potential of the primary was 
fixed so that the WD analysis preserves a ratio close to that 
value. The disentangled spectra have a low S /N and no metal- 
lic fines are exploitable. Moreover, the Balmer lines are pol- 
luted by nebular emission. The photometric temperature ratio 
is rather well defined, thanks to the large depth of the minima, 
and was fixed for the determination of the temperature of the 
primary, as usual for most detached systems. 

The mass-log g diagram shows the two components fairly 
well aligned along the ~ 15 Myr isochrone. On the HR dia- 
gram, the primary Ues almost exactly on its evolutionary track, 
while the secondary is slightly overluminous, though its posi- 
tion is quite compatible with the evolutionary tracks within the 
errors. The luminosity ratio (hence the ratio of radii) seems to 
have been slightly overestimated. 

4.32. 5 277080 

This system is a typical semi-detached binary with high quality 
light curves of the EB type. A small depression before the pri- 
mary eclipse signals the possible presence of a mass-transfer 
stream. The effect is small enough, however, that we did not 
deem it worth the eff'ort to model it through a cool spot on the 
primary. 

This system has peculiar colour indices (see Table [3] and 
Fig. |3]l. Indeed, an examination of the OGLE finding chart re- 



vealed that this binary is strongly blended with a brighter star 
lying r.'2 to the S W of the target. 

There are 1 1 out-of-eclipse spectra. Because of an orbital 
period close to 2 days, there are no observations before the 
quadratures. In spite of that, the RV curves are rather well con- 
strained. This is one of the few systems showing no nebular 
emission in the Balmer lines. 

A spectroscopic B luminosity ratio of 0.7 was found. This 
is higher than the photometric ratio (0.57), as seems to be of- 
ten the case of semi-detached systems. The disentangled spec- 
tra are fairly good. A number of metallic lines are visible in the 
spectrum of the primary: O ii 4185, C ii 4267, ii 4276-7, Mg ii 
4481 and Sim 4553. The Mgii 4481 line is the only metalfic 
line detectable in the spectrum of the secondary. The tempera- 
ture of the primary was determined after fixing the temperature 
and luminosity ratios to their photometric values. 

The positions of the stars in the mass-log^ and HR dia- 
grams are typical of Algol-type system with a secondary more 
evolved than the primary. The primary component, however, is 
much more luminous than expected from the evolutionary track 
of a single star of the same mass, and we cannot guarantee that 
its effective temperature has not been overestimated. The pe- 
culiar HR diagram might be the result of the blend mentioned 
above with a nearby bright star. 

4.33. 5 283079 

This is a typical detached system with twin components and 
a circular orbit. The target appears just isolated on the finding 
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chart, with a companion of similar brightness about 275 away 
from it to the SE. 

There are 10 out-of-eclipse spectra. The RV curves are well 
defined and constrained by the observations. There is probably 
a faint nebular emission in both Hy and H6 lines. 

With a spectroscopic B luminosity ratio close to 1 .0, a mass 
ratio of 1.003 and a temperature ratio of 0.997, this is the sys- 
tem with the most similar components in our sample. Thus it 
is probably a real pair of "twins", i.e. it has a mass ratio larger 
than 0.95. This is why the choice of the primary is undecided 
in this system: which component is the primary was chosen on 
the basis of an early iteration, and it is only in the last iteration 
that we obtained q = 1 .003 > 1 . The potential of the primary 
was fixed in the WD analysis, in order to keep the B luminosity 
ratio close to 0.99. The disentangled spectra are noisy and the 
C 11 4267 and Mg ii 448 1 are barely visible. The temperature of 
the primary was determined by fixing the temperature ratio to 
0.997 and the luminosity ratio to 0.99. 

On the mass-log g diagram, both components lie on the 3 
Myr isochrone, so this system is very young. On the HR dia- 
gram, they are slightly overluminous with respect to the evo- 
lutionary tracks, but still within the error bars. On the other 
hand, the components lie exactly on the metal-poor tracks 
(Z = 0.001). 

5. Discussion 

5.1. Common systems with HHH03/05 

Five systems are common to the HHH03/05's sample and 
ours: two detached systems (4 163552 and 5 038089) and 
three semi-detached ones (4 1 10409, 5 026631 and 5 277080). 
Comparison of a number of parameters is shown in Table |5] 
The good news is that we agree on the model, i.e. the configu- 
ration of the system in relation with the Roche lobes. The bad 
news is that, in several cases, we do not agree on most param- 
eters, except the inclination, a purely photometric parameter. 
There is nothing surprising about that, as we have the same 
light curves. 

The first system, 4 1 10409, shows a fair agreement for the 
masses (excellent for the primary, within 1.6 cr for the sec- 
ondary), the radii (within 2cr) and the surface gravities (1 cr 
or better). By contrast, the effective temperatures differ by as 
much as 7000 (primary) and 6000 K (secondary), those given 
by HHH05 being clearly underestimated. 

The second system, 5 163552, shows a good agreement 
of the temperatures, the estimates by HHH05 being only ~ 
1000 K higher than ours. On the other hand, the masses given 
by HHH05 are 3 cr larger than ours, and their radii are larger 
too. Part of the problem might arise from HHH05 neglect- 
ing the third light. In addition, HHH05 may have been more 
sensitive to crowding effects on their RV curves, since the di- 
ameter of the 2dF fibers they used is 2", while the diame- 
ter of the FLAMES/GIRAFFE fibers is only 1'.'2. But nebular 
emission may be the main reason for the discrepancy: HHH05 
overestimate the amplitude of both RV curves (relative to us), 
and indeed nebular emission can mimick a wider separation of 
B aimer lines near quadratures. 



For the semi-detached system 5 26631, there is a perfect 
agreement about the mass of the primary, while the mass of the 
secondary given by HHH05 is 4 cr larger than ours. Their solu- 
tion appears less realistic than ours, because they obtain a mass 
ratio close to unity, which would be exceptional for a semi- 
detached system. Similarly, we agree on the radius of the pri- 
mary, while HHH05 overestimate the radius of the secondary. 
HHH05 estimate temperatures lower by ~ 3000 K (primary) 
and by ~ 4000 K (secondary) than us. The amplitude of their 
RV curve is overestimated for the primary, but slightly lower 
than ours for the secondary. Again, nebular emission might be 
responsible for the discrepancy. 

The fourth system, 5 38089, shows a very large discrepancy 
in the mass estimates: for the secondary, it amounts to 63% if 
we take our value as the reference. Both masses are overesti- 
mated by HHH03, while the radius of the primary is in perfect 
agreement, and that of the secondary is significantly overesti- 
mated. On the other hand, the temperatures are in nearly per- 
fect agreement (within 2%) for the primary, and slightly dis- 
agree (by 1840 K) for the secondary. It appears that HHH03 
have strongly overestimated the amplitude of the RV curve of 
the primary, while that of the secondary is much closer to ours. 
Here, one cannot invoke any systematic effect linked with neb- 
ular emission, since the latter cannot be seen on our spectra. 
On the other hand, this system has peculiar colour indices (Fig. 
|3]l and its position is unsatisfactory on our own HR diagram. 
Crowding problems causing what could be called "spectro- 
scopic third light" are probably the main source of problems 
in both studies, and may have been enhanced in that of HHH03 
by the wider 2dF fiber. 

The last system, 5 277080, holds the record of mass dis- 
crepancy: relative to our values, HHH05 overestimate the mass 
by factors as large as 1.8 and 2.3 for the primary and the sec- 
ondary respectively. The radii are overestimated as well, but 
by a much smaller amount. The temperatures differ by more 
than 4000 K. Nebular emission cannot be the culprit, since it 
remains undetectable in our spectra. 

The problem lies in the RV curves and the disentangled 
spectra. There are three basic spectroscopic parameters on 
which all the subsequent analysis is built: the two velocity 
semi-amplitudes and one of the temperatures (the other one 
being constrained by the photometric mean surface brightness 
ratio). If the velocity semi-amplitudes are wrong, the masses, 
radii, ^-ratio, logg-values and luminosities are wrong. The lu- 
minosities are obviously highly dependent on the temperatures, 
despite a partial cancelling of a possible error by the bolomet- 
ric correction (see Section l3.13l l. We are confident on our val- 
ues for the velocity semi-amplitudes and the primary tempera- 
ture. The RV curves of the five systems shared with HHH03/05 
are shown in Figs [33l[34l and[35] They are well sampled, there 
are good constraints at the quadratures and the O - C curves 
are flat (except for 5 277080 where some systematics appears). 
Consequently, there is no reason to doubt about the quality of 
our velocity semi-amplitudes. The RV curves are not shown 
in HHH03/05, nevertheless there are certainly two reasons to 
suspect that they have a lower quality than ours. The first one 
is instrument-related. As already mentioned, our instrumental 
setup has better resolution (0.67 A instead of 2 A). Moreover, 
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Table 5. Comparison of our solutions (NGBR) for the five binaries in common with HHH03/05. 



Object 



(Ato) 



Als 

(Mo) 



Re 

(Rq) 



(Rq) 



(K) 



(K) 



loggp 
(dex cgs) 



loggs 
(dex cgs) 



4 110409 


NGBR 


sd 


13.93 ±0.59 


7.30 ± 0.36 


4.64 ±0.14 


7.81 ±0.18 


32370 ±816 


21 170 ±583 


4.25 ±0.03 


3.52 ±0.03 


77.3 ±0.1 




HHH05 


sd 


13.7+0.8 


8.9 ±1.1 


4.3 ±0.2 


8.4 ± 0.3 


25500 (fixed) 


15390 ± 160 


4.30 ±0.05 


3.54 ±0.06 


76.8 ± 0.3 


4 163552 


NGBR 


d 


9.56 ±0.61 


9.10 ±0.57 


5.06 ±0.16 


4.54 ±0.15 


24400 ± 990 


24330 ± 993 


4.01 ±0.04 


4.08 ± 0.04 


80.9 ± 0.2 




HHH05 


d 


13.3 ±1.0 


124± 1.1 


5.3 ±0.1 


5.2 ±0.2 


25500 (fixed) 


25390 ± 190 


4.11 ±0.04 


4.10 ±0.05 


78.3 ±0.1 


5 026631 


NGBR 


sd 


11.40 ±0.34 


7.86 ±0.25 


5.19 ±0.08 


4.93 ± 0.08 


28670 ± 394 


21050 ±356 


4.07 ± 0.02 


3.95 ±0.02 


61.4 ±0.1 




HHH05 


sd 


11.5 ±0.6 


11.3 ±0.8 


5.1 ±0.1 


5.7 ±0.1 


25500 (fixed) 


17130 ±300 


4.08 ± 0.03 


3.99 ±0.04 


61.5 ±0.1 


5 038089 


NGBR 


d 


13.01 ±0.20 


ll.70±0.18 


5.98 ±0.08 


4.94 ± 0.07 


30660 ±130 


31020 ±301 


4.00 ±0.01 


4.12 ±0.01 


77.1 ±0.2 




HHH03 


d 


17.1 ±1.5 


19.1 ± 1.6 


6.1 ±0.2 


6.1 ±0.3 


30100 (fixed) 


29180 (fixed) 


4.10 ±0.05 


4.15 ±0.05 


76.9 ± 0.3 


5 277080 


NGBR 


sd 


9.68 ±0.55 


4.88 ± 0.29 


4.72 ±0.09 


5.13 ±0.09 


29780 ± 480 


20410 ±372 


4.08 ± 0.03 


3.71 ±0.03 


76.6 ±0.1 




HHH05 


sd 


17.4 ±0.9 


1 1.3 ± 1.0 


5.0 ±0.1 


6.8 ±0.2 


25500 (fixed) 


15890 ±150 


4.27 ±0.03 


3.82 ±0.05 


76.8 ±0.3 



Table 6. Comparison of our results (NGBR) with the light-curves solutions of Wyithe & Wilson (I2001II2002I I (WWOl, WW02) 
and Graczyk (2003) (dG03). The following parameters are listed: the model of the system, the mass ratio, the relative radii, the 
inchnation, the eccentricity, the / luminosity ratio, the eclipse parameter (see Eq. 2 in WWOl) and the third light (normalized to 
total flux at phase 0.25). 



Object 


Ref. 


Model 


1 


Rp/a 


Rs/a 


iC) 


e 


(Ls/Lp), 


Fe 


hj 


4 117831 


NGBR 


d 


0.981 ±0.013 


0.261 ± 0.008 


0.287 ± 0.008 


78.2 ±04 





1.175 


0.599 







WW02 


sd 


0.157 ±0.054 


0.357 ± 0.032 


0.209 ±0.021 


77.9 ±2.1 


0.049 ±0.019 


0.372 


0.853 





4 163552 


NGBR 


d 


0.952 ±0.015 


0.339 ±0.011 


0.305 ±0.010 


80.9 ± 0.2 





0.798 


0.797 


0.1 19 ±0.006 




dG03 


d 


0.98 ± 0.09 


0.338 ± 0.006 


0.326 ± 0.008 


854 ±0.9 





0.951 


0.895 


0.263 ± 0.007 




WW02 


sd 


0.207 ±0.010 


0.383 ± 0.004 


0.234 ± 0.003 


82.8 ±0.7 


0.004 ± 0.002 


0.397 


1.051 





4 175333 


NGBR 


d 


0.781 ±0.014 


0.237 ± 0.008 


0.203 ± 0.007 


804 ±0.5 


0.022 ± 0.007 


0.582 


0.675 







WWOl 


d 


1.0 


0.253 ± 0.009 


0.143 ±0.007 


85.0 ± 1.5 


0.012 ± 0.003 




1.080 





5 016658 


NGBR 


d 


0.889 ±0.011 


0.320 ±0.011 


0.255 ±0.010 


80.3 ± 0.2 





0.609 


0.798 







WWOl 


d 


1.0 


0.345 ± 0.007 


0.200 ± 0.004 


85.7 ± 1.3 


0.002 ± 0.002 




1.175 





5 038089 


NGBR 


d 


0.899 ±0.011 


0.273 ± 0.004 


0.226 ± 0.003 


77.1 ±0.2 





0.693 


0.609 







dG03 


d 


1.12" 


0.246 ± 0.009 


0.255 ± 0.007 


76.8 ± 1.9 





1.053 


0.554 





5 196565 


NGBR 


d 


0.803 ± 0.007 


0.200 ± 0.005 


0.150 ±0.005 


83.5 ±0.2 


0.138 ±0.006 


0.43 


0.789 







WWOl 


d 


1.0 


0.204 ± 0.006 


0.1 16 ±0.003 


88.7 ±2.5 


0.085 ± 0.002 




1.282 





5 266015 


NGBR 


sd 


0.441 ±0.005 


0.332 ± 0.005 


0.310 ±0.005 


78.5 ±0.1 





0.494 


0.714 







WW02 


sd 


1.014±0.051 


0.183 ±0.003 


0.357 ± 0.004 


80.8 ± 04 


0.000 ± 0.002 


2.564 


1.039 





5 266131 


NGBR 


d 


0.864 ±0.013 


0.281 ±0.010 


0.242 ± 0.009 


83.5 ±0.2 


0.042 ± 0.001 


0.647 


0.847 







WWOl 


d 


1.0 


0.332 ± 0.008 


0.216 ±0.009 


85.1 ± 1.5 


0.003 ± 0.003 




1.071 





5 283079 


NGBR 


d 


1.003 ±0.013 


0.243 ± 0.007 


0.242 ± 0.007 


87.7 ± 0.3 





0.987 


0.918 







WWOl 


d 


1.0 


0.272 ± 0.006 


0.233 ± 0.007 


88.5 ± 1.8 


0.001 ±0.001 




1.028 






FromHHH03. 



for the systems in common, our data have a higher S /N (typi- 
cally a factor of about two, from the S /A^-calculators of the re- 
spective instruments), in spite of a lower S /N on average for the 
whole sample, as we observed mostly fainter systems. Thus our 
raw data are intrinsically better The second reason is related 
to the measurement of the radial velocities. Our radial veloci- 
ties were determined exclusively from the He i lines mentioned 
in Section [331 From the procedure described in HHH03/05, it 
seems clear that in most cases (i.e. with no conspicuous emis- 
sion lines) both Balmer and He i and/or the whole spectrum 
were used to determine the radial velocities. One can suspect 
that unrecognized nebular emission in the H i lines is partly re- 
sponsible for bad RV measurements. Actually, in a number of 
cases, nebular emission in the Balmer lines can easily be con- 
fused with the typical "SB2 cleavage" of blended lines. The 
disentangling of a composite spectrum being tightly linked to 
the radial velocities, the quality of the component spectra and 
their subsequent use for fixing the temperature of one com- 
ponent will be badly influenced by poorly determined radial 
velocities. 

An important point of divergence in the methodology used 
to "solve" a binary is the use (or non-use) of the spectral B lu- 
minosity ratio. Comparing the spectral and photometric values 
of the B luminosity ratio actually reveals that they often dis- 
agree for detached systems. The reason is the near-degeneracy 
of the ratio of radii with the inchnation in the case of partial 
ecUpses (Wyithe & Wilson 1200 11 1. In principle, the photomet- 



ric luminosity ratio of semi-detached systems should be better 
constrained, since they have less free photometric parameters 
than the detached ones (see Fig. 3 in Wyithe & Wilson 120011 
and l2002T l. The regression shown in Fig.[T7]clearly shows that 
the spectroscopic and photometric ratios of the population of 
semi-detached systems agree fairly well, except that the slope 
is slightly lower than one. The latter fact may be caused by 
the secondary tidal distorsion being larger than that of the pri- 
mary: the total luminosity, which is the parameter provided by 
the WD code, may be smaller than suggested by the apparent 
brightness at the quadrature. Among the detached systems. Fig. 
[T7]shows only those for which the spectroscopic and photomet- 
ric luminosity ratios were found to agree, so the good match is 
artificial. 



Thus, the spectroscopic luminosity ratios were adopted sys- 
tematically for the detached systems, except in the cases where 
they agree with the photometric one within the uncertainties 
(full dots in Fig. [TTl i. For the semi-detached and contact ones, 
the photometric ratio was preferred. In HHH05, the spectro- 
scopic ratio was only used to settle the case between two con- 
flicting photometric solutions. 



Finally, it must be emphasized that the quality of our spec- 
tral analysis owes greatly to the excellent library of synthetic 
NLTE O- and B-spectra. 
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Fig. 17. Photometric vs. Spectroscopic values of the B lu- 
minosity ratio (full symbol: detached; open symbol: semi- 
detached/contact). The dotted line corresponds to X.^g'^'^ = 
Xg °'. The best linear fits of the two samples are shown (solid: 
detached; dashed: semi-detached/contact ). Note that the de- 
tached systems shown are only those for which the photometric 
ratio agrees well with the spectroscopic one. 



5.2. Nebular emission and l<inematics 

Hy and H6 emission lines are visible in most binary spectra 
of our sample. FWHM values and heliocentric radial velocities 
are presented in Table |7] These are mean values inferred from 
spectra close to conjunction in order to separate the real emis- 
sion from the "SB2 cleavage". The strengths are not indicated, 
as in normalized spectra the nebular lines are diluted by the 
stellar continua. Therefore, it would not be possible to com- 
pare line strengths measured on two spectra of two different 
systems. Even in the case of two spectra of the same binary at 
two different epochs, the strength measurement is not reliable 
because of the dependency of the diluting stellar continuum on 
the seeing^. 

It is beyond the scope of this paper to make a detailed kine- 
matic study of the SMC. Nevertheless, it is interesting to exam- 
ine the distribution of the 20 RVs derived from the nebular Hy 
lines listed in Table |7] The distribution is illustrated in Fig.fTSJ 
The 5 km s"' -binned histogram shows a bimodal distribution 
with two modes at 145-150 and 160-165 km s"'. Since the size 
of the sample is modest, one can wonder about the significance 
of this observation. Therefore, the RVs of the nebular lines in 
the sky spectra were also investigated (21 spectra per epoch). 
The Hy data are presented in Table |8] It is readily apparent that 



^ Note that the sky subtracted from the stellar spectra was limited 
to the continuum component. Therefore, sky subtraction did not alter 
the nebular emission lines. 
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Fig. 18. Histogram of the mean heliocentric radial velocities of 
the nebular Hy lines observed in the spectra of 20 binaries. The 
sample is arranged in 5 km s ' bins. 



the RV-distribution shows two peaks as well, at 140-145 and 
160-165 km s ' (FigfT9]l. Therefore, these observations provide 
compelling evidence that the nebular H ii of the SMC displays a 
bimodal velocity distribution (at least in our small 25'-diameter 
field). 
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Fig. 19. Histogram of the mean heliocentric radial velocities of 
the nebular Hy lines observed in the 2 1 sky spectra. The sample 
is arranged in 5 km s"' bins. 

Comparing the systemic velocities of our 33 systems with 
the combined nebular Hii velocities (41 values) shows that the 
two samples have very different distributions. The nebular ve- 
locities were fitted by two Gaussian functions of equal ampli- 
tudes with modes at 145.2 and 161.4 km s ' and standard de- 
viations of 3.44 and 4. 14 km s"' , respectively. The correspond- 
ing FWHMs are about one fifth of a resolution element (~47 
km s"'). The systemic velocities are best fitted with a Gaussian 
function centred on 162 km s ' with a standard deviation of 
12.9 km s ' . That last fit is less convincing than for the H i, as 
there is less data and the velocity dispersion is 3-4 x higher. 
Further nebular and systemic velocity data from the analysis of 
the remaining SBl systems contained in the GIRAFFE field is 
expected to substantially improve the statistics. 
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Table 7. Nebular emission: mean radial velocities and mean FWHM of the Hy and H6 lines. The quoted uncertainties are 
standard deviations. 



Object 


Hy 


ffi 




Object 


Uy 


H(5 






Kad 


FWHM 


V,-ad 


FWHM 


V,-ad 


FWHM 


Vrad 


FWHM 




(km s-') 


(A) 


(kms"') 


(A) 




(kms"') 


(A) 


(kms"') 


(A) 


4 110409 


- 


- 


- 


- 


5 100485 


159.1 ±4.7 


1.01 ±0.10 


166.9 ±17.1 


0.83 ±0.19 


4 113853 


142.6 ±1.3 


0.93 ± 0.05 


147.4 ± 6.7 


0.77 ± 0.08 


5 100731 


163.7 ± 3.6 


1.02 ±0.12 


161.6 ±7.3 


0.95 ±0.11 


4 117831 


145.8 ±1.8 


0.99 ± 0.08 


152.2 ± 14.2 


0.92 ±0.13 


5 106039 


159.1 ±2.5 


0.91 ±0.03 


159.3 ±6.3 


0.83 ± 0.09 


4 121084 


147.9 ± 2.0 


1.22 ±0.04 


152.9 ±2.7 


1.05 ±0.02 


5 111649 


149.4 ± 3.0 


0.77 ± 0.09 


152.0 ± 12.9 


0.64 ± 0.08 


4 121110 


136.8 ±1.4 


1.44 ±0.07 


141.4 ±6.0 


1.20 ±0.07 


5 123390 


- 


- 


- 


- 


4 121461 


147.6 ± 2.0 


1.11 ±0.07 


150.5 ±0.6 


1.00 ±0.06 


5 180185 


164.5 ± 5.7 


0.75 ± 0.08 


- 


- 


4 159928 


146.7 ±1.1 


0.99 ± 0.03 


148.9 ±2.5 


0.94 ± 0.06 


5 180576 


157.4 ±2.9 


0.99 ±0.15 


163.8 ±3.7 


0.96 ±0.13 


4 160094 


143.6 ±1.2 


1.03 ±0.03 


144.2 ± 3.7 


0.95 ± 0.03 


5 185408 


- 


- 


- 


- 


4 163552 


143.6 ±0.3 


0.96 ± 0.02 


147.1 ±2.4 


0.95 ± 0.04 


5 196565 


- 


- 


- 


- 


4 175149 


- 


- 


- 


- 


5 261267 


160.7 ± 3.9 


0.88 ±0.14 


170.3 ±7.6 


0.78 ± 0.22 


4 175333 


146.9 ±3.3 


1.11 ±0.07 


151.1 ±1.5 


0.99 ± 0.02 


5 265970 


- 


- 


- 


- 


5 016658 


168.3 ± 9.4 


1.11 ±0.06 




- 


5 266015 


- 


- 


- 


- 


5 026631 


148.8 ±5.7 


0.93 ± 0.21 


151.8 ±4.0 


0.66 ± 0.04 


5 266131 


- 


- 


- 


- 


5 032412 


- 


- 


- 


- 


5 266513 


- 


- 


- 


- 


5 038089 


- 


- 


- 


- 


5 277080 


- 


- 


- 


- 


5 095337 


171.6±3.1 


1.07 ±0.11 


172.7 ±6.5 


0.87 ±0.12 


5 283079 


- 


- 


- 


- 


5 095557 


164.2 ±2.5 


1.01 ±0.07 


- 


- 
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Fig. 20. Velocity field of the nebular H ii. The size of symbol 
is proportional to the value of the velocity. Values greater than 
155 km s ' are in red (electronic version only). Filled sym- 
bols are for binary systems while open symbols are for skies. 
Crosses are for binary systems showing no emission. 



It is worth mentioning that the observed bimodal distribu- 
tion of H II is compatible with the kinematic study of H i in the 
SMC by Stanimirovic et al. ( 12004] ). even though the gas we see 
in our work is certainly hotter than the gas observed at 21 cm 
wavelength. Indeed, the authors mention the existence of a bi- 
modal velocity field with central velocities of 137 and 174 km 
s ' . Interestingly, Fitzpatrick ( I1985I I already mentioned two H i 
complexes with velocities 1 34 and 1 67 km s" ' , which our emit- 
ting gas is probably related with. Figure|20]shows the mean (he- 
liocentric) velocity field obtained from nebular H ii data. Filled 
and open symbols correspond to binary and sky (nebular) data, 
respectively. The size of the symbol is drawn according to the 
value of the radial velocity. A cross indicates the position of 



a binary without significant emission. The map has the same 
orientation as in FiglT] We can tentatively draw an oblique line 
from the NE to the S W corner demarcating high-velocity points 
from low-velocity points. This line has a ~ 60° slope in the EN 
(clockwise) direction. On the east (i.e. left) side of this line, 
nebular RVs are greater than ~ 155 kms"'. This velocity gra- 
dient is consistent with Fig. 3 in Stanimirovic et al., i.e. the 
velocity of H i increases from west to east in this part of the 
SMC. Therefore, the hot and cool phases of the gas seem to 
share very similar spatial and kinematic distributions. 
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Fig. 21. Comparison between two kinematic tracers: radial ve- 
locities of nebular H ii (left) and systemic velocities of binary 
systems (right). The data are arranged in 5 km s"' bins. The 
best-fitting double-Gaussian and Gaussian models are over- 
plotted. 
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Table 8. Sky positions: astrometry, sky background excesses 
(see Section 15.2b and mean radial velocities derived from the 
Hy emission line. 



id 


a (J2000) 


6 (J2000) 


Asky 


V„d 




D 


(°) 


(%) 


(kms-') 


SI 


12.20608 


-73.17275 


+3.8 ±5.4 


143.9 ±4.2 


S2 


12.68487 


-73.18533 


+ 13.4±4.1 


169.4 ±5.8 


S3 


12.43554 


-73.19611 


-5.8 ±4.7 


152.8 ±3.7 


S4 


12.49946 


-73.14372 


-4.8 ± 3.8 


160.0 ± 4.4 


S5 


12.49446 


-73.17817 


+0.5 ± 7.0 


155.8 ±6.2 


S6 


12.34712 


-73.24400 


-11.4±2.0 


156.0 ±3.7 


S7 


12.40693 


-73.25211 


-3.5 ±1.5 


145.2 ±5.5 


S8 


12.63725 


-73.31728 


+ 19.8 ±5.4 


161.9 ±2.4 


S9 


12.61700 


-73.30436 


+7.1 ±5.5 


162.2 ± 6.5 


SIO 


12.60371 


-73.21592 


-6.8 ±3.2 


161.4 ±7.3 


Sll 


12.83250 


-73.26103 


+ 15.9 ±3.8 


160.9 ±6.1 


S12 


12.48021 


-73.10711 


+4.2 + 3.1 


141.1 ±2.7 


S13 


12.85400 


-73.15156 


+34.9 ± 7.3 


168.2 ±6.6 


S14 


12.80883 


-73.29111 


-1.5 ±2.7 


161.0 ±3.9 


S15 


12.17254 


-73.13839 


-6.0 ±3.4 


142.8 ± 2.9 


S16 


12.09937 


-73.12347 


-5.7 ±3.0 


143.2 ± 1.1 


S17 


12.16858 


-73.11614 


+7.3 ±5.5 


134.9 ±2.2 


S18 


12.14775 


-73.09939 


+ 10.2 + 6.1 


127.1 ± 1.5 


S19 


11.95912 


-73.25294 


+20.2 ± 5.2 


142.1 ±4.0 


S20 


12.01675 


-73.22519 


+25.4+11.8 


146.3 ± 3.2 


S21 


12.13529 


-73.35828 


-2.3 ±6.8 


154.3 ± 3.6 



5.3. Comparison with evolutionary models and the 
mass-luminosity relation 

Comments on the individual systems were given in Section |4] 
Here, the collective properties of the sample are examined more 
thoroughly. For each of the 33 systems, Table |9] lists the dif- 
ferences AA1 between the evolutionary mass interpolated, in 
the HR diagram, in the evolutionary tracks from Charbonnel et 
al. ( 11993) . and the observed dynamical mass determined from 
the simultaneous RV-curves and light-curves analysis. Figure 
l22l shows these differences versus the observed masses for the 
23 binaries that have a detached configuration. The agreement 
between the observed and evolutionary masses is acceptable, 
though some systematics does appear. The most discrepant ob- 
jects are the secondary components of 5 38089 and 5 123390, 
and the primary of 5 38089. These systems have a peculiar po- 
sition in the HR diagram, perhaps because of third-light con- 
tamination. 

For comparison, the sample of 18 detached systems of 
HHH05 is shown on Fig.|23] Figures l22l and |23] are drawn at 
the same scale. Beside the very large scatter of HHHOS's AA1 
values (they span a range of almost 10 Mq), the most striking 
feature is the correlation between AA1 and the observed mass 
for a number of primary and secondary components, especially 
between 7 and 15 Al©. Explaining this trend is not straightfor- 
ward, but it might be related to the temperature determination. 
Indeed, HHH05 mention in their Section 3.1.2: "We found that 
the O — Bl stars were quite easy to classify, but that spectral 
classifications were more difficult for later B types due to the 
absence of He ii lines and, in these low metallicity stars, the 
lack of detectable Si ii and other metal lines". We faced the 



same difficulty, but with much less severe consequences. The 
larger telescope and better resolution clearly have brought a 
dramatic improvement in the quality of the results. 

There is a slight trend in our mass residuals: first of all, most 
AM values are positive. Second, while the mass difference is 
negligible for low masses (around 4-5 Al©), it seems to gradu- 
ally increase to almost 2 Mc at about 13 Mo- The single point 
at 17 AIq corresponds to AM ~ 1 Mq only, but it is not suffi- 
cient to define the possible relation beyond 13 Mo- It might be 
part of a parallel sequence defined by 4 other stars. One may 
speculate that, if real, the latter sequence might correspond to 
e.g. systems for which the sky background was overestimated, 
which would have caused an underestimate of Teff. The main 
trend is surprising, because the more massive systems are also 
more luminous on average and so appear brighter since they 
all lie at about the same distance. Hence the signal-to-noise ra- 
tio of their spectra is larger, and the determination of the T^ff 
should be more reliable. On the other hand, Teff is better con- 
strained in stars cooler than ~ 20 000 K, because the He i lines 
increase in strength with increasing temperature, while the H i 
Balmer lines decrease. In hotter stars, both the H i and He i lines 
decrease with increasing temperature, so that even a slight error 
on the continuum placement or (equivalently) on the sky sub- 
traction may cause a significant error on Teff . At the hot end of 
the temperature range of our sample (i.e. for Teff > 30 000 K), 
the effective temperature is again better constrained, thanks to 
the fast raise of He ii lines. Therefore, one may expect that the 
least reliable temperatures - hence the least reliable luminosi- 
ties - will occur for stars with intermediate Teff or masses. This 
might explain why the scatter of the AM values in Fig. l22]is 
large between 7 and 14 Mq. 

The fact that AM > in most cases is not easy to explain. 
Possible explanations can be grouped in three categories: 

- Systematic error in the T^ff determination: As mentioned 
above, the temperature determination is very sensitive to 
continuum placement or sky subtraction, especially when 
all lines vary in the same way with Teff. The continuum 
has been determined both automatically and manually for 
each star, and we have verified that the least-squares fit used 
to determine the temperature always gives the same result, 
unless the continuum is clearly wrong. All fits were done 
under visual monitoring, so the quality of the continuum 
placement could be judged with confidence. In addition, 
plotting the estimated temperature against the sum of the 
squared residuals generally results in a horizontal line, as 
shown in Section |4] (the system 5 180576 is a notable ex- 
ception). This is a proof that the Teff determination does not 
depend on the SNR of the spectra or, in other words, that the 
continuum was defined in the same way, whatever the SNR 
of the spectrum. The normalisation, as well as the tempera- 
ture determination, were done (in the so-called second step) 
on the original unsmoothed spectra. In order to further ex- 
plore the possible role of the SNR, we smoothed the spectra 
of the system 41211 10 with a 3-points boxcar (0.6 A wide) 
and performed the T^s determination again. The result was 
the same within ~ 100 K, which confirms that it is insensi- 
tive to the SNR of the spectra, since the latter, in this test. 



p. North et al.: VLT spectroscopy of eclipsing binaries in tiie SMC 



35 



was improved by a factor of V3. Therefore, it is doubtful 
that the explanation resides in the way the continuum was 
defined. 

The sky subtraction may be more problematic. Indeed, we 
have seen that the sky does vary slightly in the FLAMES 
field, and all we could do was to subtract an average sky 
spectrum from the stellar spectra. Underestimating the sky 
would have left some continuum flux that would have 
slightly diluted the stellar absorption lines, leading to an 
overestimate of Teff. The only cure to this problem would 
have been to use the IFUs instead of the MEDUSA fibres, 
in order to be able, at least in principle, to measure the sky 
surrounding each object. But this would have been at the 
expense of a drastic reduction of the number of targets. 
Although sky subtraction can explain some of the scatter 
we see on Fig. |22] it remains to understand why it could 
lead to a systematic overestimate of Tetf (causing the over- 
estimate of L and, ultimately, of the interpolated Al). Since 
an average sky was subtracted, one would expect as many 
negative values of AM as positive ones. However, the sky 
positions had been chosen, on the DSS chart, on the basis 
of their "darkness", so that the selected skies may be darker 
than the ones that typically surround the targets. 
Inadequacy of the stellar evolution models: In Figs. lSOl - HTl 
are plotted not only the evolutionary tracks for the metallic- 
ity Z = 0.004 of Charbonnel et al. (I1993I I, but also those at 
Z = 0.001 (Schaller et al. |1992| l. It is interesting to see that 
the more metal -poor models are slightly more luminous and 
hotter, and some systems fit them better It is possible that 
some stars have indeed a metallicity approching Z - 0.001 
(especially for the light elements, which are the main con- 
tributors to Z), as shown e.g. by Peters & Adelman (2006). 
However, it is doubtful that all our systems have such a low 
metallicity, and a look at the HR diagrams shows that this 
explanation cannot hold for all of them. 
Also plotted on Figs. |39|j4n are the evolutionary tracks of 
Claret & Gimenez d 19981 1 for a heHum content Y - 0.28 
and a metallicity Z - 0.004. These models have a helium 
enhancement AY - 0.028 relative to those of Charbonnel 
et al. ( I1993I I. since the latter have Y = 0.252. One can see 
that increasing the helium content increases the luminos- 
ity, withouth changing the temperature. Again, some stars 
would fit better such models than the standard ones, while 
others (like the components of 4 175333) would need so 
large a helium enhancement as to make this solution unten- 
able. 

Inadequacy of the synthetic spectra: If the true metallicity 
of most system is close to Z = 0.004 and the helium con- 
tent is normal, and if the standard evolutionary models can 
be trusted, then the synthetic spectra used to estimate the 
effective temperatures may be suspected. An overestimate 
of Teff would certainly have resulted from the use of LTE 
atmosphere models: as shown by Hunter et al. (I2007I I in 
their Table 13, the temperatures of stars analyzed by Kilian- 
Montenbruck et al. (1 19941 1. on the basis of LTE atmosphere 
models and non-LTE line formation calculations, are about 
29 000 K, instead of 26 000-27 000 K using NLTE models. 
Nevertheless, the synthetic spectra we used are precisely 



based on the latter, so that no obvious bias on tempera- 
tures is expected. Using the FASTWIND code, Massey et 
al. (I2009I I found temperatures of O stars ~ 1000 K hot- 
ter than those obtained by Bouret et al. ( 120031 ). who used 
the TLUSTY and CMFGEN codes. Part of the discrepancy 
might be explained, however, by moonlight continuum con- 
tamination in the spectra of Bouret et al. (2003). On the 
other hand, Mokiem et al. (120061 l2007b found temperatures 
~ 1000 K hotter than those of Massey et al. (120091 ), while 
they also used FASTWIND. Finally, Massey et al. (|2009] ) 
still find some amount of mass discrepancy for their LMC 
stars. Therefore, some debate still exists regarding the tem- 
perature scales, though the above references all deal with 
O stars, rather than with the B stars we are more concerned 
with. Unfortunately, similar studies for metal poor B stars 
in the range we are interested in (20 000 < T^^ < 30 000 K) 
are lacking, so it is difficult to assess the reliability of the 
temperature scale. One may expect, however, that it is more 
reliable than for O stars, since the NLTE eff'ects are less im- 
portant. 

The OSTAR2002/BSTAR2006 libraries we have used are 
recent and are probably the best available. The only reser- 
vation one could possibly raise is the choice of the micro- 
turbulence (2 km S-' for T^g < 30 000 K and 10 km s"' for 
Teff > 30 000 K), which may be a rough approximation of 
reality. 

Conversely, even perfectly realistic models may prove in- 
appropriate, if superficial abundance anomalies exist in the 
components of our binary systems. It is well known that 
early B-type stars may show a large overabundance of he- 
lium - the so-called helium-rich stars, of which HD 37776 
is the prototype - but no example is known among com- 
ponents of close binaries. In any case, stronger He lines 
would mimick a lower effective temperature, while it seems 
that the temperatures are overestimated. Thus, only helium- 
weak stars could explain an overestimate of temperatures, 
but none are known at such high masses. 

In conclusion, we consider the first explanation (sky sub- 
traction problem) as the most probable one, until proven other- 
wise by other studies. 

The astrophysical parameters of our sample of detached bi- 
naries allow us to give a mass-luminosity relation for the SMC 
valid for masses between ~4 and ~18 Mq. A linear regression 
fit gives, on all 46 points: 

log L/Lq = (3.04 ± 0.1 1) logM/Me + (0.90 ± 0.09) (29) 



with a RMS of 0.095 (see the red Une on Fig. |24]). This is al- 
most exactly the same relation as that obtained by Gonzales et 
al. 120051 for the LMC. The mass-luminosity diagram is plot- 
ted in Fig. |24] together with the isochrones computed from the 
Geneva models for Z = 0.004 and Y = 0.252 (Charbonnel et al. 
119931 ). All objects are, within the 1 o- uncertainty, on or above 
the theoretical zero-age main sequence. 

Most stars lie below the 50 Myr isochrone. The exceptions 
are some low-mass objects (< 6A40), i.e. both components of 5 
1 1 1649 and of 4 175333. These two systems appear to be >70 
Myr old. They clearly bias the regression line to a low value of 
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Table 9. Comparison with theoretical evolutionary models: dif- 
ference between the evolutionary mass and the observed mass. 



Object Model 



AMp 
(Mo) 



AMs 
(Me) 



4 110409 


sd 


1.076 


1.923 


4 113853 


sd 


1.500 


0.855 


4 117831 


d 


0.392 


0.366 


4 121084 


d 


1.503 


1.404 


4 121110 


d 


1.482 


0.923 


4 121461 


d 


0.731 


0.429 


4 159928 


sd 


1.089 


0.566 


4 160094 


d 


0.855 


0.261 


4 163552 


d 


0.222 


0.197 


4 175149 


sd 


4.143 


5.091 


4 175333 


d 


1.012 


0.797 


5 016658 


d 


0.296 


0.646 


5 026631 


sd 


0.924 


0.459 


5 032412 


d 


0.922 


0.208 


5 038089 


d 


1.764 


2.349 


5 095337 


d 


1.239 


0.934 


5 095557 


d 


1.215 


0.639 


5 100485 


d 


0.406 


0.628 


5 100731 


c 


0.248 


0.686 


5 106039 


sd 


0.807 


1.119 


5 111649 


d 


0.158 


0.479 


5 123390 


d 


1.236 


2.574 


5 180185 


d 


0.387 


0.045 


5 180576 


d 


1.119 


0.082 


5 185408 


d 


0.767 


0.939 


5 196565 


d 


0.446 


0.241 


5 261267 


sd 


1.921 


1.569 


5 265970 


d 


0.002 


0.257 


5 266015 


sd 


0.387 


2.744 


5 266131 


d 


0.637 


0.447 


5 266513 


d 


0.193 


0.458 


5 277080 


sd 


3.138 


2.331 


5 283079 


d 


0.704 


0.674 



the slope. Graczyk (I2003I I adopts the relation 

log L/Lo = (3.664 + 0.047) log M/ Me + (0.380 + 0.027) 

for the LMC and the SMC, which is shown as the full blue 
line in Fig.|24] The latter relation appears slightly too steep: it 
clearly runs above unevolved stars with Al > 8 - 10 Mq. 

5.4. Apsidal motion 

Four of the nine eccentric systems have a negligible or marginal 
(i.e < 2 cr) apsidal motion, while four others show a formally 
significant one (see Table [TTTi: the apsidal motion of the ninth 
system was arbitrarily fixed to a value close to the theoreti- 
cal expectation. Interestingly, the system 5 266131 has a very 
fast apsidal motion, and probably holds the record among those 
hosting non-degenerate components (Petrova & Orlov 19991 
Bulut & Demircan l2007l l. It is also worth noticing that in two 
of these systems (5 123390 and 5 266131), the relative radius 
of the primary component is very close to, or even slightly 
exceeds, the limiting radius above which orbital circulariza- 
tion occurs quickly (North & Zahn 12003! I2004I I: they have 
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Fig. 22. Detached systems: difference between evolutionary 
mass and observed mass vs. observed mass. Filled symbols: 
primary; open symbols: secondary. The point in the inset, to 
the lower left corner, shows the median error bars for the mass 
of the primary. The vertical bar does not include the errors on 
luminosity and Teff (hence on the interpolated mass), so it must 
be considered as a lower limit. 
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Fig. 23. Same as Fig. |22] for the detached systems of Hilditch 
et al. (i2005l l. 



rp - 0.238 and 0.281 respectively, while the limiting radius is 
~ 0.25. Thus, those systems are being caught in the act of cir- 
cularizing their orbit, and indeed their eccentricities are among 
the smallest of all our eccentric binaries. 

In order to compare our observed apsidal rotation rates with 
those predicted by theory, we used the equations (1) to (8) of 
Claret & Gimenez (,1993,) , which allow to obtain the average 
tidal-evolution constant ^aobs of the system from the orbital 
and stellar parameters, including the measured apsidal period 
^apsid - Injoj. On the other hand, we used the grid of stel- 
lar models for the metallicity Z - 0.004 computed by Claret 
(12005 D . After interpolation in surface gravity and mass of the 
theoretical \og{k2) values, the same formulae were used to ob- 
tain the predicted average constant £2theo- We corrected for the 
relativistic precession term, even though it is small for short pe- 
riods (Mazeh2008), since it contributes for no more than 2% 
to the total precession period. The result is shown in Fig. |25] 
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Fig. 24. The mass-luminosity relation for our 23 detached sys- 
tems. Full dots represent the primary components, while open 
dots represent the secondary ones. The dotted lines are the 
isochrones with age of 0, 5, 10, 20, 30, 40, 50 and 70 Myr 
for Z - 0.004 (Charbonnel et al. 1993 ). The red line is a simple 
least-squares fit of all 46 points, without weighing. The steeper 
blue line is the relation adopted by Graczyk (120031) for the LMC 
and SMC. 



-2.5 




l°g(ka)theor 

Fig. 25. Observed versus theoretical apsidal constant k2, aver- 
aged over the two companions of each binary system according 
to the prescriptions of Claret & Gimenez ( I1993I I. The internal 
structure models used are those of Claret (I2005I I for the metal- 
licity Z - 0.004. The oblique line is the equality. The error bars 
reflect the uncertainty on oj alone, so they represent lower lim- 
its to the real error. The open dot without error bar represents 
the system 5 265970, for which an arbitrary apsidal motion was 
assumed. 



where we plot £2obs ^s a function of ^itheo- The one-to-one re- 
lation is the straight diagonal line, and 3 systems fall on or very 
close to it. Two other systems are clearly discrepant. One of the 
latter, 5 95557, has a formally significant, but small apsidal mo- 
tion; the masses are not very precisely determined, which may 
explain part of the difference. It is quite possible that the apsidal 
motion of 5 95557 is affected by the gravitational perturbation 
of a small, unseen third companion. The other discrepant sys- 
tem, 5 123390. is only marginally so, since it is only 2 cr away 
from the equality line. In any case, one has to keep in mind 
that the error bars shown represent lower limits to the real un- 
certainties, because they include only the error on the apsidal 
period, while errors on the masses and radii of the components 
also contribute. The latter errors certainly dominate in the case 
of 5 266131, which lies below the equality line. The open dot 
represents the system 5 265970, for which an apsidal motion 
close to the theoretical value was assumed. 

It would be easy, by reobserving with photometry the few 
systems showing significant apsidal motion, to strongly reduce 
the error bar on the average fc2obs- Further spectroscopic data 
would be needed as well to obtain more precise masses and 
radii, in order to determine a meaningful value of £2ob,s- That 
would bring an additional constraint on the metallicity and on 
the extent of overshooting, for instance. 



5.5. Distribution of the mass ratios of the detached 
systems: do binaries really like to be twins ? 

From the 21 detached systems of the HHH03/05 sample, 
Pinsonneault & Stanek (2006) suggest that the proportion of 
massive detached systems with a mass ratio close to unity is 
far larger than what would be expected from a classic Salpeter- 
like (p(q) oc q^^-^^) or a flat (p(q) - const) ^-distribution. This 
statement is based on two arguments. The first one is the strik- 
ing diff'erence between the median mass ratio of the detached 
sample (q = 0.87) and of the semi-detached/contact sample 
iq = 0.65). The surprising point is not the low q-valne for 
the semi-detached sample, as this is expected for a post mass- 
transfer system, but rather the high proportion of systems with 
q > 0.85 in the detached sample. Moreover, they argue that 
this difference is real, i.e. not due to an observational bias such 
as the easier detection of systems with components of similar 
size and brightness. Their argument is that, beside the detached 
binaries with q > 0.85, there are two detached systems with a 
^-value as low as ~0.55. Indeed, it seems to be a very strong 
argument in favour of the reality of a population of twins, be- 
cause the observation of detached systems with such a low q 
value would mean that the cut-off value of the q distribution of 
the detached systems is small enough (< 0.55) to exclude an 
observational bias. Let us remind that the cut-off value is the q 
value which marks the transition between a double-lined and a 
single-lined binary (SB2-SB1 transition), i.e. a function of the 
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resolution and of the S /N, beside stellar parameters. Therefore, 
as our sample of 23 detached systems is comparable and even 
slightly larger than HHH03/05's, this is an excellent opportu- 
nity to compare both statistics and thus to shed some light on 
the controversial topic of the q distribution of detached and 
semi-detached/contact systems. 
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Fig. 26. Distribution of the 33 observed mass ratios (solid: 
detached, dashed: semi-detached/contact) with 0.05 bins. The 
best flat distribution is over-plotted for each category (detached 
and semi-detached or contact systems). Both distributions are 
truncated at a cut-ofF value of q ~ 0.7 and q ~ 0.4, respectively. 

The ^-distributions for our detached and semi- 
detached/contact binaries are given in Fig. |26] One sees 
that the two samples have quite different distributions. The 
semi-detached/contact distribution is patchy because of the 
small size of the sample (10 objects), but it seems to be 
compatible with a constant law. The cumulative probability 
distribution of the mass ratios (Fig. |27] top) confirms that it 
is consistent with a flat distribution, limited to the interval 
0.44 < q < 0.76. The statistics of semi-detached system is 
not discussed by Pinsonneault & Stanek (120061). The subject 
was computationally investigated by van Rensbergen et al. 
(I2005I I. Interestingly, they mention that 80% of the observed 
Algols have a mass ratio in the 0.4-1 range. The low-value 
for this range coincides with the observed cut-ofF of our small 
sample. The problem is how to interpret this value. Through 
our experience in processing spectral data, we assume in this 
paper that this is the detection threshold for SB2 systems, 
but van Rensbergen et al. seem to consider this value not as 
an observational bias but rather as a parameter of real astro- 
physical relevance. Indeed, they have reproduced the observed 
Algolid ^-distribution (from Budding et al.'s catalog 20041) 
assuming a "liberal" model of evolution, i.e. a lot of mass loss 
and a httle loss of angular momentum during mass transfer. 



Their simulations clearly depict a distribution with a broad 
peak in the 0.4-0.6 range and decreasing for higher (7-values 
(see Fig. 5 in van Rensbergen et al. l2005i) . That is in excellent 
agreement with our results too. The only problem is one of 
their initial assumptions: they consider that the ^-distribution 
of the detached Algol-progenitor systems obeys a Salpeter-like 
power law! That is not what we observe with our sample of 
23 detached binaries. Figures |26] and |27] (top) show that the 
^-distribution of these systems is hardly compatible with a 
decreasing power law. Indeed, the data are very well modeled 
by a flat distribution truncated at q = 0.72. The higher value 
for the SB1-SB2 cut-ofF compared to semi-detached systems is 
related to the almost ^^ -dependence of the luminosity ratio for 
main-sequence stars. Let us mention that a flat ^-distribution 
is assumed by most population synthesis studies (Pinsonneault 
& Stanek .2006j . 



P O.B - 




Fig. 27. Top: cumulative probability distribution as a function 
of the mass ratio for our 23 detached (solid line) and 10 semi- 
detached/contact (dotted line) binaries. The dashed oblique line 
plotted over the detached distribution represents a flat distri- 
bution with a cutoff at ^ = 0.72. In the case of the semi- 
detached/contact systems, a flat distribution with a cutoff at 
q = 0.43 provides a reasonable fit as well. Bottom: cumulative 
probability distribution for HHH03/05's 21 detached binaries. 
The over-plotted flat distribution with a cutoff at q - 0.5 shows 
the incompatibility of a flat distribution with the HHH03/05 
data. There is a similar figure in Pinsonneault & Stanek (I2006I I. 
They show that the data are best-fitted with a 55% flat-45% 
twin distribution (not shown in this figure). 



The cumulative probability distribution of the mass ratios 
for the HHH03/05 data is given in Fig. [27] (bottom). A simi- 
lar figure is shown in Pinsonneault & Stanek. This distribution 
differs significantly from ours on two points: the low cut-off 
value of ~0.5 and the relatively strong deviation from a flat 
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distribution. The low cut-off is rather surprising. The mean ap- 
parent magnitude of the HHH03/05 sample (/i - 15.47 mag) 
is 1.38 mag lower than ours (/i = 16.85 mag). A rough esti- 
mation of the ratio of our mean S /N to their mean S /N, taking 
into account our 2-times bigger telescope and longer exposure 
time via Eq. 3.2 of Hilditch ('2001), gives a ratio of ~0.8. That 
means that our mean spectral S /N is 20% smaller than theirs 
despite the bigger optics because of significantly fainter targets. 
Nevertheless, their 20%-better 5/A^ hardly explains a luminos- 
ity ratio cut-off of ~l/8 to be compared with our ~l/3-value 
(for detached main-sequence binaries, X. ~ q^)- Moreover, the 
instrumental resolution of the LR2 setup in MEDUSA mode is 
3-times the value of the 2dF spectrograph with the 1200B grat- 
ing. Therefore, either our sample is deficient in low-^ systems 
or there are some badly determined ^-values in one of the sam- 
ples. The first explanation cannot be entirely ruled out because 
of the small-sample statistics. Nevertheless, the comparison of 
the RV-related parameters for the five systems common to both 
HHH03/05's program and ours can certainly help to settle the 
issue. 

Can the "twin hypothesis" be excluded by our data? The 
thorough study by Lucy (I2006I I points to a negative answer He 
discusses both the "weak" hypothesis CKw) of an excess of bi- 
nary systems with q > 0.8 relative to a constant distribution, 
and the "strong" hypothesis CHs) of an excess of systems with 
q > 0.95. He shows that 'Hs is confirmed for binary systems 
in our Galaxy with an accurate enough mass ratio (o-g < 0.01). 
Therefore, it would be of great interest to test the same hypoth- 
esis in another galaxy, and our study might, at first sight, be 
considered a second step towards that goal, after HHH03/05. 
But a realistic typical error on our q values is ~ 0.05, and our 
sample of 15 systems with q > 0.84 is so small, relative to the 
102 systems used by Lucy, as to make impossible any confir- 
mation or rejection of 1-1 s- Our data might be more useful to 
constrain '7/vv, but extensive simulations would be necessary to 
estimate the detection biases, and such an effort does not appear 
justified by our small sample. Therefore, the twin hypothesis 
remains an open question as far as the SMC (and any galaxy 
other than ours) is concerned. 

5.6. Extinction and distances 

5.6.1 . Colour excess and extinction 

The histogram of the 28 colour excesses is given in Fig. |28] 
Note that one system (5 196565) does not have B and V light 
curves and four systems have an unreliable observed B - V in- 
dex (see Section [33] ). The mean value is Eb-v = 0.134 + 0.051 
mag with individual values in the range 0.052-0.252 mag. By 
comparison, for four stars in the bar, Gordon et al. ( 2003! ) 
found values ranging from 0.147 to 0.218 mag. From their 
extinction map across the SMC, Zaritsky et al. (I2002I I give 
Eg-v ~0. 05-0. 25 mag. Sasselov et al. ( 119971 1 found a mean 
value of Eb-v - 0.125 + 0.009 mag. Thus, our results are in 
very good agreement with these previous determinations. In or- 
der to investigate the spatial variation of the extinction across 
the studied field, Eb-v was plotted against 6, the distance to the 
optical center along the projected axis (not shown). No correla- 



tion was found. The GIRAFFE field is likely too small to detect 
a trend, if any. 
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Fig. 28. Histogram of colour excesses of 32 binaries in 0.025 
bins. The best Gaussian fit is over-plotted. 



5.6.2. Distance 

The parameters relevant to the determination of the distance 
moduli are shared between Tables |4] (synthetic photometry) 
and [15] The DM was computed in both V- and /-band, as 
explained in Section [3.12l One system (5 196565) has no value 
because of only one light curve available. Fortunately, for most 
systems the difference between the two computations of the 
individual DM is low and well within the uncertainties. As ex- 
pected in Section [3.12l the /-values have a slightly lower abso- 
lute uncertainty. The discrepancy between the V and / moduli 
is negligible: Av-i - -0.01 mag with a scatter cr = 0.02 mag 
for the 28 reliable systems. 

The DMi were checked for a possible dependence upon a 
number of parameters: apparent visual magnitude, absolute vi- 
sual magnitude, systemic velocity, colour excess and distance 
to the optical center along the projected axis. Statistical tests 
were performed (calculus of the Spearman rank correlation co- 
efficient, followed by a Student T test) and no significant corre- 
lation was found. The DMv vs. My and DMv vs. V diagrams 
are shown in Fig. [29] (bottom). In order to calculate an unbi- 
ased mean DM, the HR diagrams (Fig. [39pn i of the individual 
systems were used to select the systems suitable for this pur- 
pose. This is actually a more objective criterion than relying 
solely on the quality of the light curves and spectra. Therefore, 
8 objects showing strong discrepancies with theoretical models 
were conservatively discarded from the sample of 27 binaries 
(which remain from the initial sample of 33 binaries, after sup- 
pression of stars with unreliable or no colour index, or with a 
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3rd light, see below). Some of the primaries of these 8 systems 
may have luminosities really higher or lower than expected, but 
we can lucidly assume that most discrepancies are the direct 
consequences of low S /N data and/or third light contamina- 
tion. These suspicious systems are indicated by open symbols 
inFig.|29] 

From the sample of 33 objects, six at least must be dis- 
carded for the computation of the mean distance modulus: 5 
196565 because of the total lack of B and V data, 4 163552 be- 
cause of a third light contribution, and 4 systems which do not 
have reliable observed colour indices (5 038089, 5 180185, 5 
261267 and 5 277080). At that stage, there remain 27 systems 
with presumably good distance moduli, without consideration 
of their HR diagrams. Removing further the 8 systems which 
have a dubious position in the HR diagram (i.e. 4 113853, 4 
121084, 4 121110, 4 175149, 4 175333, 5 95557, 5 123390 
and 5 185408) leaves us with only 19 systems, a bit more than 
half of the initial sample. By "dubious position", we mean sys- 
tems for which the error bar of the mass does not overlap that 
of the luminosity/temperature on the HR diagram. A third sub- 
sample was defined, which contains the nine "best" systems, 
defined as those detached systems for which, either the lumi- 
nosity/temperature error bar overlaps the Z = 0.004 evolution- 
ary track, or the mass error bar overlaps the representative point 
of the corresponding component, for at least one of the com- 
ponents. These systems are 4 121461, 4 160094, 5 16658, 5 
32412,5 100731,5 111649, 5 265970,5 266131 and5 266513. 

Another important point is to identify possible biases act- 
ing upon the DM distribution. Actually, there is a small mag- 
nitude cut-off. This is the apparent magnitude cut-off, close to 
y ~ 18.1 mag. It almost coincides with the apparent magni- 
tude of 5 266513, the faintest binary of our sample, whose spec- 
tra have (smoothed) S /N in the 18-50 range. This V cut-ofFsets 
a higher limit for the absolute magnitude a binary should have 
in order to be seen across the whole depth of the SMC field. If 
we take ~19.5 as the upper limit for the DM, then we find that 
the M^cut-off^ is-1.82±0.16 mag with £b_ v = 0. 1 34 + 0.05 1 
mag. The cut-off strip for yf = 18.1 mag is indicated in Fig. |29] 
(bottom-left). Systems to the right of the oblique lines are too 
faint to be observed as SB2 with good enough spectra. 

For example, in spite of a reliable solution, the binary 4 
175333 should be discarded from the calculus of the mean 
DM in order to remove the observational bias just discussed. 
Indeed, with M^ = -1.217 and (y-Mv)o = 18.61, this system 
would not be perceived as an SB2 if it were more distant, i.e. if 
it had (V - Mv)o > 19.0 (see Fig.|29l bottom-left). In fact this 
object had already been discarded on the basis of the peculiar 
position of its components in the HR diagram. 

At the other end of the My- spectrum, there is a second 
cut-off, but this one does not bias the DM-distribution. This is 
the cut-off of the luminosity function, responsible for the de- 
pleted left-part of Fig.^(bottot7i- right). Close to MJ 4.2 

mag, this is roughly the low-limit for late O stars: intrinsically 
brighter stars should be earlier O stars, but they are very scarce, 
so the probability of finding one in a sample of 33 objects is 
very low. 

The histograms of the individual DMi in V and / are shown 
in Fig. |29] (top-left). For a given band, both the whole sam- 



Table 10. Summary of the average and median distance mod- 
uli, for the three subsamples defined in the text. The given un- 
certainties represent the RMS scatters, not the errors on the re- 
spective average distance moduli. 



Subsample 


Stat, type 


(V - Mv)o 


(/ - Mi)o 


all 27 


w. average 


19.07 ±0.19 


19.07 ±0.19 




median 


19.11 ±0.12 


19.11 ±0.12 


19 good 


w. average 


19.10 ±0.14 


19.11 ±0.14 




median 


19.13 ±0.10 


19.12±0.13 


9 best 


w. average 


19.10 ±0.15 


19.12 ±0.16 




median 


19.13 ±0.07 


19.12 ±0.07 



pie (i.e. the 27 systems remaining after the exclusion of the 
obviously problematic ones) and the "good systems" (i.e. the 
19 best ones) distributions are plotted. Although 8 systems are 
considered as less reliable, one can see that all but three (4 
175149, 4 175333 and 5 123390) discarded values are within 
the "good systems" range of DMi. This range spans 0.45 mag, 
from 18.87 to 19.32. The modes of the V and / distributions of 
the "good systems" are 19.15 and 19.20 respectively. In addi- 
tion the distribution of the 9 "best" systems is plotted in red. 

One can wonder how far the average value of the DMs is the 
most pertinent indicator of the mean distance of our stars. There 
is a ~ 0.05 to 0.1 mag deviation between the modal and mean 
values of the distributions (for both V and /), but we verified 
that these distributions, in spite of their apparent asymmetry, 
do not depart from gaussians in a significant way. Since the 
modal value is not a robust one for such small samples, we 
rather considered the weighted average and the median values. 

The average and median distance moduli for the three sub- 
samples defined above are given in Table [10] for the V and / 
bands. The average values are weighted by the reciprocal vari- 
ance of the individual moduli. The RMS scatter is given be- 
side each average value; beside each median value is given the 
width of the distribution, defined as the half difference between 
the upper and lower quartiles. 

We see that the V and / distance moduli of any given sub- 
sample are perfectly consistent. More surprisingly, the sample 
of all 27 systems has a slightly shorter distance modulus than 
the two other subsamples, while it contains many systems with 
large positive AM values (see Section l53] l. Likewise, the sam- 
ple of the 9 best systems (i.e. those which match the standard 
Z = 0.004 evolutionary tracks best) have the same distance 
modulus as the sample of the 19 good systems. These results 
suggest that the components which are overluminous with re- 
spect to their Z = 0.004 track are not necessarily overluminous 
in reality. Otherwise, if their luminosity were overestimated, 
they should also lie at a larger apparent distance. 

Finally, combining both V and / results, we adopt 

<(m - M)o) = 19.1 1 + 0.03 (66.4 ± 0.9 kpc) (30) 

as the most probable average distance modulus for our sample. 
That value is 0.20 + 0.03 mag (5.8 + 0.9 kpc) higher than the 
mean DM derived by HHH05, that is to say 18.912 ± 0.035 
(s.e.) (60.59 + 0.98 kpc) from two samples totalizing 39 bi- 
naries. Our DM value does not agree HHH05's, but contrary to 



p. North et al.: VLT spectroscopy of eclipsing binaries in tiie SMC 



41 



_ 1 1 1 1 1 1 


""\ 




- 


_ 1 1 1 1 1 1 

^N-27 

^N=19 

- N=g 


1 

^ - 


- 


\ 


\ - 




\ 




- •• \^' 










- , 1 \ \y , 1 




V, , - 



19 



19 
DM, 



2° covered by HHHOS's sample! Such a huge distance gradient, 
if real, would have been detected by Groenewegen l2000l Thus 
our larger distance cannot be reconciled with that obtained by 
HHH05. It is most probably due to the better observational ma- 
terial used here, as well as to the excellent synthetic spectra of 
Lanz & Hubeny (l2003ll2007l l. On the other hand, our DM per- 
fectly agrees with that obtained by Groenewegen ( I2000I I using 
cepheids (19.1 1 ±0.1 1 from the W index, and 19.04+0.17 from 
the K band). Taken at face value, our distance modulus clearly 
favors the so-called long distance scale. 
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Fig. 29. Top: distribution of the V- {left) and /-band {right) 
distance moduli. The distribution of the whole sample of 27 
systems with good colours and presumably without 3rd Ught 
(dashed line) is shown. Only binaries with good HR diagrams 
are selected for the smaller statistics. The best fitting Gaussian 
curves for the 19-binary distributions are over-plotted. Bottom: 
distance modulus vs. absolute visual magnitude {left) and dis- 
tance modulus vs. apparent visual magnitude {right). The es- 
timated magnitude cut-off for a mean extinction is indicated 
by the diagonal solid line; the dashed lines correspond to the 
cut-ofF lines for Eb-v - 0.134 ± 0.051 mag. Ten binaries with 
dubious luminosities are shown with open symbols. 



theirs, it is not representative of the whole SMC. Our DM refers 
to the 25'-diameter GIRAFFE field at -0.45° of the SMC op- 
tical center in a SW direction. Is this result coherent with the 
current knowledge of the 3-D structure of the SMC? The SMC 
bar is currently thought to be NE to SW-oriented, with the NE 
part closer to us than the SW part and the optical center occu- 
pying a mid position (Groenewegen l2000l HHH05). According 
to HHH05, there is a ~2-3 kpc variation of the mean DM over 
their 2° field. Given that gradient, our sample can be expected 
to be ~0.5 kpc further than the optical center, since our field is 
~0.45°to the SW from it; but this difference is actually smaller 
than the uncertainty on our mean distance. If we took at face 
value the average distances moduli obtained by HHH05 and 
by us respectively, and assumed that HHH05 give a distance 
modulus representative of the SMC optical centre, one could in 
principle compute the real distance between the optical centre 
and the region corresponding to our field of view. Thus, taking 
into account the angular separation between the optical center 
and the center of the GIRAFFE field of view, the separation 
would correspond to 5.8 ±1.3 kpc along the SMC bar (see Fig. 
[TJ. That would mean that the bar is almost directed towards us, 
with a mean distance varying by as much as 25.8 kpc over the 



5.6.3. Distance according to tine detaclied systems 
versus distance according to tine otiiers 

In the context of a search for systematic effects, one may won- 
der whether detached systems give a different distance modu- 
lus than the semi-detached or contact ones. Separating our 19 
"good" systems into a group of 13 detached ones and a group 
of 6 semi-detached ones, we have computed the weighted av- 
erage distance modulus (in the / band) for each group. For the 
detached systems, we obtain (/ - Mj)o{d) - 19.07 ± 0.04 mag, 
the error being the RMS scatter divided by VT3. For the 
semi-detached and contact ones, we obtain (/ - Mi)o{sd, c) - 
19.15 ± 0.05 mag. Thus, our 6 semi-detached and contact sys- 
tems are, on average, 0.08 mag more remote than our 13 de- 
tached systems. Given the small statistics and the error bars, 
the difference is not significant since it amounts to only about 
1.6 cr. 

Nevertheless, we have looked at the distance moduli de- 
termined by HHH03/05, in order to see whether a similar dif- 
ference exists in their data. Interestingly, almost the same dif- 
ference appears: for 21 detached systems, one has DM{d) - 
18.85 ± 0.06, while for the 29 others one has DM{sd, c) = 
18.91 ± 0.05. So, here too, the semi-detached and contact sys- 
tems are more remote than the detached ones on average, and 
by a similar amount, of 0.06 mag. The effect is only of the or- 
der of 1 cr, so it is formally not significant. The coincidence be- 
tween our results and those of HHH03/05 remains intriguing, 
however 



5.6.4. Deptil 

The observed dispersion of the DMi is the convolution of 
the true depth of our population of binaries with the aver- 
age uncertainty on the DM determined for one binary sys- 
tem. Hence, under the assumption of a Gaussian distribution 
of both the cosmic and the error standard deviations, the true 
1 cr depth is given by the quadratic difference between the ob- 
served standard deviation of the moduli and the typical error 
on the modulus of an individual system. The observed stan- 
dard deviation being 0.144 mag, and taking the average of the 
squares of the estimated errors as representative of the typi- 
cal eiTor, the 1 cr depth becomes VO.1442 - 0.00479 ^ o.l26 
mag RMS (~3.9 kpc), since the average variance of an individ- 
ual / distance modulus is 0.00479 mag^. Groenewegen gives 
0.11 mag (3.4 kpc) as the intrinsic 1 cr spread for his sample of 
SMC Cepheids, so our result is in perfect agreement with his. 
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The HHH05 sample (29 binaries) displays a higher dispersion, 
~0.28 mag (~7.8 kpc), but the authors consider that value as an 
overestimate. Our depth estimate is only 15% larger than that 
of Groenewegen (120008 . 

5.6.5. Temperature scale and distance modulus 

We have seen that a suspicion remains, that the effective tem- 
peratures may have been overestimated in a number of systems. 
What is the effect of such an overestimate on the resulting dis- 
tance modulus? For the sake of simplicity, let us consider a sin- 
gle star rather than a binary system. Its distance modulus can 
be written (under the assumption that we use the V magnitude): 

m-M = V,, + 2.5 log[lO-^-^'^'')- 3.1 Eb-v 
= const + 2.5 log(lO-°-^<^"'"-*^>) 
-3.1 [(B-V)g-(B-V)Q\ 

since Vg - const. Invoking Eqns. [T] and [16] which allow to 
express (B - V)o and BC as a function of Teff, expressing the 
bolometric magnitude as a function of radius R and Tetf , and 
remembering that (B - V)q is constant for a given object, one 
obtains 

m-M ^ const +2.5 log (l0-°-4[-^'°s«-'°'°s^-+='^^""°s7-.r]-) 

-1.395 log Teff 
= const +3.095 log TefF 

where we have taken into account that the radius is constant. 
Introducing a second star will essentially not change the result, 
because the ratio of the temperatures will appear (as well as the 
ratio of radii), which can be considered constant. Therefore, by 
differentiating the last equation above, one finally gets a very 
simple relation between a change AT^ff of effective tempera- 
ture, and the resulting change of the distance modulus: 



A(m-M)o- 1.344 



Areff 



(31) 



'eS 



In the case of 4 175333, for instance, the temperature excess 
is especially large since it reaches ~ 2000 K. If the tempera- 
ture is indeed overestimated, then the system is A(m - M)o - 
1.344 X 0.092 = 0.123 mag too far Strangely enough, its mod- 
ulus is only 18.60, so that it would be reduced to 18.48 under 
this assumption. The system 4 121110 has AT^eff ~ 1580 K and 
A(m - M)o ~ 0.08, so its modulus should be 19.00 instead of 
19.08. Likewise, 4 121461 and 4 160094 have A(m - M)q ~ 
0.08 and 0.06 respecively, while they are listed as "good" sys- 
tems. 

A quick estimate of the average temperature excess of the 
9 "best" systems leads to an average modulus excess of about 
0.04 mag, so that the corrected distance modulus would be re- 
duced to 19.07, which still remains rather long compared to the 
frequently quoted value of 18.9. In any case, it appears diffi- 
cult to reduce the distance modulus to less than 19.0. Thus, our 
data do favour the "long" rather than the "short" distance scale 
to the Magellanic Clouds. 



6. Conclusions 

The first goal of the present paper was to extract the fundamen- 
tal stellar parameters characterizing 33 SB2 systems and then 
to test the reliability of stellar evolution models at low metal- 
licity. We have shown that most components of detached sys- 
tems have properties in fair agreement with the predictions of 
the Geneva single-star models at Z = 0.004, although several 
systems appear slightly overluminous for their mass. From the 
subsample of detached systems, we give a mass-luminosity re- 
lation vaUd for masses between ~ 4 and ~ 18A^0. Despite the 
lack of evolutionary cross checks, the semi-detached binaries 
appear to satisfy the typical morphology of Algol-like systems: 
a main-sequence primary with a large, less massive and more 
evolved secondary companion. 

Our binaries were then used as primary standard candles. 
Average distance moduli obtained from the /-band data show 
no systematic difference compared to values obtained from the 
V data. For each band, the frequency distribution of the DMi is 
not significantly different from a Gaussian. Combining V and / 
results and considering both the mean and the median values, 
we adopt {(m - M)o) = 19.11 + 0.03 as the average DM for the 
19 most reliable systems. We have shown that, even if the ef- 
fective temperatures are overestimated, so as to make the com- 
ponents of our systems overluminous relative to their mass, the 
average DM would still remain above 19.00, so that our results 
support rather the long distance scale to the Magellanic Clouds 
than the short one. The true depth of our sample is estimated 
to be 2 X 0.126 - 0.25 mag (~7.6 kpc), if we consider the 2 cr 
interval as representative. 

We have found a significant apsidal motion for at least three 
systems. For one system, the apsidal motion might be affected 
by an unseen third companion, but at least two other systems 
most probably owe it to the stellar structure of their compo- 
nents alone. The observed and predicted precession rates agree 
perfectly for the latter two systems, which are now circulariz- 
ing their orbit, and the agreement extends to one or two other 
systems showing much less precise apsidal motion. The system 
5 266131 shows probably the largest precession motion of all 
systems hosting main sequence stars, which is due exclusively 
to the non-point mass nature of the components. 

Beside the astrophysical parameters and the distance mod- 
uli, other interesting results have been obtained. We have 
shown that the twin hypothesis for massive binaries is not sup- 
ported by our results, although the latter are not strong enough 
to reject it either Our sample of 23 detached systems agree 
quite well with a standard flat mass ratio spectrum. There is 
a clear cut-off of the ^-distribution at ~ 0.72. Therefore, the 
high proportion of detached binaries having a q value close to 
unity is more probably the result of an observational bias. The 
discrepancy observed by Pinnsonneault & Stanek (I2006I I may 
well be due to a number of underestimated mass ratios in the 
HHH03/05 sample. The distribution of the mass ratios of the 
semi-detached/contact binaries is consistent with a flat distri- 
bution restricted to the interval 0.4 < q < 0.8. The sample is 
too small (10 objects) to tell more. 
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Another important result is the bimodal distribution of the 
radial velocity of the Hii gas, with peak values at 145.2 + 0.5 
and 161.4 ±0.6 km s'. 

Finally, we must emphasize that the accomplishment 
of this project was possible thanks to the public avail- 
ability of the OGLE-II catalog of EB light curves, the 
OSTAR2002/BSTAR2006 Ubraries of synthetic spectra, the 
KOREL code and the PHOEBE interface to the WD codes. 
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Table 11. Orbital parameters: number of spectra used in the solution, minimum and maximum signal-to-noise ratio of the 
smoothed spectra, eccentricity, argument of the periastron for the epoch fo given in Table [3] and first derivative of it, inclina- 
tion, semi-major axis, mass ratio, phase shift and systemic velocity. 
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0.767 ± 0.006 


0.00014 ±0.00019 


177.4 ±0.7 


5 038089 


11 


128- 


-216 





90 


77.1 ±0.2 


21.88 ±0.16 


0.899 ±0.011 


-0.00047 ± 0.00039 


152.9 ± 1.2 


5 095337 


10 


48- 


-89 





90 


88.8 ±0.4 


9.98 ± 0.08 


0.874 ±0.013 


0.00278 ± 0.00033 


183.8 ± 1.8 


5 095557 


11 


33- 


-75 


0.218 
±0.003 


155.7 ±2.5 
2.6 ±0.5 


84.7 ± 0.3 


17.54 ±0.32 


0.854 ± 0.026 


0.00837 ± 0.00062 


203.3 ± 2.6 


5 100485 


13 


46- 


-86 





90 


84.3 ± 0.2 


13.12 ±0.04 


0.972 ± 0.006 


-0.00109 ± 0.00026 


159.2 ±0.6 


5 100731 


8 


35- 


-74 





90 


59.7 ± 1.1 


11.36±0.16 


0.762 ±0.012 


-0.01465 ± 0.00173 


173.0 ±1.6 


5 106039 


9 


52- 


102 





90 


75.7 ±0.1 


16.81 ±0.21 


0.527 ± 0.012 


-0.00182 ± 0.00041 


153.3 ±2.1 


5 111649 


10 


52- 


-89 





90 


74.4 ± 0.2 


19.14±0.11 


0.936 ±0.010 


-0.00162 ± 0.00040 


149.4 ±0.8 


5 123390 


14 


91- 


179 


0.042 
±0.006 


62.3 ± 5.2 
4.8 ± 1.6 


74.3 ± 0.2 


18.39 ±0.19 


0.804 ±0.010 


0.00470 ± 0.00046 


141.7 ± 1.8 


5 180185 


12 


37- 


-64 


0.036 
±0.005 


88.6 ±2.8 
0.0 (fixed) 


84.9 ±0.2 


29.91 ±0.14 


0.901 ± 0.009 


-0.00182 ±0.00051 


165.2 ±0.6 


5 180576 


12 


25- 


-64 





90 


75.0 ±0.4 


13.31 ±0.08 


0.784 ± 0.009 


0.00350 ± 0.00079 


168.1 ± 1.0 


5 185408 


12 


38- 


-66 





90 


76.2 ±0.3 


12.69 ±0.06 


0.849 ± 0.008 


-0.00141 ± 0.00057 


174.9 ± 0.8 


5 196565 


13 


44- 


-86 


0.138 
±0.006 


231.9 ±2.7 
1.2 ±0.6 


83.5 ± 0.2 


24.65 ±0.12 


0.803 ± 0.007 


-0.03068 ± 0.00038 


162.3 ± 0.5 


5 261267 


10 


41- 


-88 





90 


89.7 ± 0.2 


11.92 ±0.08 


0.525 ± 0.005 


0.00289 ±0.00018 


176.0 ± 1.3 


5 265970 


10 


67- 


HI 


0.054 
±0.001 


362.5 ± 4.7 
3.8 (fixed) 


82.5 ±0.1 


24.67 ± 0.24 


0.620 ± 0.008 


0.01600 ±0.00030 


159.5 ± 1.3 


5 266015 


10 


69- 


118 





90 


78.5 ±0.1 


17.60 ±0.09 


0.441 ± 0.005 


-0.00215 ± 0.00028 


188.0 ± 1.1 


5 266131 


10 


37- 


-71 


0.042 
±0.001 


227.2 ± 2.3 
50.2 ± 2.0 


83.5 ± 0.2 


12.82 ±0.15 


0.864 ±0.013 


-0.00962 ± 0.00028 


183.1 ±2.6 


5 266513 


10 


18- 


-50 





90 


84.0 ± 0.4 


9.91 ± 0.06 


0.845 ± 0.010 


-0.00073 ± 0.00056 


164.4 ± 1.1 


5 277080 


11 


74- 


•130 





90 


76.6 ±0.1 


15.96 ±0.13 


0.504 ± 0.009 


-0.00020 ± 0.00034 


155.7 ± 1.3 


5 283079 


10 


28- 


-68 





90 


87.7 ± 0.3 


12.28 ± 0.08 


1.003 ± 0.013 


0.00232 ± 0.00030 


179.4 ± 1.3 
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Table 12. Wilson-Devinney analysis, as obtained on the basis of both photometry and spectroscopy: system type (detached, semi- 
detached or contact), temperature ratio, normalized surface potentials and passband luminosity ratios. Third light contribution is 
given for the system 4 163552 (normalized to total flux at phase 0.25). The given potentials and the luminosity ratios are those 
obtained including the spectroscopic constraints (final WD/PHOEBE analysis). 



Object 


Model 


T^JT^., 


n. 


fis 




(Ls/Lp)v 

ti.v 


(Ls/ip); 


iUIL^f^ 


4 110409 


sd 


0.654 ± 0.007 


5.732 ± 0.092 


2.922 (fixed) 


1.292 ±0.023 


1.353 ±0.024 


1.367 ±0.024 


1.392 ±0.036 


4 113853 


sd 


0.737 ±0.015 


3.868 ± 0.090 


3.209 (fixed) 


0.680 ±0.010 


0.713 ± 0.009 


0.753 ± 0.009 


0.681 ±0.058 


4 117831 


d 


0.981 ±0.023 


4.855 (fixed) 


4.479 ±0.108 


1.175 ±0.052 


1.177 ±0.050 


1.180 ±0.049 


1.009 ±0.103 


4 121084 


d 


0.940 ± 0.007 


3.770 (fixed) 


3.817 ±0.043 


0.707 ±0.061 


0.713 ±0.010 


0.720 ±0.010 


0.704 ±0.061 


4 121110 


d 


0.894 ± 0.007 


3.820 ± 0.035 


4.184 ±0.059 


0.424 ± 0.005 


0.433 ± 0.005 


0.442 ± 0.005 


0.415 ±0.047 


4 121461 


d 


0.949 ± 0.023 


6.980 ±0.133 


6.822 ± 0.099 


0.907 ± 0.030 


0.914 ±0.029 


0.923 ± 0.029 


0.954 ± 0.063 


4 159928 


sd 


0.723 ± 0.010 


3.138 ±0.024 


2.983 (fixed) 


0.357 ± 0.003 


0.376 ± 0.003 


0.397 ± 0.004 


0.453 ± 0.029 


4 160094 


d 


0.854 ±0.018 


5.517 (fixed) 


5.545 ±0.100 


0.476 ± 0.050 


0.489 ± 0.007 


0.502 ± 0.007 


0.428 ±0.101 


4 163552 


d 


0.997 ± 0.004 


3.983 ± 0.026 


4.214 ± 0.044 


0.798 ± 0.024 
0.142 ±0.014 


0.799 ±0.016 
0.112 ±0.012 


0.800 ±0.014 
0.1 19 ±0.006 


0.926 ±0.171 


4 175149 


sd 


0.849 ± 0.006 


4.411 ±0.056 


3.171 (fixed) 


1.176 ±0.025 


1.200 ±0.023 


1.229 ±0.023 


1.415 ±0.132 


4 175333 


d 


0.879 ± 0.022 


5.054 ±0.151 


5.010 ±0.196 


0.582 ± 0.034 


0.595 ± 0.035 


0.609 ± 0.036 


0.562 ±0.103 


5 016658 


d 


1.007 ±0.009 


4.086 (fixed) 


4.590 ± 0.068 


0.609 ± 0.074 


0.612 ± 0.008 


0.614 ± 0.008 


0.604 ± 0.074 


5 026631 


sd 


0.734 ± 0.007 


3.507 ± 0.030 


3.226 (fixed) 


0.495 ± 0.003 


0.519 ±0.003 


0.547 ± 0.003 


0.736 ± 0.049 


5 032412 


d 


0.880 ± 0.008 


7.284 (fixed) 


6.949 ± 0.073 


0.578 ± 0.050 


0.584 ± 0.007 


0.583 ± 0.007 


0.546 ± 0.005 


5 038089 


d 


1.012 ±0.009 


4.606 ± 0.047 


5.065 ± 0.067 


0.693 ±0.019 


0.692 ±0.018 


0.691 ±0.018 


0.678 ± 0.032 


5 095337 


d 


0.936 ± 0.007 


3.685 ± 0.022 


3.835 ± 0.040 


0.676 ± 0.008 


0.684 ± 0.008 


0.691 ±0.008 


0.659 ± 0.085 


5 095557 


d 


0.908 ±0.013 


5.991 ±0.107 


6.602 ±0.165 


0.523 ± 0.015 


0.532 ±0.015 


0.541 ±0.015 


0.625 ±0.190 


5 100485 


d 


1.008 ±0.007 


5.193 ±0.076 


5.216 ±0.080 


0.966 ± 0.034 


0.965 ± 0.034 


0.963 ± 0.034 


0.930 ± 0.052 


5 100731 


c 


0.828 ± 0.035 


3.851 ±0.058 


3.351 (fixed) 


0.517 ±0.014 


0.533 ±0.013 


0.552 ±0.012 


0.438 ± 0.041 


5 106039 


sd 


0.638 ± 0.006 


5.178 ±0.057 


2.928 (fixed) 


1.007 ±0.009 


1.076 ±0.009 


1.168 ±0.009 


1.062 ±0.063 


5 111649 


d 


1.025 ±0.011 


4.600 (fixed) 


4.736 ± 0.067 


0.889 ± 0.056 


0.886 ±0.015 


0.882 ±0.014 


0.880 ± 0.052 


5 123390 


d 


1.017 ±0.020 


5.077 (fixed) 


5.631 ±0.106 


0.587 ± 0.060 


0.587 ±0.011 


0.585 ±0.010 


0.584 ± 0.061 


5 180185 


d 


1.036 ±0.029 


7.724 ±0.169 


8.571 ±0.249 


0.707 ± 0.042 


0.705 ± 0.025 


0.700 ± 0.024 


0.611 ±0.042 


5 180576 


d 


0.813 ±0.027 


4.918 ±0.170 


5.071 ±0.207 


0.450 ± 0.027 


0.465 ± 0.028 


0.483 ± 0.029 


0.45 ± 0.07 


5 185408 


d 


0.944 ±0.018 


5.052 (fixed) 


4.957 ± 0.094 


0.757 ± 0.047 


0.764 ±0.017 


0.772 ±0.016 


0.619 ±0.045 


5 196565 


d 


0.867 ±0.016 


5.956 (fixed) 


6.639 ±0.120 


0.43 ± 0.04" 


b 


0.453 ± 0.009 


0.413 ±0.040 


5 261267 


sd 


0.728 ± 0.004 


3.353 ±0.013 


2.923 (fixed) 


0.439 ± 0.002 


0.461 ±0.003 


0.487 ± 0.001 


0.451 ±0.039 


5 265970 


d 


0.901 ±0.012 


4.157 (fixed) 


5.370 ± 0.056 


0.219 ± 0.024 


0.264 ± 0.001 


0.227 ± 0.001 


0.233 ± 0.024 


5 266015 


sd 


0.757 ± 0.005 


3.512 ±0.026 


2.759 (fixed) 


0.494 ± 0.002 


0.515 ±0.002 


0.539 ± 0.002 


0.617 ±0.045 


5 266131 


d 


0.936 ± 0.007 


4.510 (fixed) 


4.727 ± 0.065 


0.647 ± 0.060 


0.655 ± 0.008 


0.661 ± 0.007 


0.644 ± 0.060 


5 266513 


d 


0.959 ±0.018 


4.390 (fixed) 


4.530 ± 0.084 


0.692 ± 0.070 


0.697 ± 0.016 


0.703 ± 0.015 


0.733 ± 0.076 


5 277080 


sd 


0.685 ± 0.004 


3.932 ± 0.047 


2.884 (fixed) 


0.567 ± 0.015 


0.600 ± 0.008 


0.641 ± 0.008 


0.690 ± 0.022 


5 283079 


d 


0.997 ± 0.010 


5.162 (fixed) 


5.184 ±0.074 


0.987 ± 0.086 


0.987 ±0.015 


0.988 ±0.015 


0.936 ± 0.088 



" B data missing, ratio given by the WD code on the basis of the / data, via the atmosphere models 
'' V data missing 
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Table 13. Astrophysical parameters for the primary components: masses, radii, effective temperatures, surface gravities, lumi- 
nosities and equilibrium equatorial velocity (i.e. synchronous velocity for circular systems and pseudo-synchronous velocity for 
eccentric systems). 



Object 


Mp 


Rp 


P 

■'cff 


log^p" 


log Lf/L^ 


V'' 




(Me) 


(Ro) 


(K) 


(dex cgs) 




(kms"') 


4 110409 


13.93 ±0.59 


4.64 ±0.14 


32370 ±816 


4.25 ± 0.03 


4.326 ±0.051 


79 ±2 


4 113853 


6.03 ± 0.38 


3.50 ±0.12 


23150 ±761 


4.13 ±0.04 


3.499 ± 0.065 


134 ±5 


4 117831 


4.81 ±0.25 


2.71 ± 0.08 


19190 ±389 


4.26 ± 0.04 


2.949 ± 0.044 


118±4 


4 121084 


11.72 ±0.44 


3.64 ± 0.09 


31580 ±552 


4.39 ± 0.03 


4.072 ± 0.038 


223 ±6 


4 121110 


9.55 ± 0.33 


3.84 ± 0.07 


28390 ± 767 


4.25 ± 0.02 


3.933 ± 0.050 


175 ±3 


4 121461 


6.05 ± 0.26 


2.59 ±0.10 


22970 ± 783 


4.40 ± 0.04 


3.222 ± 0.068 


67 ±3 


4 159928 


11.27 ±0.46 


4.85 ±0.12 


29080 ± 870 


4.12 ±0.03 


4.178 ±0.056 


213 ±5 


4 160094 


9.05 ± 0.58 


3.23 ±0.15 


27550 ± 1106 


4.38 ± 0.05 


3.730 ±0.081 


96 ±4 


4 163552 


9.56 ±0.61 


5.06 ±0.16 


24400 ± 990 


4.01 ± 0.04 


3.910 ±0.075 


166 ±5 


4 175149 


11.78 ±0.47 


4.85 ±0.16 


33100 ±272 


4.14 ±0.03 


4.402 ± 0.032 


123 ±4 


4 175333 


5.09 ± 0.23 


2.41 ± 0.08 


21800 ±644 


4.38 ± 0.03 


3.071 ±0.059 


98 ±3 


5 016658 


6.38 ±0.39 


3.57 ±0.12 


21390 ±519 


4.14 ±0.04 


3.378 ±0.051 


145 ±5 


5 026631 


11.40 ±0.34 


5.19 ±0.08 


28670 ± 394 


4.07 ± 0.02 


4.212 ± 0.028 


186 ±3 


5 032412 


17.12 ±0.34 


4.74 ±0.11 


35410 ± 174 


4.32 ± 0.02 


4.500 ± 0.022 


66 ±2 


5 038089 


13.01 ± 0.20 


5.98 ± 0.08 


30660 ± 130 


4.00 ± 0.01 


4.452 ± 0.014 


127 ±2 


5 095337 


8.72 ± 0.36 


3.68 ±0.12 


26990 ± 1250 


4.25 ± 0.03 


3.808 ± 0.085 


206 ±7 


5 095557 


6.68 ± 0.71 


3.60 ±0.10 


23710 ±363 


4.15 ±0.05 


3.564 ± 0.036 


75 ±2 


5 100485 


6.67 ±0.14 


3.14 ±0.06 


22670 ± 498 


4.27 ± 0.02 


3.367 ± 0.042 


104 ±2 


5 100731 


8.73 ±0.51 


3.76 ±0.18 


24550 ± 473 


4.23 ± 0.05 


3.664 ± 0.054 


168 ±8 


5 106039 


8.69 ± 0.37 


3.63 ± 0.08 


26320 ± 335 


4.26 ± 0.03 


3.754 ± 0.029 


84 ±2 


5 111649 


5.56 ±0.16 


5.30±0.11 


17000 ± 168 


3.74 ± 0.02 


3.322 ± 0.024 


91 ±2 


5 123390 


9.83 ± 0.26 


4.38 ±0.14 


27840 ± 549 


4.15 ±0.03 


4.013 ±0.044 


102 ±3 


5 180185 


6.28 ±0.14 


4.41 ±0.14 


18500 ±400 


3.95 ± 0.03 


3.311 ±0.046 


41 ± 1 


5 180576 


7.30 ±0.31 


3.25 ±0.14 


25000 ± 1057 


4.28 ± 0.04 


3.568 ± 0.083 


105 ±5 


5 185408 


7.03 ± 0.21 


3.05 ± 0.09 


24170 ±475 


4.32 ± 0.03 


3.453 ± 0.042 


106 ±3 


5 196565 


7.20 ± 0.28 


4.93 ±0.13 


21530 ±368 


3.91 ±0.03 


3.670 ±0.038 


63 ±2 


5 261267 


9.17 ±0.39 


4.32 ±0.10 


27930 ± 904 


4.13 ±0.03 


4.008 ± 0.060 


171 ±4 


5 265970 


10.21 ±0.39 


7.14 ±0.16 


23280 ± 153 


3.74 ± 0.03 


4.128 ±0.023 


103 ±2 


5 266015 


15.58 ±0.31 


5.84 ± 0.08 


32120 ±545 


4.10 ±0.02 


4.512 ±0.032 


163 ±2 


5 266131 


8.96 ± 0.23 


3.60 ±0.13 


26510 ±774 


4.28 ± 0.03 


3.760 ± 0.060 


140 ±5 


5 266513 


5.79 ± 0.35 


2.84 ±0.13 


20800 ± 1167 


4.30 ± 0.05 


3.131 ±0.106 


130 ±6 


5 277080 


9.68 ± 0.55 


4.72 ± 0.09 


29780 ± 480 


4.08 ± 0.03 


4.196 ±0.032 


123 ±2 


5 283079 


7.56 ±0.41 


2.98 ± 0.08 


25160 ±475 


4.37 ± 0.03 


3.504 ±0.041 


118±3 



Computed from g = ^ 
Computed from log ^ 



21ogf +41og^ 



48 



P. North et al.: VLT spectroscopy of eclipsing binaries in the SMC 



Table 14. Astrophysical parameters for the secondary components: masses, radii, effective temperatures, surface gravities, lumi- 
nosities and equilibrium equatorial velocity (i.e. synchronous velocity for circular systems and pseudo-synchronous velocity for 
eccentric systems). 



Object 


Ms 


^s 


T^S a 
■'cff 


logos' 


logLs/io' 






(Mo) 


(Re) 


(K) 


(dex cgs) 




(kms"') 


4 110409 


7.30 ±0.36 


7.81 ±0.18 


21170 ±582 


3.52 ± 0.03 


4.041 ± 0.052 


133 ±3 


4 113853 


4.10 ±0.28 


3.79 ±0.13 


17060 ± 664 


3.89 ± 0.04 


3.038 ± 0.074 


145 ±5 


4 117831 


4.72 ± 0.25 


2.77 ± 0.08 


18830 ±588 


4.23 ± 0.04 


2.936 ± 0.060 


120 ±4 


4 121084 


9.98 ± 0.39 


3.28 ± 0.09 


29680 ± 560 


4.41 ± 0.03 


3.875 ± 0.041 


202 ±6 


4 121110 


7.13 ±0.26 


2.85 ± 0.07 


25010 ± 703 


4.38 ± 0.03 


3.455 ± 0.053 


130 ±3 


4 121461 


5.82 ± 0.26 


2.59 ±0.10 


21790 ± 909 


4.38 ± 0.04 


3.131 ±0.080 


67 ±3 


4 159928 


6.27 ± 0.30 


4.07 ±0.12 


21030 ± 686 


4.02 ± 0.03 


3.462 ± 0.062 


179 ±5 


4 160094 


6.81 ± 0.45 


2.63 ±0.14 


23500 ± 1059 


4.43 ± 0.05 


3.276 ±0.091 


78 ±4 


4 163552 


9.10 ±0.57 


4.54 ±0.15 


24330 ± 992 


4.08 ± 0.04 


3.812 ±0.077 


149 ±5 


4 175149 


7.76 ±0.31 


6.18 ±0.20 


28110 ±306 


3.75 ± 0.03 


4.330 ± 0.034 


156 ±5 


4 175333 


3.98 ±0.14 


2.06 ± 0.07 


19170 ±741 


4.41 ± 0.03 


2.713 ± 0.074 


84 ±3 


5 016658 


5.67 ±0.35 


2.84 ±0.11 


21510 ±553 


4.29 ± 0.04 


3.190 ±0.057 


115±5 


5 026631 


7.86 ± 0.25 


4.93 ± 0.08 


21050 ±355 


3.95 ± 0.02 


3.631 ±0.033 


177 ±3 


5 032412 


13.14 ±0.28 


4.06 ±0.11 


31160±315 


4.34 ± 0.02 


4.144 ±0.029 


57 ±2 


5 038089 


11.70±0.18 


4.94 ± 0.07 


31020 ±301 


4.12 ±0.01 


4.307 ± 0.021 


105 ±2 


5 095337 


7.62 ± 0.32 


3.23 ±0.12 


25270 ± 1185 


4.30 ± 0.04 


3.581 ±0.087 


181 ±6 


5 095557 


5.70 ± 0.64 


2.86 ± 0.08 


21530 ±446 


4.28 ± 0.05 


3.198 ±0.044 


60 ±2 


5 100485 


6.49 ±0.13 


3.06 ± 0.06 


22840 ± 527 


4.28 ± 0.02 


3.359 ± 0.044 


102 ±2 


5 100731 


6.65 ± 0.40 


3.23 ±0.16 


20300 ±931 


4.24 ± 0.05 


3.201 ± 0.091 


144 ±7 


5 106039 


4.58 ± 0.22 


5.46 ±0.11 


16800 ± 265 


3.63 ± 0.03 


3.328 ± 0.032 


126 ±2 


5 111649 


5.21 ±0.15 


4.89 ±0.10 


17420 ± 260 


3.78 ± 0.02 


3.295 ± 0.032 


84 ±2 


5 123390 


7.91 ± 0.22 


3.31 ±0.13 


28320 ± 784 


4.30 ± 0.04 


3.800 ± 0.059 


77 ±3 


5 180185 


5.66 ±0.13 


3.60 ±0.13 


19170 ±657 


4.08 ± 0.03 


3.197 ±0.067 


33 ± 1 


5 180576 


5.73 ± 0.25 


2.65 ±0.13 


20310 ± 1090 


4.35 ± 0.05 


3.030 ±0.103 


86 ±4 


5 185408 


5.97 ±0.19 


2.80 ± 0.08 


22830 ± 640 


4.32 ± 0.03 


3.280 ± 0.055 


97 ±3 


5 196565 


5.78 ± 0.23 


3.69 ±0.13 


18670 ±461 


4.07 ± 0.03 


3.170 ±0.052 


47 ±2 


5 261267 


4.81 ±0.23 


3.86 ±0.10 


20330 ± 666 


3.95 ± 0.03 


3.359 ± 0.061 


153 ±4 


5 265970 


6.33 ± 0.30 


3.72 ±0.15 


20980 ± 340 


4.10 ±0.04 


3.381 ± 0.044 


54 ±2 


5 266015 


6.86 ±0.16 


5.46 ± 0.08 


24320 ± 438 


3.80 ± 0.02 


3.971 ± 0.034 


153 ±2 


5 266131 


7.74 ± 0.21 


3.10±0.12 


24820 ± 750 


4.35 ± 0.04 


3.513 ±0.063 


120 ±5 


5 266513 


4.90 ± 0.30 


2.45 ±0.13 


19940 ± 1184 


4.35 ± 0.05 


2.931 ±0.113 


112±6 


5 277080 


4.88 ± 0.29 


5.13 ±0.09 


20410 ± 372 


3.71 ±0.03 


3.612 ±0.035 


134 ±2 


5 283079 


7.58 ± 0.41 


2.97 ± 0.08 


25160 ±534 


4.37 ± 0.03 


3.501 ± 0.044 


117±3 



^ Computed from T'^ 
^ Computed from g 



' '-'cff' -'cff)pli')l-'eff 



GM 



Computed from log 4^ =2 log -^ + 4 log r^r^^ 



eif,o 
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Table 15. Distance determination: bolometric absolute magnitudes, bolometric corrections, visual absolute magnitudes, visual 
extinction and reddening-free distance modulus in the V- and /-band. The /-band distance modulus is available only for systems 
for which the B and V magnitudes are available. The "q" superscript denotes values at quadrature. 



Object 



"bol 
(mag) 



M? ? 

bol 

(mag) 



(mag) 



BC^V 
(mag) 



(mag) 



(mag) 



(mag) 



(mag) 



51og[d]-5j; 
(mag) 



5 log[d] - 5/ 
(mag) 



4 110409 


-6.06 ±0.13 


-5.35 + 0.13 


-3.10 + 0.06 


-2.12 + 0.07 


-2.96 ± 0.08 


-3.23 + 0.07 


0.207 ± 0.009 


0.642 + 0.029 


19.06 + 0.061 


19.06 + 0.057 


4 113853 


-4.00 + 0.16 


-2.85 + 0.18 


-2.36 + 0.08 


-1.60 ±0.09 


-1.64±0.10 


-1.25 + 0.10 


0.194 ±0.009 


0.601 + 0.029 


19.00 + 0.078 


18.99 + 0.075 


4 117831 


-2.62 ±0.11 


-2.59 + 0.15 


-1.90 + 0.05 


-1.85 ±0.07 


-0.73 ± 0.08 


-0.74 + 0.09 


0.064 ± 0.006 


0.198 + 0.018 


18.99 + 0.062 


18.97+0.061 


4 121084 


-5.43 ±0.09 


-4.94 + 0.10 


-3.06 + 0.04 


-2.93 ± 0.05 


-2.37 ± 0.07 


-2.00 + 0.07 


0.187 ±0.005 


0.580 + 0.017 


19.28 + 0.051 


19.27+0.049 


4 121110 


-5.08 ±0.12 


-3.89 + 0.13 


-2.83 + 0.06 


-2.55 ±0.07 


-2.25 ± 0.07 


-1.34 + 0.08 


0.158 ±0.008 


0.490 + 0.023 


19.08 + 0.057 


19.08+0.054 


4 121461 


-3.30 ±0.17 


-3.08 + 0.20 


-2.34 + 0.08 


-2.22 ±0.10 


-0.96 ±0.11 


-0.86 + 0.12 


0.147 ±0.009 


0.456 + 0.027 


19.05 + 0.083 


19.09 + 0.080 


4 159928 


-5.69 ±0.14 


-3.90 + 0.15 


-2.88 + 0.07 


-2.12 ±0.08 


-2.81 ±0.08 


-1.78 + 0.09 


0.168 ±0.007 


0.521+0.021 


19.29 + 0.066 


19.30 + 0.064 


4 160094 


-4.57 ±0.20 


-3.44 + 0.23 


-2.77 + 0.10 


-2.40 ±0.11 


-1.80±0.13 


-1.04 + 0.14 


0.075 ± 0.009 


0.233 + 0.027 


18.96 + 0.102 


18.97+0.100 


4 163552 


-5.02 ±0.19 


-4.78 + 0.19 


-2.48 + 0.10 


-2.47 ±0.10 


-2.55 ±0.10 


-2.30 + 0.11 


0.187 ±0.008 


0.580 + 0.026 


18.35 + 0.079 


18.36 + 0.077 


4 175149 


-6.26 ±0.08 


-6.07 + 0.09 


-3.15+0.02 


-2.79 ± 0.03 


-3.11 ±0.07 


-3.29 + 0.07 


0.053 ± 0.007 


0.164 + 0.022 


18.52 + 0.057 


18.54 + 0.054 


4 175333 


-2.93 ±0.15 


-2.03 + 0.19 


-2.22 + 0.07 


-1.89 ±0.09 


-0.71 ±0.09 


-0.14 + 0.11 


0.070 ± 0.007 


0.217 + 0.022 


18.61+0.074 


18.60 + 0.072 


5 016658 


-3.70 ±0.13 


-3.23 + 0.14 


-2.17 + 0.06 


-2.19 ±0.06 


-1.53 ±0.08 


-1.04 + 0.10 


0.080 ± 0.006 


0.248 + 0.019 


19.13 + 0.068 


19.15+0.066 


5 026631 


-5.78 ±0.07 


-4.33 + 0.08 


-2.85+0.03 


-2.12 ±0.04 


-2.93 ±0.04 


-2.20 + 0.05 


0.105 ±0.004 


0.326 + 0.013 


19.13 + 0.036 


19.17+0.034 


5 032412 


-6.50 ± 0.05 


-5.61+0.07 


-3.32 + 0.01 


-3.03 ±0.02 


-3.18 ±0.05 


-2.58 + 0.06 


0.250 ± 0.006 


0.775 + 0.019 


19.19 + 0.044 


19.22 + 0.041 


5 038089 


-6.38 ±0.04 


-6.02 + 0.05 


-2.98 + 0.01 


-3.01 ±0.02 


-3.40 ±0.03 


-3.01+0.04 


0.057 ±0.011 


0.177 + 0.034 


18.89 + 0.042 


18.83+0.032 


5 095337 


-4.77 ±0.21 


-4.20 + 0.22 


-2.72 + 0.11 


-2.57 ±0.11 


-2.05 ±0.11 


-1.64 + 0.12 


0.124 ±0.016 


0.384 + 0.050 


19.17 + 0.097 


19.21+0.089 


5 095557 


-4.16 ±0.09 


-3.24 + 0.11 


-2.42 + 0.04 


-2.19 ±0.05 


-1.74 ±0.07 


-1.06 + 0.07 


0.1 16 ±0.004 


0.360 + 0.013 


19.19 + 0.052 


19.20 + 0.051 


5 100485 


-3.67 ±0.10 


-3.65 + 0.11 


-2.31+0.05 


-2.32 ± 0.06 


-1.36 ±0.06 


-1.32 + 0.06 


0.077 ± 0.009 


0.239 + 0.028 


18.84 + 0.052 


18.87+0.047 


5 100731 


-441 ±0.13 


-3.25 + 0.23 


-2.50 + 0.05 


-2.04 ±0.11 


-1.91 ±0.11 


-1.21 +0.14 


0.106 ±0.005 


0.329 + 0.016 


19.28 + 0.090 


19.32+0.089 


5 106039 


-4.63 ± 0.07 


-3.57 + 0.08 


-2.66 + 0.03 


-1.56 ±0.04 


-1.97 ±0.05 


-2.01 + 0.05 


0.154 ±0.010 


0.477 + 0.032 


18.95 + 0.050 


18.95+0.042 


5 111649 


-3.56 ±0.06 


-3.49 + 0.08 


-1.59 + 0.02 


-1.65 ±0.04 


-1.97 ±0.05 


-1.84 + 0.05 


0.177 ±0.009 


0.549 + 0.028 


18.86 + 0.046 


18.88+0.040 


5 123390 


-5.28 ±0.11 


-4.75 + 0.15 


-2.78 + 0.05 


-2.83 ± 0.07 


-2.50 ±0.08 


-1.92 + 0.10 


0.092 ±0.015 


0.285 + 0.048 


18.75 + 0.078 


18.76 + 0.068 


5 180185 


-3.53 ±0.12 


-3.24 + 0.17 


-1.80 + 0.05 


-1.89 ±0.08 


-1.73 ±0.08 


-1.35 + 0.10 


0.067 ±0.012 


0.208 + 0.036 


19.45 + 0.073 


19.31+0.067 


5 180576 


-4.17 ±0.21 


-2.82 + 0.26 


-2.54 + 0.10 


-2.04 ±0.13 


-1.63 ±0.13 


-0.79 + 0.15 


0.150 ±0.013 


0.465 + 0.039 


19.14 + 0.106 


19.13+0.101 


5 185408 


-3.88 ±0.11 


-3.45 + 0.14 


-2.47 + 0.05 


-2.32 ± 0.07 


-1.41 ±0.07 


-1.13 + 0.09 


0.100 ±0.005 


0.310 + 0.015 


19.12 + 0.057 


19.13+0.056 


5196565 8 


-4.42 ± 0.09 


-3.17 + 0.13 


-2.18 + 0.04 


-1.82 ±0.06 


-2.25 ± 0.07 


-1.35 + 0.09 










5 261267 


-5.27 ±0.15 


-3.65 + 0.15 


-2.79 + 0.08 


-2.04 ± 0.08 


-2.48 ± 0.08 


-1.61 +0.08 


0.085 ± 0.009 


0.263 + 0.027 


19.35 + 0.068 


19.29 + 0.064 


5 265970 


-5.57 ±0.06 


-3.70 + 0.11 


-2.36 + 0.02 


-2.12 ±0.04 


-3.21 ±0.05 


-1.58 + 0.09 


0.078 ± 0.005 


0.242 + 0.016 


19.25 + 0.048 


19.28+0.046 


5 266015 


-6.53 ±0.08 


-5.18 + 0.08 


-3.08 + 0.04 


-2.46 ± 0.04 


-3.45 ±0.05 


-2.71+0.05 


0.179 ±0.005 


0.555 + 0.017 


19.23 + 0.038 


19.25+0.036 


5 266131 


-4.65 ±0.15 


-4.03 + 0.16 


-2.68 + 0.07 


-2.53 ±0.07 


-1.97 ±0.10 


-1.50 + 0.11 


0.136 ±0.012 


0.422 + 0.036 


19.11 +0.081 


19.12 + 0.075 


5 266513 


-3.08 ±0.26 


-2.58 + 0.28 


-2.10 + 0.13 


-1.99±0.14 


-0.98 ±0.15 


-0.58 + 0.16 


0.159 ±0.012 


0.493 + 0.038 


19.13 + 0.118 


19.11 +0.114 


5 277080 


-5.74 ±0.08 


-4.28 + 0.09 


-2.93 + 0.04 


-2.04 ± 0.04 


-2.81 ±0.05 


-2.24 + 0.05 


0.225 ±0.013 


0.698 + 0.042 


18.52 + 0.056 


18.71 +0.046 


5 283079 


-4.01 ±0.10 


-4.00 + 0.11 


-2.56 + 0.05 


-2.56 ±0.05 


-1.45 ±0.07 


-1.44 + 0.07 


0.161 ±0.005 


0.499 + 0.017 


19.11+0.054 


19.07+0.052 



Computed from A/^oi = -5 log jj^ — 10 log ^ 



^'^'bol 



Interpolated from Lanz & Hubeny I2003I2007I 
Computed from My = Mj,yi - BCy 
Computed from Eg^y = (S - V) - (B - V)o 
Computed from Ay =KvEb-v wlthRv/ = 3.1 ±0.3 

Computed froin 5 log[rf] - 5 = C'l - mJ, - A J where mJ, = -2.5 log! 10"' 

Computed froin 5 log[rf] - 5 = /'I - AfJ + (V - 7)2 - 0.600 AJ 
B and V light curve.s missing 



Table 16. Light curves: ratio of the primary minimum to the RMS (/-band), RMS scatters and minimum chi-squared values 
from WD/PHOEBE code. 



Object 


A/minl/'^; 


o-l 


o-y 


D"8 


Object 


A'minl/o"; 


a-l 


o-y 


o"i 






(mag) 


(mag) 


(mag) 






(mag) 


(mag) 


(mag) 


4 1 10409 


-65 


0.016 


0.028 


0.042 


5 100485 


30.1 


0.017 


0.018 


0.020 


4 113853 


9.3 


0.017 


0.019 


0.025 


5 100731 


19.2 


0.021 


0.013 


0.022 


4 117831 


10.8 


0.028 


0.018 


0.025 


5 106039 


23.7 


0.016 


0.107 


0.114 


4 121084 


34.6 


0.019 


0.013 


0.016 


5 111649 


19.7 


O.OII 


0.017 


0.029 


4 121110 


31.7 


0.015 


0.016 


0.024 


5 123390 


13.8 


0.010 


0.050 


0.063 


4 121461 


10.5 


0.031 


0.025 


0.017 


5 180185 


18.6 


0.019 


0.040 


0.038 


4 159928 


18.7 


0.013 


0.008 


0.009 


5 180576 


9.7 


0.021 


0.021 


0.030 


4 160094 


12.5 


0.019 


0.008 


0.011 


5 185408 


10.0 


0.023 


0.012 


0.013 


4 163552 


46.1 


0.009 


0.0 1 1 


0.017 


5196565" 


21.0 


0.017 






4 175149 


36.1 


0.024 


0.021 


0.028 


5 261267 


49.7 


0.017 


0.052 


0.024 


4 175333 


12.1 


0.027 


0.020 


0.022 


5 265970 


334 


0.012 


0.015 


0.016 


5 016658 


224 


0.020 


0.009 


0.020 


5 266015 


49.2 


0.012 


0.018 


0.015 


5 026631 


24.5 


0.010 


0.011 


0.013 


5 266131 


7.8 


0.019 


0.026 


0.101 


5 032412 


334 


0.009 


0.012 


0.028 


5 266513 


15.5 


0.034 


0.020 


0.025 


5 038089 


22.4 


0.014 


0.042 


0.049 


5 277080 


43.1 


0.013 


0.027 


0.100 


5 095337 


17.9 


0.018 


0.186 


0.206 


5 283079 


26.5 


0.020 


0.017 


0.024 


5 095557 


21.1 


0.026 


0.022 


0.019 













3 and V light curves missing 
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Table 17. Radial velocity curves: RMS scatters. 



Objecl 


(Tp 


•^s 


Objecl 


o-p 


o-s 




(kms-') 


(km S-' ) 




(kms"') 


(kms"') 


4 110409 


11.5 


4.2 


5 100485 


6.1 


4.4 


4 113853 


7.8 


12.4 


5 100731 


7.0 


17.4 


4 117831 


11.0 


12.6 


5 106039 


8.1 


12.2 


4 121084 


7.6 


13.8 


5 111649 


4.3 


5.5 


4 121110 


7.6 


9.1 


5 123390 


8.4 


14.1 


4 121461 


7.6 


10.1 


5 180185 


4.0 


2.9 


4 159928 


8.3 


11.9 


5 180576 


8.0 


12.8 


4 160094 


13.7 


29.4 


5 185408 


7.1 


12.7 


4 163552 


10.5 


8.8 


5 196565 


5.7 


9.5 


4 175149 


8.6 


10.5 


5 261267 


7.3 


13.5 


4 175333 


19.0 


23.4 


5 265970 


12.8 


11.8 


5 016658 


7.7 


17.2 


5 266015 


5.2 


9.8 


5 026631 


7.8 


7.4 


5 266131 


9.9 


8.4 


5 032412 


3.3 


4.9 


5 266513 


14.6 


9.8 


5 038089 


3.1 


2.9 


5 277080 


8.3 


8.0 


5 095337 


8.3 


12.0 


5 283079 


124 


11.7 


5 095557 


31.0 


20.6 










16658 



0.5 
phase 



Fig. 30. /-band light curves with O - C residuals. See text for comments on individual stars. 
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Fig. 31. Same as Fig.|30] for 12 more stars. 
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Fig. 32. Same as Fig.|30] for nine more stars. 
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Fig. 33. Radial velocity data and best-fitting RV curves. The O - C residuals are shown with a different arbitrary off'set for each 
component. See text for comments on individual stars. 
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Fig. 34. Same as Fig.|33] for 12 more stars. 
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Fig. 35. Same as Fig.|33] for nine more stars. 
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Fig. 36. Mass-surface gravity diagrams: the positions of the primary and the secondary components are indicated by filled and 
open symboles, respectively. The dotted lines ares isochrones from Charbonnel et al. (I1993I I at Z = 0.004, with ages of 3, 5, 10, 
20, 30, 40, 50 and 100 Myr. See text for comments on individual stars. 
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Fig. 37. Same as Fig.|36] for 12 more stars. 
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Fig. 38. Same as Fig.|36] for nine more stars. 
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Fig. 39. HR diagrams: the positions of the primary and the secondary components are indicated by filled and open symbols, 
respectively. The primary is the component with the larger mass. The evolutionary tracks from Charbonnel et al. ( I1993I I at 
Z = 0.004, corresponding to the observed masses, are indicated by solid (primary) and dashed black lines (secondary). Since 
these authors adopt a helium content Y - 0.24 + 3 x AZ, the helium content of these models is F = 0.252. The oblique line 
corresponds to the ZAMS. The bold segments on the ZAMS, at the departure point of the evolutionary tracks, indicate the +1 cr 
error on the mass. The red tracks and ZAMS correspond to a poorer metallicity Z - 0.001 (Schaller et al. |1992| l and a helium 
content Y = 0.243. The green tracks are interpolated from the Z - 0.004 models of Claret & Gimenez ( 119981 ) for a helium content 
Y = 0.28; thus, they show the effect of a helium enhancement AY - 0.028 relative to the tracks of Charbonnel et al. ( I1993I I. See 
text for comments on individual stars. 
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Fig. 40. Same as Fig.|39] for 12 more stars. 
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Fig. 41. Same as Fig.|39] for nine more stars. 



